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Summary 



The present thesis, which is mainly based on my research papers I-IV hsted 
overleaf, is devoted to some forefront issues in the field of leptogenesis. As 
well as nicely explaining the origin of the matter-antimatter asymmetry of 
the Universe, this mechanism is intimately related to the nature of neutrino 
masses in that it is the cosmological consequence of the see-saw mechanism. 
This connection makes possible to relate parameters measured or to be mea- 
sured in neutrino experiments to a cosmological parameter, the baryon-to- 
photon ratio, r^^ — 6 x 10^^'^. 

According to the mass M of the heavy neutrino that decays, the gener- 
ation of asymmetry has to be described by different sets of classical Boltz- 
mann equations: in the unflavored regime, when M > 10^^ GeV, the lepton 
flavors are summed over, and the asymmetry is estimated from two equa- 
tions, one of which tracks Nb-l, the total B~L number, where B and L 
denote baryon and lepton numbers, respectively. In the fully flavored regime, 
when M < 10^^ GeV, the lepton flavor asymmetries B/3 — Lq,, a — e + /j, or 
T, must be tracked instead. The fact that two descriptions coexist is due to 
the r-lepton Yukawa interactions, which become important roughly when the 
temperature in the primordial plasma drops below 10^^ GeV. As we pointed 
out in paper III, the transition between the two regimes requires a quan- 
tum kinetic treatment, where correlations in flavor space and partial losses 
of coherence are properly accounted for. 

In the unflavored regime, when the heavy neutrino masses are hierarchi- 
cal, one typically obtains a scenario where only the contribution from the 
lightest right-handed neutrino, A^i, matters. A lower bound on Mi of about 
10^ GeV for successful leptogenesis can then be derived. The associated lower 
bound on the reheat temperature T^eh of the Universe after inflation being of 
the same order, a tension arises between leptogenesis and supergravity the- 
ories because of a possible overproduction of gravitinos. Furthermore, it is 
interesting that one finds a stringent upper bound on the absolute neutrino 
mass scale, mi < 0.1 eV. Finally, it should be noted that, in the unflavored 
regime, the final asymmetry does not depend on the lepton mixing matrix. 
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which contains the potentially measurable Dirac and Majorana CP-violating 
phases. The necessary source of CP violation is instead provided by unmea- 
surable phases in the high-energy sector. 

In the fully flavored regime, the lowest bounds on Mi and on the reheat 
temperature do not change, as we stressed in paper II. On the other hand, 
the parameters accessible in neutrino experiments, notably the angle ^13 and 
the Dirac and Majorana CP-violating phases, acquire importance. According 
to the value they take, the predictions from leptogenesis can differ by orders 
of magnitude, as we showed in paper II. Moreover, as investigated in detail 
in paper IV, the source of CP violation required for leptogenesis may be 
exclusively provided by the Dirac phase, which is the low-energy phase with 
the best prospects of being measured in the near future. 

There are different ways to go beyond the "vanilla" picture where the 
asymmetry is produced by the lightest right-handed neutrino and the lower 
bounds introduced above hold. Following paper I, three possibilities are dis- 
cussed. First, a quasi-degenerate mass spectrum for the heavy neutrinos is 
considered, implying that two or three heavy neutrinos contribute to the 
generation of asymmetry. It is then possible to relax the bounds on Mi and 
Treh to the TeV scale, and the upper bound on mi can be evaded. Second, 
specific situations can be found where the asymmetry is produced by the 
next-to-lightest right-handed neutrino, N2. The lower bound on Mi is then 
replaced by a lower bound on M2, still implying a lower bound on Treh similar 
to the vanilla case. Finally, even within the A'^i-dominated scenario, a special 
choice of high-energy parameters, which requires large cancellations, can be 
chosen in order to relax the lower bounds on Mi and T^eh- 

Leptogenesis has recently entered a new era in that some of what has been 
discussed above has only been realized in the last couple of years. Notably, 
let me mention the importance of flavor effects, the possible relevance of 
quantum effects, as well as the role played by the next-to-lightest right- 
handed neutrino. 

In the future, the case for leptogenesis will be weakened or strengthened 
depending on the outcome of the next-generation neutrinoless double-beta 
decay searches, the planned long-baseline neutrino experiments, as well as 
the LHC. 



ii 



The content of this thesis is mainly based on the following works: 

I S. Blanchet and P. Di Bari, 

Leptogenesis beyond the limit of hierarchical heavy neutrino masses, 
JCAP 0606 (2006) 023 [arXiv:hep-ph/0603107]. 

II S. Blanchet and P. Di Bari, 

Flavor effects on leptogenesis predictions, 
JCAP 0703 (2007) 018 [arXiv:hep-ph/0607330]. 

III S. Blanchet, P. Di Bari and G. G. Raffelt, 

Quantum Zeno effect and the impact of flavor in leptogenesis, 
JCAP 0703 (2007) 012 [arXiv:hep-ph/0611337]. 

IV A. Anisimov, S. Blanchet and P. Di Bari, 
Viability of Dirac phase leptogenesis, 
JCAP 0804 (2008) 033 [arXiv:0707.3024]. 



iii 



iv 



Acknowledgments 



It is a very difficult task, if not impossible, to list all people who contributed 
directly or indirectly to the achievement of my Ph.D. studies. Therefore, I 
shall undertake the more modest task of underlining the major contributions. 

I would like first to express my gratitude to the Max-Planck-Institut 
fiir Physik (MPl) for providing funding and optimal research possibili- 
ties through the newly born International Max Planck Research School 
for Elementary Particle Physics. Thanks also to the Technische Universitat 
Miinchen (TUM), in particular to Prof. Michael Ratz, for providing the aca- 
demic framework to my Ph.D. studies. 

Next, I would like to thank my supervisor within the institute, Gcorg 
Raffelt, for his kindness, support, availability, and the many useful advice he 
gave me both at the scientific and non-scientific levels. 

For the research presented in this thesis, I am mostly indebted to Pasquale 
Di Bari. I benefited a lot from his expertise in the fields of leptogencsis and 
cosmology in general. Moreover, I enjoyed the long discussions about our 
work and beyond. 

Then, my thanks go to all members of the "astroparticle" theory group 
who overlapped with me since October 2005 at the MPI for a nice, relaxed 
and inspiring atmosphere. Special thanks to Yvonne Wong for suggestions 
about English issues and to Pasquale Serpico for useful information about 
the process of terminating a Ph.D. at the MPI and TUM. Thanks also to 
my office mate Stefano Pozzorini for the enjoyable "Swiss" atmosphere inside 
the office. Finally, thanks to the other IMPRS students, especially Florian 
Hahn-Wornle, Edoardo Mirabella, and Joseph Pradler. 

Last but not least, I am grateful to my parents for the irreplaceable role 
they played from the very beginning until the present achievement of my 
studies. 



V 



vi 



Contents 



Summary i 

Acknowledgments v 

1 Introduction 1 

1.1 The matter- antimatter puzzle 1 

1.2 The puzzle of neutrino masses 7 

1.2.1 Theory of neutrino oscillations 7 

1.2.2 Experimental evidence for neutrino oscillations 8 

1.2.3 Absolute neutrino mass scale 11 

1.3 The see-saw mechanism and leptogenesis 12 

2 Vanilla leptogenesis 17 

2.1 General scenarios of unfiavored leptogenesis 17 

2.2 The A^i-dominated scenario 24 

2.3 Leptogenesis and supersymmetry 35 

3 Adding flavor to vanilla leptogenesis 39 

3.1 When does flavor matter and why? 39 

3.2 Flavored Boltzmann equations and spectator processes .... 41 

3.3 In practice, what changes? 44 

3.4 Dependence on the initial conditions and lower bounds .... 48 

3.4.1 Alignment 49 

3.4.2 Democratic and semi-democratic cases 51 

3.4.3 One-flavor dominance 52 

3.5 Study of a specific example 54 

4 From classical to quantum kinetic equations 61 

4.1 Validity of the different pictures 62 

4.1.1 When are flavor effects important? 62 

4.1.2 Maximum ffavor effects 65 



vii 



4.2 Neutrino mass bound 67 

4.3 Limitations of a simple rate comparison 70 

4.4 Density matrix equation 71 

5 Going beyond vanilla leptogenesis 79 

5.1 Degenerate limit for the heavy neutrinos 80 

5.2 A^2-dominated scenario 83 

5.3 Effects of 11^221 85 

6 Leptogenesis from low-energy CP-violating phcises 87 

6.1 CP violation in neutrino physics 88 

6.1.1 Neutrino oscillations and the Dirac phase 88 

6.1.2 Neutrinoless double-beta decay and the Majorana phases 89 

6.2 Dirac phase leptogenesis 90 

6.2.1 The hierarchical limit 90 

6.2.2 The degenerate limit 101 

6.3 Leptogenesis from the Majorana phases 110 

6.4 Discussion Ill 

7 Conclusion 115 
A Neutrino mixing parameters 121 
B The see-saw mechanism with three RH neutrinos 123 



viii 



Chapter 1 
Introduction 



1.1 The matter-antimatter puzzle 

One can surely say that our understanding of the Universe has made a huge 
leap forward in the last few years. This is partly because the amount of data 
has dramatically increased, so that an epoch of "precision cosmology" has 
started. But this is also due to an accumulation of evidence for concepts that 
were still considered exotic not so long ago, such as dark matter and dark 
energy. Even though their nature is still unknown, at least there seems to be 
a consensus about their existence. 

Nowadays, one speaks about a "Standard Cosmological Model" , in anal- 
ogy with its very successful counterpart of particle physics. The Standard 
Cosmological Model tells us that the Universe is in a phase of accelerated 
expansion and that the total energy in the Universe is shared among at least 
four components (see Fig. 1.1) which sum to Qtot — 1, meaning that the 
Universe is flat to a good precision. The dominant component (about 73%) 
is called dark energy, dark matter makes about 23%, ordinary matter (both 
luminous and dark) only 4% and neutrinos 0.2-2%, the uncertainty here 
stemming from the unknown absolute neutrino mass scale, as we shall see in 
Section 1.2.3. 

It is well known that the nature of dark matter is still mysterious. The 
particle interpretation seems to be widely supported, and the candidates 
are numerous: axion, lightest supersymmetric particle (neutralino, gravitino), 
sterile neutrino, and many others. 

Dark energy is probably a much bigger issue still than dark matter. It is 
supposed to drive the accelerated expansion, but its nature is very unclear. 
It also raises a fundamental question about how to treat the vacuum energy 
in quantum fleld theory. 
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Figure 1.1: The mass-energy bud- 
get of the Universe. 



Figure 1.2: The antiproton-to- 
proton ratio at the top of the 
atmosphere, as observed (points) 
and predicted from the models 
(lines) [1]. 



Ordinary matter, which constitutes our bodies as well as the Earth and 
the stars, does not seem at first to introduce any challenge to our under- 
standing. However, this naive perception is wrong because two very puzzling 
questions remain: 

1) Why is antimatter essentially absent in the observable Universe? 

2) Why is the number density of baryons so small compared to photons 
or neutrinos? 

These two questions are puzzling because, according to the Standard Big- 
Bang Theory, matter and antimatter evolved in the same way in the early 
Universe. On the other hand, today the observable Universe is composed 
almost exclusively of matter. Antimatter is only seen in particle physics ac- 
celerators and in cosmic rays. Moreover, the rates observed in cosmic rays 
are consistent with the secondary emission of antiprotons, rtp/np ~ 10~^ (see 
Fig. 1.2). 

Ordinary matter is composed of baryons (protons, neutrons) and leptons 
(electrons). One can assign an experimentally conserved number to baryons 
and leptons. Baryons and leptons carry one unit of these numbers, whereas 
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antibaryons and antileptons carry one negative unit. In this way one can 
say that the predominance of matter over antimatter is equivalent to the 
existence of a net baryon number. 

Following the Standard Big-Bang Theory and relying on the Standard 
Model (SM) of particle physics, the relic density of baryons, i.e. nucleons 
here, can be easily estimated. One has the usual Boltzmann equation for the 
number density of baryons ub (or antibaryons) 



dz 



+ 3HnB 



11% 



[n 



cqx2 
B) 



(1.1) 



eq 
B 



where z — Mb/T, Mb is the baryon mass, H the Hubble expansion rate, n 
the equilibrium number density of baryons, and T the cosmic temperature. 
The collision term on the right-hand side is given in terms of an thermally- 
averaged annihilation cross-section (cyilf |), whose definition can be found 
in [2] for example. Eq. (1.1) can be conveniently rewritten in terms of the 
variable n^/n^, where is the number density of photons, allowing one to 
factor out the effects of the expansion of the Universe. One obtains 



d{nB/n^) 
dz 



n^z 



H{Mb) 



riB 



eq 



(1.2) 



In our case, the important annihilation channel is into pions (e.g. 
p + p ^ 7r~^ + 7r^). Taking the averaged cross-section to be ((xlf |) = CiM~^, 
with — 135 MeV and where Ci is a numerical factor of order unity, the 
freeze-out occurs at T ~ 22 MeV. Neglecting the entropy production in 
e+e~-annihilations, one finds for today's abundance (subscript '0') [2] 



Ub 



7 X 10-^°Cf ^ 



(1.3) 



Note that the ratio of baryon number density to photon number density 
today is usually referred to as the haryon-to-photon ratio, 



Ub 

Vb= — 



(1.4) 



The result of our simple computation, Eq. (1.3), is clearly a small number, 
which would perhaps explain the question 2) above. But one notices imme- 
diately a first problem, namely because the abundances of baryons and an- 
tibaryons are predicted in this way to be the same. Baryons and antibaryons 
do not evolve in distinctive ways so that one expects today the same amount 
of each of them. So, the argument we have just described leaves open the 
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question 1) above. But let us for the moment ignore this point and try to see 
if the abundance of baryons matches observation. 

The photon density follows directly from the measurement of the Cosmic 
Microwave Background (CMB) temperature and from Bose- Einstein statis- 
tics: n-y ~ T"^. Determining the baryon content of the Universe is more diffi- 
cult. Direct measurements are not accurate, because only a small fraction of 
baryons formed stars and other luminous objects (see Fig. 1.1). However, we 
can rely on two different indirect probes. 

The first probe is Big-Bang Nucleosynthesis (BBN). The abundances of 
hght elements such as ^He, D, ^He and ^Li predicted by the standard theory 
of BBN crucially depend on rjB- Comparing predictions with observations, as 
shown in Fig. 1.3, the following baryon-to- photon ratio is inferred [3]: 

77b ~ (5.5 ± 1.0) X 10-^°. (1.5) 

The CMB temperature anisotropics, very well measured by the WMAP 
satellite, offer the second probe. These anisotropics reflect the acoustic os- 
cillations of the baryon-photon fluid which happened around photon decou- 
pling. A precise computation can be done evolving Boltzmann equations for 
anisotropics, assuming that they are generated by quantum fluctuations dur- 
ing inflation. Fig. 1.4 illustrates how the amount of anisotropics with angular 
scale ~ l/£ depends on rjB- The baryon-to-photon ratio obtained from 3 years 
of WMAP data is [4] 

~ (6.1 ±0.2) X 10"^°. (1.6) 

The synthesis of light elements occurred during the first 3 minutes in 
the history of the Universe, whereas the photon decoupling occurred when 
the Universe was 400 thousand years old. The fact that these two completely 
different probes of the baryon content of the Universe give compatible results 
is one of the great successes of modern cosmology.^ 

It should now be clear that the result from the "classical" computation 
of the baryon density given in Eq. (1.3) is at odds with the observed value, 
Eq. (1.6). In order to avoid the "baryon annihilation catastrophe" leading to 
the value in Eq. (1.3), one has to generate a primordial asymmetry between 
baryons and antibaryons. This small asymmetry, at the level of one part in 
one billion, would imply that after the annihilation process has occurred at 
full strength, one remains with the small excess of baryons over antibaryons. 
The problem of generating this small excess of baryons over antibaryons is 
often called the baryogenesis problem. 

^Throughout the thesis, we shah exlusively use the value of r/s obtained from CMB 
temperature anisotropies, Eq. (1.6), which has much smaller errors. 
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Figure 1.3: The observed abundances of light elements compared to the stan- 
dard BBN predictions [5]. The smaller boxes indicate 2a statistical errors, 
the larger ones ±2cr statistical and systematic errors. 




Figure 1.4: The dependence of CMB temperature anisotropies on the baryon 
abundance fib (or r^^), compared with data [6]. 
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The solution to the baryogenesis problem requires the generation of a 
small baryon asymmetry primordially. Sakharov, in 1967, enunciated the 
three necessary conditions for such a process to be possible at some stage 
in the history of the Universe [7] : 

1. Baryon number violation. 

2. C and CP violation. 

3. Departure from thermal equilibrium. 

In principle, the SM contains all these ingredients. Indeed, 

1. Due to the chiral nature of weak interactions, B and L are not con- 
served [8] . At zero temperature, this has no observable effect due to the 
smallness of the weak coupling. However, as the temperature reaches 
the critical temperature Tew of the electroweak phase transition, B and 
L violating processes come into thermal equilibrium [9,10]. The rate 
of these processes is related to the free energy of the sphaleron-type 
field configurations which carry topological charge [11]. In the SM they 
lead to an effective interaction operator of all left-handed fermions, 
Ob+l = WiQiiQuQii^Lii which violates baryon and lepton number 
by three units. On the other hand, B — L remains conserved. 

2. The weak interactions of the SM violate C maximally and violate 
CP via the Kobayashi-Maskawa mechanism [12]. The latter originates 
from the quark mixing matrix, often called the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix, which contains one CP- violating phase. This 
phase is known to be non-zero since CP violation has been observed 
in the K and B mesons systems (see e.g. the review on CP violation 
in [5]). 

3. A strongly first-order electroweak phase transition in the early Universe 
could provide the out-of-equilibrium condition. A first-order phase tran- 
sition proceeds via nucleation and growth of bubbles [13]. 

This scenario is called electroweak baryogenesis in the Standard Model. How- 
ever, it fails for two reasons. First, it turns out that, for the electroweak phase 
transition to be strongly first order, the mass of the Higgs particle should be 
smaller than about 45 GeV [14] (see also [15]). However, LEP II gives the 
well-known bound > 114 GeV [5]. Second, the source of CP violation in 
the quark sector is far too small, due to the smallness of some of the quark 
masses [16]. 
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In conclusion, successful baryogenesis requires physics beyond the SM, 
just as the dark matter and dark energy problems! One intriguing solution 
to the problem of baryogenesis is deeply connected with the neutrino sector 
and in particular with neutrino masses. This is the topic of the next section. 



1.2 The puzzle of neutrino masses 



1.2.1 Theory of neutrino oscillations 

Even with tiny masses, massive neutrinos can behave very differently from 
massless ones. In particular, massive neutrinos naturally lead to neutrino 
mixing and to neutrino oscillations, which have been recently observed, as 
we shall see below. Let us sketch how this happens. 

Consider a neutrino beam created in a charged current interaction along 
with the antilepton a"'", a = e, ^, r. By definition the neutrino created is 
called z/q,. In general, this is not a physical particle, but rather a superposition 
of physical fields Vi with different masses mi: 



[1.1) 



where U is the lepton mixing matrix, also known as the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) matrix [17-19]. By analogy with the CKM ma- 
trix in the quark sector, the lepton mixing matrix can be conveniently 
parametrized as 

(1.8) 



*2 



U X diag (e' 2 , e' 2 , 1 



V 



( ^ ° 

C23 
V -S23 




Cl3 


-si3e 



i<5 




1 




si3 e 


Cl3 



-i<5 



v 



Cl2 

-Sl2 





Sl2 

Cl2 

1 



where Sij = sin^^, Qj = cos 9ij, S and are the Dirac and Majorana CP- 
violating phases, respectively. The Majorana phases, which are not present 
in the quark sector, are related to the possible Majorana nature of neutrinos 
(see Section 1.3). 

For a simple-minded approach to the propagation of the state {va), we 
assume that the 3-momentum p of the different components in the beam are 
the same. However, since their masses are different, the energies of all these 
components cannot be equal. Rather, for the component z/j, the energy is 
given by the relativistic energy-momentum relation Ei = \/pM-^n|. After a 
time t, the evolution of the initial beam, Eq. (1.7), assuming that neutrinos 
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are stable particles, gives 

K{t)) = Y.e--^'^^'UUv,). (1.9) 

i 

Since all EiS are not equal if the masses are not, Eq. (1.9) represents a 
different superposition of the physical eigenstates Vi compared to Eq. (1.7). 
In general, the state in Eq. (1.9) can therefore show properties of other flavor 
states. The amplitude of finding a in the original beam is 

{l^a'K{t)) = J2^~'''''KiUa'i, (1.10) 
i 

using the fact that = Sij. The probability of finding a in the orig- 

inal Ua beam at any time t is then the modulus squared of the amplitude, 
Pva-^u^iif) = \{t^a'Wa{t))\^- In all practical situations, neutrinos are extremely 
relativistic, so that one can approximate the energy-momentum relation as 
Ei ~ |p| -|- m^/(2|p|) and replace t by the distance x traveled by the beam. 
After a few manipulations one can finally write the vacuum oscillation prob- 
ability as 

Pu^-.u^, (x) = 5aa' - 4 ^ Re {U*iUa'iUajU*,j) Siu^ 

i>j 

+2 J]lm {u:,U^,,U^jU:,j) sin (^x) , (1.11) 
i>j ^ ^ 

where £^ ~ |p| and Amfj = rnf — m|. Prom this formula it is apparent that 
neutrino oscillations require non-zero neutrino masses and mixings. 




1.2.2 Experimental evidence for neutrino oscillations 

The last 10 years have been extremely successful for the field of neutrino 
physics. In 1998, the Super-Kamiokande experiment in Japan reported the 
first compelling evidence for neutrino oscillations as a way to explain the 
anomaly in atmospheric neutrinos. Super-Kamiokande not only confirmed 
the previously found deficit in z/^-type events but also measured a zenith an- 
gle dependent f^-deficit which was inconsistent with expectations based on 
calculations of the atmospheric neutrino flux. The neutrino oscillation expla- 
nation z/^ — with a quasi-maximal mixing angle [20] appeared therefore 
as the most convincing one. 

Two dedicated laboratory experiments have been conceived in order to 
check this picture. The experiment K2K in Japan, which collected data until 
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Figure 1.5: Best fit regions at 90% 
CL. for atmospheric and accelera- 
tor neutrinos [23]. 



Figure 1.6: Best fit regions at 90, 99 
and 99.73% CL. for solar and reac- 
tor neutrinos [23]. 



November 2004, used a pulsed beam of muon-neutrinos produced at KEK 
and detected at Super-Kamiokande (distance of 250 km). The currently run- 
ning MINOS experiment uses a pulsed beam of muon-neutrinos produced at 
NuMI (Fermilab), and the far detector is located at a distance of 735 km 
in the Soudan mine, Minnesota. Both experiments point to a neutrino oscil- 
lation interpretation of their data, with mixing parameters compatible with 
those explaining the atmospheric anomaly [21,22]. A summarizing plot for 
"atmospheric" neutrinos can be found in Fig. 1.5. The best-fit parameters 
are [23]: 

\Aml^\=Aml^ = (2.5 ± 0.2) x 10"^, (1.12) 
sin2 2^23 = 1.02 ±0.04, (1.13) 

where the ranges indicated are at la. 

Neutrinos from the Sun were first studied by Ray Davis in the Homestake 
mine in South Dakota in the 1960's [24]. This pioneering work made use of 
the radiochemical technique, Ve + ^''Cl — > e~ + '^^Ar. He quickly found that 
the observed fiux was smaller than the one predicted by Bahcall and collab- 
orators from the Standard Solar Model [25]. The solar neutrino puzzle was 
born. Later, other experiments were conceived to check this deficit in solar 
neutrinos: Kamiokande [26] and later Super-Kamiokande [27], which used the 
water Cherenkov technique, also found a deficit, as well as SAGE [28] and 
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GALLEX/GNO [29, 30], which used the radiochemical method with GaUium. 

It should be noted that the deficits found in the different experiments 
were different, highlighting the energy dependence, as well as the fact that the 
experiments were based on different techniques. The radiochemical method 
allows only for the measurement of charged-current (CC) events, which is 
exclusively sensitive to electron neutrinos. Experiments relying on water 
Cherenkov detection register only elastic scattering (ES) events, z/+e —>■ v+e. 
This technique can measure in principle all neutrino flavors, but the efficiency 
in the and r channels is much lower, due to the lower cross-sections. 

Only the SNO experiment in Sudbury, Canada, used a heavy water detec- 
tion which allowed for the detection via three channels simultaneously: CC, 
ES as well as neutral current (NC) thanks to the reaction u + d ^ u -{-n + p, 
where d denotes dcuteron. This made possible a reliable measurement of the 
total flux of neutrinos from the Sun. Hence, in addition to conflrming the 
deficit in solar neutrinos in the CC and ES detection channels [31], this ex- 
periment allowed for the first direct test of the Standard Solar Model. In 
2002, the data were found to be actually consistent with the prediction [32], 
finally resolving the longstanding solar neutrino puzzle. This can be seen as a 
huge success in the history of neutrino physics. The Borexino experiment, lo- 
cated in the Gran Sasso laboratory, has recently added another milestone to 
the understanding of neutrinos from the Sun by detecting in real time a fiux 
of ^Be neutrinos consistent with the Standard Solar Model prediction [33] . 

The neutrino oscillation explanation z/g v^^r to the deficits in neutri- 
nos from the Sun in different experiments was popular from the beginning, 
because it appeared as the simplest and most elegant solution. But the first 
experiment which really left neutrino oscillations as the only possible expla- 
nation for the solar neutrino puzzle is KamLAND, rejecting, for instance spin 
flavor oscillations [34,35], which were still allowed after SNO. KamLAND is 
a scintillation detector located in Japan which observed the Ve from neighbor- 
ing nuclear reactors thanks to the process t/g +P ~^ ^ + e"*". This experiment 
saw a 6(7 evidence for the disappearance of Pg [36, 37]. A summarizing plot of 
the solar neutrino and KamLAND data can be found in Fig. 1.6. The best-fit 
parameters are [23]: 

|Arr4| = Am,2„i = (8.0 ± 0.3) x 10-^eV^ (1.14) 
tan^^i2 = 0.45 ±0.05, (1.15) 

where the ranges indicated are at la. 

For completeness, we should add that a short-baseline reactor neutrino 
experiment named CHOOZ, in France, was performed to get a handle on 
the third mixing angle, ^13 [38] . They report no signal, which, together with 
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atmospheric and K2K data, yields the la upper bound [23] 

sin^ 2^13 < 0.05. 



(1.16) 



1.2.3 Absolute neutrino mass scale 

The oscillation formula Eq. (1.11) shows that neutrino oscillation experiments 
are only able to provide information on neutrino mass-squared differences. 
They are insensitive to the absolute neutrino mass scale. Fortunately, there 
are other ways to probe it. 

First, the so-called "direct" measurement makes use of the very precise 
determination of the upper end of the spectral distribution of the electron in 
the tritium /3-decay, '^H — > '^He + i^e + e~ . The Mainz experiment obtained 
the bound m,y^ < 2.2 eV [39] and the Troitsk experiment rrii,^ < 2.5 eV [40], 
both at 95% C.L., where 



with U denoting the lepton mixing matrix we introduced in Section 1.2.1 
[cf. Eq. (1.8)]. The upcoming KATRIN experiment [41], which is under con- 
struction at the moment, is expected to have a discovery potential down to 
rrii,^ ~ 0.35 eV. 

Second, if neutrinos are Majorana particles (see next section), then it 
might be possible to observe neutrinoless double-beta (Oi//?/3) decay for nuclei 
hke ^^Ge. Double-/? decay, (A, Z) (A, Z + 2) + 2e- +29^, has been already 
observed, even though the rate is of second order in the Fermi coupling 
constant G-p, implying for instance a lifetime of about 10^^ years for '''^Ge. 
The Oi//3/3-decay rate is even more suppressed, and depends on yet another 
combination of the neutrino masses as in Eq. (1.17) [42], 



Up to now, no compelling signal^ for Oi//3/3-decay has been observed, and the 
currently running CUORICINO experiment, which uses ^^'^Te, obtains now 
a bound [45] 




(1.17) 




(1.18) 




(90% C.L.) 



(1.19) 



^Part of the Hcidclbcrg-Moscow collaboration claims an evidence for Oi^/?/3-decay, lead- 
ing to the measurement |(m)| = 0.2-0.6 eV (90% C.L.) [43,44], where the uncertainty is 
due to the poorly known nuclear matrix elements. 
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which is shghtly stronger than the one obtained in the previous Heidelberg- 
Moscow [46] and IGEX experiments [47, 48], where ^^Ge was used. Note that 
the range shown in Eq. (1.19) stems from the uncertainty in the nuclear 
matrix elements. Future planned experiments such as GERDA [49], MAJO- 
RANA [50] or CUORE [51] aim at a sensitivity in the 50 meV range. 

There is a third probe of the absolute neutrino mass scale, and this is 
cosmology. More precisely, the CMB data and the data from the large scale 
structure (LSS) have a sensitivity to the sum of the neutrino masses ^ifrii. 
Looking at the literature in the last 3 years that dealt with this issue, it is 
clear that the bound on the sum of the neutrino masses depends on which 
data sets are included. We simply show here the result from a somewhat 
conservative calculation in [52]: 



The three probes we have presented involve different combinations of the 
neutrino masses, where even cancellations are possible in the case of Ou(3f3- 
decay [cf. Eq. (1.18)]. But the bottom line here is that all three probes tend 
to the same conclusion: the neutrino mass scale has to be in the sub-e V range. 
In this sense, this makes neutrinos a very peculiar particle in the SM, with 
a mass that is six orders of magnitude smaller than the electron! Fig. 1.7 
illustrates nicely this gap between neutrinos and the other SM particles. 
Clearly, such a gap must be explained on theoretical grounds, and this is the 
topic of the next section. 

1.3 An elegant solution: 

the see-saw mechanism and leptogenesis 

As we have just seen, neutrinos are very special particles due to their extreme 
lightness. This might suggest that neutrinos have a unique mass generation 
mechanism that leads to small masses in a natural way. One of the first 
attempts in the late 1970's and early 1980's was the use of the see-saw mech- 
anism [53-58], which we shortly describe now. 

Let us imagine that we would like to write down the most simple and 
general extension of the Standard Model that accomodates neutrino masses. 
For simplicity we consider here only one generation of neutrinos. In the SM, 
in order to generate massive fermions, one introduces a right-handed (RH) 
component for each particle. So let us introduce one RH neutrino in a 
usual Dirac-type term m-oVLNji -\- h.c. Contrary to the other SM particles. 



3 




(95% C.L.). 



(1.20) 



1=1 
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Figure 1.7: Mass spectrum of fermions. 



this is not the only term allowed for neutrinos. Since neutrinos are neutral 
particles, the RH component can only be a singlet of all SM gauge interac- 
tions, including hypercharge. This unique property of neutrinos implies that 
the Majorana mass term ^Mm{NjiYNr + h.c. is allowed, where the super- 
script c denotes charge conjugation. Actually, in general there is no reason 
why this term should be absent. The mass term for neutrinos can then be 
conveniently written in a matrix form as 



1 
2 



{Nnr ) 



mD 
TJiB Mm 



Nr 



+ h.c. 



[1.21) 



Note that upper left component of the mass matrix is zero because it is not 
possible to write down a Majorana mass term for left-handed fields without 
breaking gauge invariance. 

The see-saw mechanism then assumes that the Majorana mass is much 
larger than the Dirac mass, i.e. Mm ^ mo- Diagonalizing the mass matrix 
with this assumption yields the two mass eigenvalues 



m 



m 



D 



M^ 



M, 



M 



(1.22) 



where the first eigenvalue corresponds to the mass of a light Majorana neu- 
trino, and the second eigenvalue corresponds to the mass of a heavy Majorana 
neutrino. It is interesting to notice that plugging a Dirac mass at the elec- 
troweak scale, mo = 100 GeV, with a heavy neutrino mass of Mm = 10^^ GeV 
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maybe related to the underlying Grand Unified Theory (GUT), one obtains 
a light neutrino mass m = 0.1 eV. This can be seen as a strong support for 
this model. For the generalization of the previous discussion to three massive 
neutrinos and more details, see Appendix B. 

The version of the see-saw mechanism we have just described, with the 
addition of RH neutrinos to the SM, is nowadays referred to as the type-I 
see-saw mechanism. For completeness, let us mention that other see-saw 
mechanisms were proposed to explain the smallness of neutrino masses: the 
addition of SU{2) triplet Higgses to the SM leads to the type-11 mecha- 
nism [58-60], and the addition of SU (2) triplet fermions leads to the typc-III 
mechanism [61, 62]. Throughout the thesis, we shall exclusively deal with the 
type-I mechanism. 

As we have just shown, the see-saw mechanism elegantly solves the prob- 
lem of generating small neutrino masses in a natural way. But this is not the 
only virtue of the model: its cosmological consequence is indeed that it can 
also solve the puzzle of the matter- antimatter asymmetry of the Universe 
through thermal leptogenesis, as first proposed by Fukugita and Yanagida 
20 years ago [63]. We saw in Section 1.1 that one needs to satisfy all three 
Sakharov's conditions in order to solve the baryogenesis problem. Let us see 
how they are satisfied in the case of thermal leptogenesis. 

1. Lepton number is violated in the decays of the heavy neutrinos into 
lepton doublets and Higgs doublets. Additionally, B and L are violated 
by the non-perturbative sphaleron processes, which however conserve 
B—L. This is thus the same source of B violation as the one described in 
the case of electroweak baryogenesis (see Section 1.1). The difference is 
that the sphaleron transitions are here transferring L number produced 
in the decays of the heavy neutrinos into the required B number. These 
transitions should be in equilibrium for temperatures [64-66] 

Tew ^ 100 GeV <T< 10^^ GeV. (1.23) 

2. CP violation occurs in the decays of the heavy neutrinos. The relevant 
quantity to evaluate is the CP asymmetry parameter, defined as 

r - r- 

£. = -i^^, (1.24) 

where = r{Ni and F^ = J^a ^i^i ^ ^a'^')- As shown in 

Fig. 1.8, one can calculate such an asymmetry by computing the inter- 
ference between the tree-level diagram and the two relevant one-loop 
diagrams, namely the self-energy diagram and the vertex correction. 
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Figure 1.8: The diagrams necessary to compute the CP asymmetry in lepto- 
genesis. 



We will show the result of such a calculation in the next chapter. This 
source of CP violation is typically sufficient to explain the matter- 
antimatter problem. 

3. The out-of-equilibrium condition will be satisfied thanks to the expan- 
sion of the Universe. Convenient quantities to describe when the decays 
of the heavy neutrinos freeze out are the decay parameters, Ki, defined 
as the ratio of the total decay width of the heavy neutrino to the 
expansion rate at T = Mj, 

rD,i(T = 0) 

where T-Q i = Fj + Fj. This is the key quantity for the thermodynamical 
description of the decays of heavy particles in the early Universe [2] . In 
leptogenesis it can be conveniently expressed in terms of the effective 
neutrino mass [67] 

rrii = — ^Kivn^, (1.26) 

where is the equilibrium neutrino mass, given by 

~ 1.08 X 10~^ eV. (1.27) 

In this thesis, when referring to "leptogenesis" we will exclusively mean 
the production of a lepton asymmetry by the decays of heavy singlet neutri- 
nos. Other "leptogenesis" scenarios, such as the Affleck-Dine mechanism in 
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supersymmetric theories [68-70] or leptogenesis via neutrino oscillations [71, 
72] will not be discussed here. 

After the clear evidence for non-zero neutrino masses appeared in 1998, 
the see-saw mechanism, and hence leptogenesis, became very attractive, solv- 
ing simultaneously two important puzzles of modern physics within a simple 
and minimal extension of the SM. Another nice feature of the model is that 
the heavy fermionic singlets (RH neutrinos) necessary for the (typc-I) sec-saw 
mechanism have a natural connection with grand unification; e.g. 5*0(10) [73] 
predicts the existence of such fermionic singlets. It is therefore not a surprise 
if a huge activity was registered in this field since 1998, with an ever growing 
precision in the computations. 

We present in Chapter 2 the picture of leptogenesis that emerges when the 
fiavor content of the leptons coming from the decays of the heavy neutrinos 
is ignored. This is what we call the "unfiavored" treatment. In Chapter 3 
we discuss why fiavor effects are important in a considerable part of the 
parameter space, and what their main implications are. This will allow us to 
build up a "fiavored" picture of leptogenesis. In Chapter 4 the validity of the 
classical computation of the final baryon asymmetry in the case of fiavored 
leptogenesis is analyzed. We shall argue that the use of a more complicated 
density matrix equation might be relevant in some region of the parameter 
space. Then, in Chapter 5, an overview of different ways to go beyond the 
conventional picture of leptogenesis, where only the lightest RH neutrino is 
taken into account and the mass spectrum of the heavy neutrinos is very 
hierarchical, is given. In Chapter 6 we consider a by-product of fiavor effects, 
namely the fact that the low-energy phases in the PMNS matrix have a more 
important role than previously thought. In particular, the role played by the 
Dirac phase, source of CP violation in neutrino oscillation experiments, will 
be thoroughly studied. Finally, we conclude with a summary of the main 
results and a discussion about what one can expect from the experimental 
side in the next few years regarding a possible test of leptogenesis. 



Chapter 2 

Vanilla leptogenesis 



In this chapter, following a purely unflavored^ treatment, we first introduce 
the most general Boltzmann equations which need to be solved in order to 
estimate the baryon asymmetry of the Universe produced through leptoge- 
nesis. This will allow us to introduce very useful quantities and definitions 
which will be used in the remainder of the thesis. Then, we describe how the 
general picture simplifies considerably under natural assumptions, such as a 
hierarchical mass spectrum for the heavy neutrinos. This will enable us to 
arrive at a "vanilla" leptogenesis, i.e. a typical scenario of leptogenesis, where 
the lightest RH neutrino, A^i, provides the main contribution. Wc will then 
derive two important constraints on the parameters of the model: a lower 
bound on the mass of the lightest RH neutrino, Mi, and an upper bound 
on the absolute neutrino mass scale, mi. Finally, we shortly comment on 
leptogenesis within a supersymmetric framework and discuss the so-called 
gravitino problem. 



2.1 General scenarios of unflavored leptoge- 
nesis 

Wc would like here to describe how to calculate the baryon asymmetry pro- 
duced through leptogenesis in full generality, i.e. for an arbitrary RH neutrino 
mass spectrum and no restrictions on the parameters of the model. 

The Lagrangian in Eq. (B.l) contains all the terms relevant for leptoge- 
nesis. Since we are concerned here with a purely unflavored analysis, it will 

^In the literature the expressions "flavor-independent" or "one- flavor" are often used 
to describe the same concept. 
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prove useful to define the states 



and 



Ki)^^^i^E^-l^«) (2-1) 



K)^^7TfTrE^-l^"«)' (2-2) 



which can be thought of as the states produced in decays and inverse decays 
(and AL = 1 scattering processes, see below) in the process of leptogenesis. 
So, to each heavy neutrino corresponds one lepton state and one anti- 
lepton state Note that is not the CP conjugate of although \Ia) 
is the CP conjugate of \la). We shall discuss this point in more detail in the 
next chapter and see which important consequences follow. 

Let us hst now the processes that are relevant for the computation of the 
baryon asymmetry through leptogenesis [67, 74-76] : 

• Decays and inverse-decays: Ni and the conjugate processes 

• AL = 1 Higgs-mediated scattering processes, such as the s-channel 
£iNi Q^t and the t-channels liQ^ ^ Nit, and iit NiQ^, where 
Qs and t are the third generation quark doublet and the top SU{2) 
singlet, respectively, as well as those involving gauge bosons, such as 
liNi $U (with A = W^'^ or B). 

• The s-channel AL = 2 scattering processes i^^^ with on-shell Ni 
are already accounted for by decays and inverse decays. However, the 
off-shell s-channel contribution, as well as the li-channel and t-channel 
scatterings {££ ^ ^^^^), must be included. 

For simplicity, we shall neglect in the following the so-called spectator 
processes [77,78], which arc processes that arc not directly related to the 
generation of the asymmetry but that are fast and could affect it indirectly. 
For example, the sphaleron transitions belong to spectator processes, since 
they are supposed to be in equilibrium during the leptogenesis process and 
since they are responsible for transmitting the L asymmetry produced in 
the decays into a B — L asymmetry in a non-trivial way. For the time being, 
we neglect spectator processes, and simply assume that the asymmetry is 
directly produced in B — L instead of L. In any case, spectator processes 
altogether are not expected to induce corrections larger than 30% [78]. We 
shall discuss in detail spectator processes in the next chapter (Section 3.2). 
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Thermal corrections will also be neglected. They may lead to relevant 
corrections, though with big theoretical uncertainties, in the weak washout 
regime {Ki < 5), but negligible corrections are expected in the more relevant 
strong washout regime {Ki > 5) [75]. 

Following an unflavored treatment, the set of Boltzmann equations can 
be worked out in the form [67, 74, 79] 

= -{Di + Si)iN^,-N'^^), z = 1,2,3 (2.3) 



dz 

dNB-L 

dz 



J2 ^i{Di + Si) {N^^ - N'^) - Nb-l W , (2.4) 



where z = Mi/T. With Nx we denote the abundance of X per RH neutrino in 
ultra- relativistic thermal equilibrium. Defining Xi = Ml/Ml and Zi = z -s/xl, 
the decay factors are given by 

Di = ^^K,x,z l-\ , (2.5) 
Hz \7i/ 

where H is the expansion rate and the decay parameters Ki were introduced 
in Eq. (1.25). The total decay rates, F^.i = F,; + Fj, arc the product of the 
decay widths times the thermally averaged dilation factors (l/7i), given by 
the ratio lCi{zi) / 1C2{zi) of the modified Bessel functions. 

As we have seen, the AL = 1 scatterings Si and the related washout 
contribution W,f^^^^ arise from two different classes of Higgs- and lepton- 
mediated inelastic scatterings involving the top quark (Q3) and gauge bosons 
(^4) . Their main effect is to enhance the heavy neutrino production and thus 
the efficiency factor in the weak washout regime, and, since they also con- 
tribute to the washout, they lead to a correction of the efficiency factor in 
the strong washout regime as well. Top quark and gauge boson scattering 
terms are expected to be of similar size. However, the reactions densities for 
the gauge boson processes are presently controversial [75,76]. Therefore, we 
will not discuss these processes here. A simple analytic approximation for 
the sum Di -\- Si was obtained in [80] , 



Di + Sic^ 0.1 Ki 

where 



10 



(2.6) 



(2.7) 



ln(Mi/M,)- 

The Higgs mass Mh was introduced to cut off the infrared divergence of the 
t-channel process. It turns out that Di + Si ~ Di in the strong washout 
regime. The generalization of Eq. (2.6) for i 7^ 1 can be derived from [76]. 
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The equihbrium abundances of the heavy neutrinos and their rates are 
also expressed through the modified Bessel functions, 

1 dA^^*^ 1 

The RH neutrinos can be produced by inverse decays and AL — 1 scatterings. 
Nevertheless, in the relevant strong washout regime, inverse decays alone are 
sufficient to make the RH neutrino abundance reach its thermal equilibrium 
value prior to the decays which produce the asymmetry. The asymmetry 
generated together with the RH neutrino production will then be efficiently 
washed out, so that the details of the RH neutrino production will not affect 
the final asymmetry, thus greatly reducing theoretical uncertainties. This is 
one of the nice features of the strong washout regime on which we shall focus. 
For this reason, the AL = 1 scattering terms Si will play a subdominant role. 

The washout factor W in Eq. (2.4) can be written as the sum of two 
contributions [81], 

W = J2^iiKi) + ^W. (2.9) 

i 

The first term is the sum of the contributions from inverse decays and AL — 1 
scatterings, 

Wi = Wl^ + W,^^=\ (2.10) 

For the lightest RH neutrino, Ni, it was shown in [80] that Wi can be con- 
veniently rewritten as Wi{z) = j{z)Wl^{z), where 



.(^) = 0.l(l + f 



Mh V z) z 



(2.11) 



where /i = 1 (2/3) in the strong (weak) washout regime. The corresponding 
results for i 7^ 1 can again be obtained from [76]. 

In the strong washout regime, inverse decays, where the resonant AL = 2 
contribution has to be properly subtracted [75,82], dominate AL = 1 scat- 
terings [75,80], so that 

Wiiz) ~ Wl^'iz) = ^Ki ^iKM) zf . (2.12) 

The second term in Eq. (2.9) arises from the non-resonant AL = 2 
processes and gives typically a non-negligible contribution only in the non- 
relativistic limit, for 2; ^ 1 [74,80,81]. For hierarchical light neutrinos, it 
can be safely neglected for reasonable values Mi <^ 10^^ GeV. We shall come 
back to this contribution when discussing the upper bound on the absolute 
neutrino mass scale at the end of Section 2.2. 



2.1 General scenarios of unflavored leptogenesis 



21 



The effects of production and washout are simultaneously accounted for 
by the efficiency factors associated with each Nj,. Let us indicate with N^^-l 
a possible pre-existing asymmetry at the initial temperature of leptogenesis 
Tin- The final B — L asymmetry can then be written as [2,80], 

N'b-l = N'i-L exp (|- E / ^^(^')) +Y.^^<^ (2-13) 
with the final efficiency factors k\ = Ki{z ^ oo) given by 

kI = - dz' ^ exp J2 £ dz" W,{z")^ , (2.14) 

where we defined Zi^-^ = Mi /Tin. Notice that each efficiency factor depends in 
general on all decay parameters, i.e. k,1 = k,1{Ki, K2, K3). 

The baryon-to-photon ratio at recombination can then be calculated as 

VB = asph^^ ^ 0.96 X lO-^TV^.^ , (2.15) 

where N^^^ ~ 37, and Ogph = nB/riB-L accounts for the sphaleron con- 
version of a B — L asymmetry into a B asymmetry. If the electroweak 
sphalerons go out of equilibrium before the electroweak phase transition, 
one has Osph = 28/79 [83]. If, instead, electroweak sphalerons remain in equi- 
librium until slightly after the electroweak phase transition (as would be 
the case if, as presently believed, the electroweak phase transition was not 
strongly first order), this factor would be Ogph = 12/37 [84]. Both coefficients 
are of order 1/3 and, for definiteness, we shall use Ugph — 28/79 throughout 
the thesis. 

If the mass differences satisfy the condition for the applicability of per- 
turbation theory \Mj-Mi\/Mi > max[(/i'l' h)ij]/ {16 n"^) with j ^ i [85], then 
a perturbative calculation from the interference of tree-level with one-loop 
self-energy and vertex diagrams (see Fig. 1.8) gives [86, 87] 

_ 3 ^ Im [(/.t h)l] ^{x,/x,) 



where the function ^(x), shown in Fig. 2.1, is defined as [88] 

1 + x\ 2 — X 



2 



(1 + x) In 



X / L — X 



(2.17) 
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Figure 2.1: The function ^{x) defined in Eq. (2.17). 



A particularly useful paramctrization of the Yukawa coupling matrix in the 
context of leptogenesis involves the orthogonal matrix fl [89]^, 



where we defined Dm = diag(Mi, M2, M3) and Dm = diag(mi, m2, itls). The 
matrix U diagonalizes the light neutrino mass matrix rriiy given in Eq. (B.5), 
i.e. W nil, U* = Dm, and it can be identified with the lepton mixing matrix in 
a basis where the charged lepton mass matrix is diagonal (see Appendix A) . 
Moreover, neglecting the effect of the running of neutrino parameters from 
high energy to low energy [90, 91], one can assume the U matrix to be iden- 
tified with the PMNS matrix, partially measured in neutrino oscillation ex- 
periments. In the following, we shall refer to the parametrization Eq. (2.18) 
as the "orthogonal" or the "Casas-Ibarra" parametrization. 
The fl matrix can be conveniently parametrized as 



(2.18) 





where 




(2.20) 



^Compared to the R matrix in [89], one has the simple relation 0. = E) . 
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i?i3 = 1 , (2.21) 

\ c^si a/I - '^ii / 
/I \ 

i?23 = x/l - CU|2 -CU32 . (2.22) 

\ UJ2.2 \/l - a;|2 / 

This parametrization for an orthogonal complex matrix corresponds to the 
transposed form of the CKM matrix in the quark sector or of the PMNS 
matrix in neutrino mixing [cf. Eq. (1.8)], with the difference that here one 
has complex rotations instead of real ones. 

Notice that, using the orthogonal parametrization, Eq. (2.18), the de- 
cay parameters can be expressed as linear combinations of the neutrino 
masses [80, 92], 

Ki-Y.—\^l\- (2-23) 

3 

The parametrization Eq. (2.19) is especially useful to understand the 
general structure of different scenarios of leptogenesis. 

• For Q = it!i3, one has £2 = 0, while Si is maximal if [93] 

m^Re{ujl^)/\ujl^\=mr[l-Re{ujl^)]/\l-ujl^\. (2.24) 

In the hierarchical limit for the heay^ neutrinos (HL), i.e. 
Ml ^ M2 <^ M3, one obtains the A^i-dominated scenario (A^iDS), 
where the final asymmetry is the result of only A^i-related processes. 

• For Q = i?23, one has £1 = 0, while £2 is maximal if 

m2 Re{ujl^)/\ujl^\ = m^[l - R^{ujl2)]/\l - ^32!- (2-25) 

At the same time, one has mi — mi, so that the washout from A^i can be 
neglected if mi is small enough. Therefore, in the HL and for hierarchi- 
cal light neutrinos, one obtains the A^2-dominated scenario ( A^2DS) [93] , 
which will be the topic of Section 5.2. 

• If = i?i2, then El undergoes a phase suppression compared to its 
maximal value but |£2| (M1/M2) This implies that in the HL 
one again recovers the iViDS [93]. On the other hand, if Mi ~ M2, both 
A^i and N2 should be taken into account since both CP asymmetries 
are expected to be of the same order. 
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Before ending this section, we would like to make two remarks. First, using 
the Casas-lbarra parametrization, Eq. (2.18), it can be easily seen that the 
PMNS matrix U cancels out in combinations like {h'^h)ij (or [m\^m-Q) ij) . This 
implies that both the CP asymmetry, Eq. (2.16), and the decay parameters 
or effective neutrino masses, Eq. (1.26), do not depend on U . Hence, the final 
asymmetry is completely independent of the PMNS matrix. In particular, if 
the CP asymmetry is insensitive to the CP-violating phases in [/, then the 
source of CP violation necessary for leptogenesis must come from the "high- 
energy" sector, i.e. the Vl matrix, which is not probed in neutrino experiments. 
We shall see in the next chapter that the situation is completely different 
when flavor effects are taken into account. 

Second, the border between the A^iDS and the iV2DS will also be affected 
by flavor effects. Furthermore, contrary to the unflavored case, where £3 is 
always suppressed by factors Mi/M2,3, the flavored CP asymmetries £3^ are 
not suppressed, opening the way to a potential A^3-dominated scenario. We 
shall nonetheless see, for instance in Section 6.2.1, that, even though the 
CP asymmetry may not be suppressed, the washout from the other two 
heavy neutrinos is difficult to avoid. 

2.2 The ATi-dominated scenario 

Let us now discuss the iViDS, where the asymmetry is dominantly generated 
by the hghtest RH neutrino, A^i. This typically (but not necessarily) occurs 
when a hierarchical heavy neutrino spectrum is considered, 3 Mi < M2 [94]. 
The general expression for the final asymmetry, Eq. (2.15), reduces to 

r/B^10-'£i4, (2.26) 

where k\ can be calculated solving a system of just two kinetic equations [2, 
75,81], 

-(£>! + 5i)(7V^,-7V^J, (2.27) 

Si {Di + Si) {Nn, - TV^J - Wi Nb-l . (2.28) 

These equations are obtained from the general set, Eqs. (2.3) and (2.4), 
neglecting the asymmetry generation and the washout terms from the two 
heavier RH neutrinos. 

For Ml ^ 10"*^^ GeV (^riatm/ X] )> the term ^W{z) in the washout term 
[cf. Eq. (2.9)] is negligible and the solutions depend just on A'l, since this is 
the only parameter in the equations. The B — L asymmetry can be worked 



diV 
dz 

dNB-L 

dz 
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out in an integral form [2] , and one obtains a special case of the more general 
Eq. (2.13), 



Nb-l{z; z) = Nb-l exp 



dz'Wi{z') ) +eiKi{z]z) , 



(2.29) 



where now a possible asymmetry produced by the two heavier RH neutrinos 
and frozen at z > is included in Nb-l- The efficiency factor ki{z; z) can 
be expressed through a Laplace integral, 



Ki{z-Z) = - j 



az — — — exp 
az' 

dz' exp [—■^(2;', z)] 



dz"W^{z") 



(2.30) 



In the strong washout regime we are interested in, for > 5, one can 
safely use the approximations dN^^/dz' ~ dA^^/dz', Di -\- ~ Di, 
Wi{z') ~ Wl^{z') [cf. Eq. (2.12)], and one finds that the final value for the 
efficiency factor, when 2; — > 00, is given by [80] 



KiZb{K,) 



1 - exp ( --KiZb{Ki) 



(2.31) 



if ^ ^ -2b — 2. The value z' — zb{Ki) is where the quantity '^(2;', 00) has 
a minimum and the integral in Eq. (2.30) receives a dominant contribution 
from a restricted 2;'-interval centered around it. In the strong washout regime, 
it can be calculated as a solution of 



W^izB) 



d'N'^Jdz' 



Z=ZB 



(2.32) 



For very large Ki, the right-hand side of this equation tends to unity and 
Zb — Zos, the value of z when the washout from inverse decays switches 
off, i.e. Wib{z > Zos) < 1- Fig. 2.2 shows (dashed lines) that Eq. (2.31), 
with zb{Ki) given by Eq. (2.32), reproduces the numerical result (solid line) 
within 10% for Ki>3 [80]. 

Even though the approximation dNj^^/dz' ~ dN^^/dz' works rigorously 
only in the strong washout regime, Eq. (2.31) describes also the correct weak 
washout regime for a thermal initial A^i-abundancc (iV™^ = 1), because k,{ 
depends only on the value of the initial abundance and not on the decay rate. 
However, in the intermediate regime {Ki ~ 1) the error is about 30%. For 
Ki < 1, the approximation dNi^-^^/dz' ~ dN'^^/dz' docs not work well and 
zb{Ki), evaluated with Eq. (2.32), incorrectly saturates to a constant value 
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Figure 2.2: Efficiency factor in the A'^iDS. The solid lines are the numerical 
solutions of Eqs. (2.27) and (2.28) with AW = 0. The analytical expres- 
sion (2.31) yields the dashed line if zb{Ki) is given by Eq. (2.32) and the 
circled line if it is given by Eq. (2.34). 



1.33 in the limit Ki 0, corresponding to the maximum of dN'^_^/dz 
(see upper panel of Fig. 2.4). As a matter of fact, in the weak washout regime 
the maximum of the asymmetry production does not occur at z^^^ but at 
higher values 



weak 
''max 



(K,) ^ + 15/8, 



(2.33) 



roughly when the age of the Universe is equal to the RH neutrino lifetime. 
Indeed, in the weak washout regime, the iVi's decay far from equilibrium, 
when inverse decays can be neglected. In this case one has approximately 



A /"weak / 



2), where [80] 



/"weak / 



iV- exp 



iVi^, exp <^ -K, 



152; /15 ^ _^ 

2 8 ■ V 8 / ^ V 15^ 



with N^j!^^ = 1 for a thermal initial A'^i-abundance. The maximum of 
^jyweak/^^ giveS the value of z where the production of the asymmetry is 
maximum in the weak washout regime, and one can easily find that this value 
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is given approximately by z^^{Ki) [cf. Eq. (2.33)]. An example is shown in 
the upper panel of Fig. 2.3 for Ki = 0.01. An improvement of the analytical 
expression in Eq. (2.31) is obtained replacing Zb{Ki) from Eq. (2.32) with an 
expression that coincides with it at large Ki and with ^^max — 2 at i^i ~ 1, 
such as 

ZBiK,) ~ 2 + 4Xi°-i^ exp (-'^] . (2.34) 



Fig. 2.2 shows (circles) that this expression, plugged into Eq. (2.31), repro- 
duces the numerical solution (solid line) with a precision always better than 
10%. 

Notice that, in the particularly relevant range 5 < i^i < 100, Eq. (2.31) 
is well approximated by a power law [95] , 

k{{K,)c^^. (2.35) 

In the case of a vanishing initial iVi-abundance {N^^i — 0)' to 
take into account two different contributions, a negative and a positive one, 

k{{K^) = kI{Ki) + 4(^1) . (2.36) 

Defining Zeq by the condition A^iVi(-Zeq) = -^^(-^eq), the negative contribu- 
tion arises from a first stage when N^^ < N^^, for z < Zeq, and is given 
approximately by [80] 



kI{Ki) ~ -2 exp 



exp ( ^iV^^(^eq) 



(2.37) 



8 

The A^i-abundance at Zeq is well approximated by 

Nn, (^eq) ^ mi) ^ , ^^^'^ ,2 , (2-38) 

interpolating between the limit Ki 3> 1, where Zeq ^ 1 and N]s[j^{zeq) = 1, 

and the limit Ki <^ 1, where Zeq ^ 1 and A^iVi(-Zeq) = N{Ki) = SnKi/A. 

The positive contribution arises from a second stage when N^^ > N^^, for 
z > Zeq, and is approximately given by [80] 



4(^i) - 



1 - exp ( -^i^i^B(i^l)A^JVi(^eq; 



(2.39) 



As can be seen in Fig. 2.2, the use of Eq. (2.34) yields an improvement also 
for a vanishing initial A^i-abundance. 
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Figures 2.3 and 2.4 show, for a thermal initial A^i-abundance, the dynam- 
ics of the asymmetry generation, comparing one example of weak washout 
with one example of strong washout [94] . In the top panels we show the func- 
tion dni/dz' = exp[—ip{z',z)], defined for z' < z, for different values of z. 
The difference between the two cases is striking. In the weak washout regime, 
each decay contributes to the final asymmetry for any value of z' at which 
the asymmetry is produced. In the strong washout regime, all the asymmetry 
produced at z' < zb — 2 is efficiently washed out by inverse decays, so that 
only decays occurring at z' ~ zb give a contribution to the final asymmetry. 

After having thoroughly studied the efficiency factor, we would like to 
turn now to the other important piece in the calculation of the baryon asym- 
metry Eq. (2.26), namely the CP asymmetry ei. The general expression (2.16) 
for El can be re-cast through the fl matrix as [93] 

Si = ^{X2) ef^imi, Ml, Q21, ^31) + [{(a^s) - {(2^2)] A£i(mi, Mi, Q21, ^^31, ^^22) ■ 

(2.40) 

In the hierarchical limit (HL), for .X3,X2 ^ 1, one has ^(.X2) — ^{xs) ~ 1, 
yielding Si ~ £:^^(mi. Mi, Q21, ^31)- Therefore, the dependence on four of the 
see-saw parameters, namely M2, M3 and ^22, disappears in the HL, and one 
is left with only six parameters. Notice moreover that k,{ = K,{{mi, Mi, Ki), 
where Ki — Kilrrii, D,2i, ^31), and thus the final asymmetry depends on the 
same six parameters. In particular, let us emphasize once more that the final 
asymmetry is independent of the PMNS matrix, which contains three mixing 
angles and three phases [cf. Eq. (A.5)]. 
Let us define 

and 

V, lm[(]2 ] 

^atm Z^j 

The HL for Si can then be written as [96] 

£fL(mi, Ml, Q21, ^^3i) = £(Mi) /3(mi, Q21, ^^31)- (2.43) 

It is interesting that /9 (mi, 1^21; ^31) ^ 1) so that in the HL one has the 
upper bound \s^^\ < s{Mi) [97]. A more precise bound, which is known as 
the Davidson-Ibarra bound, was later derived [96], 



(2.44) 
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Figure 2.3: Dynamics in the weak washout regime for a thermal initial A^^i- 
abundance (iVj^^ = 1). Top panel: rates. Bottom panel: efficiency factor ki 
and A^i-abundance. The maximum of the asymmetry production rate occurs 
at z' ~ z^^^iKi = 0.01) ~ 12 [cf. Eq. (2.33)]. 
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Figure 2.4: Dynamics in the strong washout regime. Top paneh rates. Bottom 
paneh efficiency factor ki and A^i-abundance. The maximum of the final 
asymmetry production rate occurs at ~ zb- 
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where we used the fact that {rris — mi)/matm = ?7?.atm/(?7ii + ms). Further 
refinements were added, leading to the even more accurate mi-dependent 
bound, 

kf^l <£(Mi)^^^^/(mi,mi), (2.45) 
mi + nis 

where [88, 98] 

~\ f l-mi/nii ifmi<m3, 
/(mi,mi)^< /- 7 — ,~ .„ (2-46) 

obtained by maximizing over the f2-parameters at fixed mi. From this bound 
one gets that /(mi, mi) and hence the CP asymmetry are maximal in the 
limit mi — > 0. One notices as well that the CP asymmetry vanishes when 
mi = mi. 

The presence of the bound on the CP asymmetry Eq. (2.45) motivates 
the introduction of the so-called "effective leptogenesis phase" 5l as 

/3(mi, ^21, ^^3i) = /3max(mi, mi) sin 5i,{mi, Q21, ^31) , (2.47) 

where 

/3max(mi,mi) = ^-^/(mi,mi), (2.48) 

mi + m3 

such that the upper bound in Eq. (2.45) corresponds to sin^L = 1. 
It is also useful to express the function /(mi, mi) as [93] 

/(mi, mi) = Y^^{mi, mi) , (2.49) 

mi 

where ymax(^i, ^1) represents the configuration of fl parameters which max- 
imizes the CP asymmetry at fixed mi. It turns out that such a configuration 
always occurs for Q21 = [93], and hence l^ax is the maximum of Im[Q|^] at 
fixed mi. In other words, for ^21 = and ymax("^i, mi) = Im[fi3]^], the phase 
5l is maximal, while for a generic choice of fl, the CP asymmetry undergoes 
a phase suppression 

sin5L(mi, Q21, ^^3i) = ^^1+^3 (ini[Q2 ] _^ (7^lm[Q^J) 

mi /(mi, mi) 

Ynmx{mi,mi) 

where a = \/v[i^ — m\/m.^x.^. It can be readily seen that sint^L = 1 for 
Im[r22i] = and ImfQg^^] = Fmax- 

It is instructive to calculate sin 5l for each of the three elementary complex 
rotations that can be used to parametrize Q, [cf. Eq. (2.19)]: 
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• For fl = i?i3, one has sin — Imffigj^J/Ymax and the phase is maximal if 
Im[n|^] — Fmax = 'nT'i/rria^tm', notice that there is no difference between 
normal and inverted hierarchy. 

• For Q = i?i2, one has sin^L = c^Iin[Jl|i]/yinax < cr, larger for in- 
verted hierarchy than for normal; however, for fully hierarchical light 
neutrinos (mi <^ mso\), one has that Ki = Ksoi |^2il ^'^^ normal hierar- 
chy and Ki ~ Katm. I^iil ^r inverted hierarchy, and since k,{ oc K^^''^ 
[cf. Eq. (2.35)], the final asymmetry is slightly higher for normal hier- 
archy compared to inverted at fixed |Q|i| [93]. 

• For Vt = i?23, one has sin^L = = 0; one can check that ei = 
applies independently of M2 and M3 and therefore not only in the HL. 
Notice that the conclusions in the previous two cases are still valid if 
one multiplies R13 or R12 with R23 respectively, since it does not affect 
sin Si,. 

Interesting constraints follow if one imposes that the asymmetry produced 
from leptogenesis explains the value of the baryon-to-photon ratio inferred 
from CMB observations, Eq. (1.6), [4,81] 

r/B(mi, Ml, Q21, ^^3i) = vT^ = (6.1 ± 0.2) x 10"^° . (2.51) 
If Ml < 10^^ GeV {mlJ ^. m^), then 

Ml 

Ml = = 

K'{Ki) Prn3x{mi,Ki) sin5L(f^2i,f^3i) 

> Mr{Ki) = —--^ — , (2.52) 

where we introduced the quantity 
Mi = ^ = (6.25 ± 0.4) X 10« GeV > 5 X lO^GeV . (2.53) 

3 Ogph ?7T.a,tm 

The last inequality gives the 3a value that we used to obtain all the results 
shown in the figures.^ Eq. (2.52) is quite general and shows the effect of 
the phase suppression [93] and of a higher absolute neutrino mass scale [95] 
in making Mi higher. In Fig. 2.5 we show M™" (thick solid line) for fully 
hierarchical light neutrinos (mi = 0) and maximal phase (sin^L = 1)- It 



^Notice that Mi gives the lower bound on Mi for a thermal initial iVi-abundance in 
the limit Ki ^ 0. 
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is convenient to introduce a value such that, for Ki > K^^, the final 
asymmetry calculated for a thermal initial A'^i-abundance {N^^_^ = 1) differs 
from the one calculated for a vanishing initial iVi-abundance {N^^^ = 0) by 
less than some quantity S. When S — 10%, ~ 3.5 and one obtains the 
lowest value [93], 

Ml > 4 X lO'^ GeV . (2.54) 

The lower bound on Mi also translates into a lower bound on Tin, the initial 
temperature of leptogenesis, 

( M™^'^1 

T.n > (T/Dhl ^ > 1-5 X 10^ GeV {K, > 3.5) . (2.55) 

A plot of this lower bound is shown in Fig. 2.5 (thick dashed line). The 
relation between M™"^ and 7]™"^ can be understood from the top panel of 
Fig. 2.4, which shows that the final asymmetry is the result of the decays 
occurring just around zb, when inverse decays switch off, whereas all the 
asymmetry produced before is efficiently washed out. 

Assuming a period of inflation at early stages, the minimal initial tem- 
perature Tj™™ that allows for successful leptogenesis can be identified with 
the minimal reheat temperature T™™ after infiation. In locally supersym- 
metric theories, the high temperature required poses a problem known as 
the gravitino problem, as we shall see in the next section. 
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For increasing values of the absolute neutrino mass scale, mi, there is 
a joint effect of the suppression of the CP asymmetry, more specifically of 
/3max("^i? ""^i) [cf- Eq. (2.48)], plus the extra washout from AL = 2 processes 
mediated by off-shell RH neutrinos, which can be conveniently factored out 
of the efficiency factor in the following way [80] : 

where uj ~ 0.186 and = ml + + m§. This suppression leads to a 
stringent upper bound on the lightest neutrino mass, 

mi<0.12eV, (2.57) 

as derived analytically in [80] and numerically in [75, 88, 99]. The upper bound 
can be seen in the upper part of Fig. 4.3, corresponding to the unfiavored 
regime. It is obtained for high values Mi ~ 10^'^ GcV. 

There have been earlier attempts to get a bound on the absolute neu- 
trino mass scale simply from the fact that mi < mi and then imposing an 
upper limit mi < 10~^ eV, justified as a generic 'out-of-equilibrium' condi- 
tion. However, it turns out that such a restrictive upper limit on m-i does 
not hold for reasons that are clear from Fig. 2.4: the 'out-of-equilibrium' 
condition is realized also in the strong washout regime when z > ~ z^- 
Therefore, within the A^iDS, it is not a problem to have mi as large as 1 eV. 
Incidentally, the upper bound can only be understood when the washout from 
AL — 2 processes and the upper bound on the CP asymmetry are jointly 
taken into account [99]. 

We would like now to summarize the results of this section. As we have 
seen, in the limit where the heavy neutrino mass spectrum is hierarchical 
(HL), Ml <^ M2 <^ M3, the picture of leptogenesis in the unfiavored case sim- 
plifies considerably. A typical scenario of leptogenesis, which we call "vanilla 
leptogenesis" , emerges and it turns out that one has to study only the pro- 
duction of asymmetry by the lightest RH neutrino, Ni, reducing the number 
of parameters from 10 to 6. This is possible because the theoretically favored 
range for Ki lies between K^oi = mgoi/^* = 8.2, given by the solar scale, and 
Xatm = ?7T.atm/"T.* = 48, givcu by the atmospheric scale. In this range, the 
washout from the lightest RH neutrino is strong enough to make the contri- 
butions from the heavier two RH neutrinos negligible, and the final result is 
independent of the initial number of Ni [80, 88, 93]. This is an extremely nice 
feature of the strong washout regime. Note that the A^iDS can actually be 
ensured by setting R23 — 1, which exludes a large contribution from N2 [93]. 
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Within vanilla leptogenesis, it is possible to derive from successful leptoge- 
nesis a lower bound on Mi and on the initial temperature Tin, as shown in 
Eqs. (2.54) and (2.55). Moreover, there is a stringent upper bound on the 
absolute neutrino mass scale, as shown in Eq. (2.57). 

We shall discuss the implications of going beyond the assumptions leading 
to vanilla leptogenesis in Chapter 5. 



2.3 Leptogenesis and supersymmetry 

In the Minimal Supersymmetric Standard Model (MSSM) complemented 
with three RH neutrinos and the corresponding superpartners, the picture 
of leptogenesis is qualitatively quite different from the non-supersymmetric 
case, but it turns out that, quantitatively, they are very similar. 

The interactions of the heavy (s)neutrino field can be derived from the 
leptonic superpotential 

W = ^M,N,N, + K^KH^N, + flKH^E^, (2.58) 

where and Ea arc the SU{2) lepton doublets and singlets chiral super- 
fields, respectively, and Hu and Ha are the Higgs chiral superfields. The 
scalar components of both Higgs bosons, which we denote ^u,d, have vacuum 
expectation values: = Vu — vsm(3 and {^d) = Vd — vcosfi, where 

V — 174: GeV. Their ratio is then given by 

tan/3 = ^. (2.59) 

When flavor effects (see next chapter) are included in supersymmetric lepto- 
genesis, the value of tan (3 is relevant because /^^ = (1 -|- tan^ /3)/J, where fa 
is the SM Yukawa coupling. 

Typically, supersymmetry breaking terms are of no relevance for the 
mechanism of lepton number generation, and we are left with the follow- 
ing trilinear couplings in the Lagrangian, written in terms of four-component 
spinors. 



-hr 



MiLlNi^u + iaPRNi^u + LPR4>"Ni + iiPn^'Ni + /i.e., (2.60) 



where L, and N denote sleptons, higgsinos and singlet sneutrinos, respec- 
tively. 

Prom these couplings one obtains the tree-level relations 

{h^h).. 

^Nii + ^Nil - ^NiL + ^NiLt - ^N*e " ^ NiL " ^i- (2-61) 
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There are now new diagrams contributing to the CP asymmetry. On 
top of the usual contributions shown in Fig. 1.8, there are three additional 
sources coming from the decay of the heavy neutrinos into sleptons, from the 
decay of RH sneutrinos into leptons and from the decay of RH sneutrinos 
into sleptons. One can then define a CP asymmetry for the decay of RH 
neutrinos into leptons and sleptons, and another one for the decay of RH 
sneutrinos into leptons and sleptons, as follows: 

= _ (r^z + r^^)-(r^Lt + rM) ^ ^2.62) 

Tat 

£ff = - , [Z.06) 

N 

where Fjv and denote the total decay rate of RH neutrinos and RH 
sneutrinos, respectively. 

These CP asymmetries were computed in [87] to be 

where we recall that Xj — Mj/M^, and 



g{x) = ^/x 



+ ln 



X — 1 \ X 



X 



a;>l 3 



(2.65) 



In the HL (x^ ^1), the CP asymmetry in the MSSM is therefore twice as 
large as the one in the SM. 

Then, since in the MSSM there are two Higgses, the coefficient Ogph is 
slightly different from its value in the SM, and one obtains for the baryon- 
to-photon ratio 

77b ~ 1.03 X IQ-^SivKi- (2.66) 

As we have seen, there are new decay channels in the MSSM, which yield 
an enhancement of the CP asymmetry by a factor 2. On the other hand, 
there is also an enhancement of the washout by a factor of 2, which implies 
that the constraints on Mi, Tjn and m\ derived in the last section remain 
essentially unchanged [75]. 

Assuming a period of inflation and reheating before leptogenesis occurs, 
the lower bound on the initial temperature of leptogenesis Tin can be iden- 
tified with a lower bound on the reheat temperature Tj-ch of the Universe 
after infiation. Within locally supersymmetric theories, it is well known that 
gravitinos are produced during the reheating phase. The point is actually 
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that they may be overproduced, i.e. their abundance may overdose the Uni- 
verse, leading to the so-called gravitino problem (for early discussions, see 
[100-103]). There are two situations to be distinguished: stable gravitinos 
and unstable ones. 

If the gravitino is the lightest supersymmetric particle, it is stable and 
therefore represents a good dark matter candidate. In order for gravitinos 
not to exceed the dark matter abundance, the reheat temperature has to 
satisfy [104-107] 

Treh < lO'-lO^ GeV. (2.67) 

On the other hand, if gravitinos are unstable, they may lead to a large 
entropy production when they decay during or after big-bang nucleosynthesis, 
spoiling the nice agreement between theory and observations (see Fig. 1.3). 
This leads to the bound (see [108] and references therein) 

T,eh < 10^ GeV, (2.68) 

unless the gravitino mass is larger than about 20 TeV. 

Consequently, whatever specific scenario of supergravity one consid- 
ers, there is a clear tension with the lower bound from leptogenesis given 
in Eq. (2.55). Different ways to relax this tension have been proposed in the 
literature. Let us give three well-known examples. 

One possibility is to produce RH neutrinos non-thermally in the decays of 
the inflaton [75, 109-113]. A recent study [114] shows that the lower bound on 
Treh from leptogenesis can be relaxed in this way by two orders of magnitude. 

Another possibility is provided by "soft leptogenesis" [115,116], which 
is a supersymmetric scenario which requires only one heavy RH neutrino. 
The interference between the CP-odd and CP-even states of the heavy scalar 
neutrino resembles very much the neutral kaon system. The mass splitting 
as well as the required CP violation in the heavy sneutrino system comes 
from the soft supersymmetry breaking A and B terms, associated with the 
Yukawa coupling and mass term of A^i, respectively. The lower bound on the 
reheat temperature in this scenario can go as low as 10^ GeV [75] . 

Finally, it is also possible to use the enhancement of the CP asymmetry 
for quasi-degenerate heavy neutrinos Mi ~ M2 — M3 [87] (see Fig. 2.1) in 
order to relax the lower bound on the reheat temperature. This inspired the 
scenario of "resonant leptogenesis" [76, 117, 118], in which case the scale can 
be lowered to TeV. We shall come back to quasi-degenerate heavy neutrinos 
in Section 5.1. 
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Chapter 3 

Adding flavor to vanilla 
leptogenesis 

We described in the last chapter a picture of leptogenesis where the flavor con- 
tent of the lepton doublets produced in the decays of the heavy neutrinos is 
neglected. This is done by summing over the flavor in the CP asymmetry pa- 
rameter [cf. Eq. (1.24)], in the decay, inverse-decay and scattering rates, and 
the properly normalized number density in the Boltzmann equation (2.28) 
is Nb-l — Ylia^B/z-La- Summing over the flavor implies that all quantities 
in unflavored leptogenesis depend on some combination of {h'^h)ij, where the 
sum over a was explicitly carried out. Thus, even though the Yukawa cou- 
pling in the Lagrangian has a flavor index a [cf. Eq. (B.l)], the latter never 
shows up in the calculation. 

This picture of leptogenesis was thought to be the correct one, or at least 
a good approximation of it, until the beginning of 2006, when two groups 
published independently results that showed the importance of flavor effects 
in leptogenesis [119, 120]. In this chapter we would like to explain in detail 
why flavor effects may be important, when they are expected to matter and 
how they modify in practice the results presented in the previous chapter. 

3.1 When does flavor matter and why? 

The SM Lagrangian contains a term that gives rise to the masses of the 
charged leptons. Before spontaneous symmetry breaking, this term can be 
written as faf-LaCRa^ + h.c., where Crq denote the right-handed lepton flelds, 
which are singlets under SU{2), but with hypercharge -2. Note that we chose 
a lepton basis such that this term is flavor diagonal. 

In the early Universe the Yukawa coupling can be strong enough to 
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maintain processes like la^-a ^ or lofia ^ where A — W^'^, B are 
the SU{2) X f/(l) gauge bosons, in equilibrium, i.e. > H. With a vacuum 
Higgs mass of 120 GeV, the interaction rate was estimated in [121] to be 

~ 5 X 10-' fX (3.1) 

Obviously, when the temperature drops due to the expansion of the Universe, 
the first interactions that will enter equilibrium are the ones involving the 
T-lepton, simply because its Yukawa coupling is larger. It turns out that the 
temperature at which the r-lepton Yukawa interactions enter equilibrium is 
~ 10^^ GeV. For the muon, this will happen at T ~ 10^ GeV, and for the 
electron, at still much lower temperature.^ 

If the charged-lepton Yukawa interactions are in equilibrium (F^ > H) 
and faster than inverse decays (see Section 4.1.2), 

Ta > E Tlr, , (3.2) 

i 

during the relevant period of the asymmetry generation, then the lepton 
quantum states [cf. Eq.(2.1)] lose coherence between the production at 
decay and the subsequent absorption in inverse processes. When the loss of 
coherence is complete, the Higgs bosons will interact with incoherent lepton 
flavor eigenstates instead of the coherent superposition produced in the 
decays. In the limit where the quantum state becomes completely incoherent 
and is fully projected onto one of the flavor eigenstates, each lepton flavor 
can be treated as a statistically independent particle species. This is what we 
call the "fully flavored regime". One has to distinguish a two-flavor regime, 
for 10^ GeV < Mi < 10^^ q^v, such that the condition Eq. (3.2) is satisfied 
only for a = r, and a three-flavor regime, for Mi < 10^ GeV, where the 
condition Eq. (3.2) applies also to a — /j,. 

In the fully flavored regime, there are two new effects compared to the 
unflavored regime [119]. These can be understood introducing the projectors 
and writing them as the sum of two terms, 

Pia ^ mK)\' = PL + ^ (3.3) 
— - - AP° 

The first effect is a reduction of the washout compared to the unflavored 
regime and is described by the tree-level contribution P^^ = {Pia + Pia)/2, 

^In supersymmetric leptogenesis, the rate in Eq. (3.1) is multiplied by (1 +tan/?)^ (see 
Section 2.3), so that flavor effects can start to matter at higher temperatures [122]. 



3.2 Flavored Boltzmann equations and spectator processes 



41 



which sets the fraction of the total asymmetry produced in A^j-decays that 
goes into each single flavor a. In the fully flavored regime, the Higgs will make 
inverse decays on flavor eigenstates instead of the linear superposition 
and hence the washout rate is reduced by the projector P^^. 
The second effect is an additional CP-violating contribution due to a 
different flavor composition between and CP|/^). This can be described in 
terms of the projector differences A PjQ, = Pia — Pia, such that A P^a — 0. 
Indeed, deflning the flavored CP asymmetries, 

r — r 

^ i ' ^ i 

where Via = Pia^i and Via ^ Pia^i and the total decay rates Fj and Fj were 
introduced after Eq. (1.24), these can be now expressed as 

e,^ = SiP^^ + ^^, (3.6) 

showing that the first term is the usual contribution due to a different decay 
rate into leptons and antileptons, and the second is the additional contribu- 
tion due to a possible different flavor composition between and CP\£'j). 
Note that when the flavored CP asymmetries are summed over the flavor, 
one recovers the total CP asymmetry used in the unflavored regime, i.e. 

It is interesting to notice at this point that one can imagine a scenario 
where the total CP asymmetry Si is 0, i.e. no asymmetry would have been 
produced in the unflavored regime, but, thanks to the AP contribution, the 
flavored CP asymmetries Eia do not vanish. More speciflcally, one can imagine 
that the only source of CP violation comes from low-energy phases in the 
PMNS matrix. Since the total CP asymmetry is insensitive to the PMNS 
matrix, £j = in this case. However, APj^ and e^a do explicitly depend on 
the PMNS matrix [see for instance Eq. (3.28)], so that they can be non-zero 
only thanks to the low-energy phases. In great contrast to the unflavored 
picture, fully flavored leptogenesis can then be successful solely thanks to 
the CP violation coming from low-energy phases. This represents a very nice 
possibility which will be the subject of Chapter 6. 



3.2 Flavored Boltzmann equations and spec- 
tator processes 



We follow here the approach presented in [78,119]. Heavy neutrino decays 
produce lepton flavor asymmetries, Nl^, that are computed using Boltz- 
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mann equations similar to Eqs. (2.27) and (2.28). Rigorously, one should 
include the effect of electroweak sphalerons, which constitute an additional 
source of lepton flavor violation. This can be symbolically done by adding 
a washout term dNf^ /dz. Then, for consistency, we also need to add the 
equation dA''^ /dz = dN^^ / dz to account for baryon number violation by the 
sphaleron processes. Given that the sphaleron interactions preserve the three 
charges Aq, = B /3 — La associated to anomaly- free currents, it follows that 
jyEW^g _ jyEW, gy subtracting the equations for the lepton flavor densities 
from the equation for baryon number weighted by a suitable factor 1/3, one 
obtains the following network of flavored Boltzmann equations [119]: 

^ = -D,{N^,-NPJ (3.7) 

^ = eiaD,{N^,-N^^\)-PlWl^{Ne^ + N^), (3.8) 

where, a = e, yU, r and, for simplicity, we included only decays and inverse 
decays with proper subtraction of the resonant contribution from AL = 
2 and AL = processes, and only focused on the lightest RH neutrino 
A^i. Notice that A^a^ is the number density, properly normalized, and 
Nl^ — + , where denote the RH lepton fields. 

Wc are neglecting non-resonant AL = 2 and AL = processes, a good 
approximation for Mi ^ 10^^ GeV, as we will always consider. We are also 
neglecting AL = 1 scatterings, which give a correction to a level less than 
~ 10% in the most interesting strong washout regime [123], and thermal 
corrections. 

The number density asymmetries for the particles X entering in Eq. (3.8) 
are related to the corresponding chemical potentials through 

9xT-^ { l^x/T fermions, 
6 1^ 2jj,x/T bosons, ^ ' 

where gx is the number of degrees of freedom of X. For any given temperature 
regime, the speciflc set of reactions that are in chemical equilibrium enforce 
algebraic relations between different chemical potentials [83]. In the entire 
range of temperatures relevant for leptogenesis, the interactions mediated by 
the top-quark Yukawa coupling ht, and by the gauge interactions, are always 
in equilibrium. Moreover, at the intermediate-low temperatures where flavor 
effects can be important, strong QCD sphalerons [124] are also in equilibrium. 
This situation has the following consequences: 

• Equilibration of the chemical potentials for the different quark colors 
is guaranteed because the chemical potentials of the gluons vanish, 
Hg = 0. 



3.2 Flavored Boltzmann equations and spectator processes 



43 



• Equilibration of the chemical potentials for the two members of a SU (2) 
doublet is guaranteed by the fact that — ~l^w- = above the elec- 
troweak phase transition. This condition was implicitly implemented in 
Eq. (3.8), where we used /i^ = ii^l = I^vl ^'^^ A** = A*</>+ = A*</>o to write 
the particle number asymmetries directly in terms of the number den- 
sities of the SU{2) doublets. 

• Hypercharge neutrality implies 

XI (/^Qi + 2/^«i - l^di - l^ei - l^ei) + 2/x$ = 0, (3.10) 

i 

where Ui, di and ej denote the SU{2) singlet fermions of the i-th gen- 
eration. 

• The equilibration condition for the Yukawa interactions of the top 
quark /it = I^q^ + A«$. 

• Because of their larger rates, QCD sphalerons equilibration occurs at 
higher temperatures than for the corresponding electroweak processes, 
presumably around Tg ~ 10^'^ GcV [125,126], and in any case long 
before equilibrium is reached for the r- Yukawa processes. This implies 
the additional constraint 

X(2/XQ,-/x„,-//,J = 0. (3.11) 

i 

To express the asymmetries A^^^ and N<j, in terms of the Na^, we define 
two matrices, and C*, through the relations [119, 127]: 

Ne. = -Yl C'a/3^^A„ N^--J2 ^|^A„ (3.12) 

13 13 

SO that Eq. (3.8) can be now rewritten as follows: 

^ = D, {N^, - 7V^^) - Wl"" Y.{Ci, + C|)7Va,. (3.13) 

The numerical values of the entries in and C* are determined by the con- 
straints among the various chemical potentials enforced by the fast reactions 
that are in equilibrium in the temperature range {T ~ Mi) where the Aq, 
asymmetries are produced. 

In the temperature range 10^ GeV ^ T < 10^^ GeV, we consider the 
bottom-, tau- and charm- Yukawa interactions to be in equilibrium, implying 
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that the asymmetries in the SU{2) singlets b, Cr and c degrees of freedom 
are populated. The corresponding chemical potentials obey the equilibrium 

constraints fib = /iQg — /ic = /iQj + f^<s> and fir = fii^ — fi^- Moreover, the 
electroweak sphaleron processes are also in equilibrium, implying 



As concerns Icpton number, each electroweak sphaleron transition creates all 
the doublets of the three generations, implying that individual lepton fla- 
vor numbers are no longer conserved. An asymmetry will then be generated 
along the £r and ig+n directions in flavor space. Even though the electroweak 
sphalerons induce L± ^ 0, where by the subscript _L we mean the direction 
in flavor space perpendicular to i^- and ie+ij., the condition A_|_ = is not vio- 
lated, and hence Eq. (3.13) consists of just two equations for N^^ and N^^^^. 
This is what we call the two-flavor regime. As concerns baryon number, elec- 
troweak sphalerons are the only source of B violation, implying that baryon 
number is equally distributed among the three quark generations, i.e. B/3 
in each generation. This modifles the detailed equilibrium conditions for the 
quark chemical potentials. Solving the corresponding set of linear equations 
for the chemical potentials, one obtains [128] 



In the temperature range T < 10^ GeV, one has to include the equih- 
bration constraints from the strange-quark Yukawa interactions and, more 
importantly, from the muon- Yukawa interactions. Given that the electron 
remains the only lepton with a negligible Yukawa coupling, the Yukawa in- 
teractions completely deflne the flavor basis for the leptons as well as for the 
antileptons (that are now the CP-conjugate states of the leptons). Corre- 
spondingly, the lepton asymmetries are also completely defined in the flavor 
basis. In this regime the coefficients C^^ and C* are given by [119, 129] 




(3.14) 





(3.15) 




(3.16) 



3.3 In practice, what changes? 



We have seen in the last section that, if the relevant temperature for lep- 
togenesis is below roughly 10^^ GeV, a rigorous description of the asymme- 
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try evolution has to be performed in terms of the individual flavor asym- 
metries Aq = B/3 — La rather than in terms of the total asymmetry 
Nb-l = ^a„! cis usually done in the unflavored treatment. 
We would like to make two remarks: 

• From Eq. (3.13) it can be seen that the two diagonal entries of 
sum up with the entries in C*. Actually, this sum turns out to be close 
to 1 both using Eq. (3.15) and Eq. (3.16); moreover, the addition of the 
coefficients in C* tends to equahze the flavored diagonal elements [128] . 

• It was shown in [130] that the off-diagonal elements in the matrix 
have an effect smaller than 40% on the final asymmetry. Taking into 
account the matrix C*, this effect is even reduced. 

Therefore, for simplicity, we shall use in the following + C* = 1, 
which takes into account spectator processes in an approximate way. The 
generalization of the flavored Eqs. (3.7) and (3.8) to three RH neutrinos 
(i = 1, 2, 3) is then given by [119, 120, 131, 132] 

^ = - A (A^iv, - A^^^) (3.17) 
^ = E^-A(iV^,-iV^^)-E^«^/°^A„. (3.18) 

i i 

Notice that, in the two-flavor regime, Aq, = B /2 — La, where a = e + n oi 
T, and the two equations for the individual electron and muon asymmetries 
are replaced by one kinetic equation for the sum Na^+^, where the individual 
flavored CP asymmetries and projectors have to be replaced by their sum, 
namely eie+^j = £ifj.+£ie and Pie+^i = Pi^ + Pie [129]. The total asymmetry is 
then given by Nb-l = -^Ae+^+-^Ar- The calculation is therefore intermediate 
between the unflavored case and the three-flavor regime, though the results 
are very similar to the three-flavor regime [123] . 

The evolution of the N^^s can be worked out in an integral form, 

NAaiz) = N'^a exp ^ia [ ^l"" i^')^ + E '^-(^) ' (^-W) 

with the 6, in the two-flavor case, or 9, in the three-flavor case, efficiency 
factors given by 



{z; K,, PL) = - [ d.' ^ exp (- E PL £ d." Wr{z")\ . (3.20) 
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The final B — L asymmetry is then given by 

A^i-L = E<.' (3-21) 

from which one obtains the baryon-to-photon ratio rjs using Eq. (2.15). 

Using the convenient Casas- Ibarra parametrization, Eq. (2.18), the tree- 
level projectors can be written as 

It will also prove useful to introduce the fiavored decay parameters, 

2 

(3.23) 



X = P .° K 




obtained using Eqs. (3.22) and (2.23). 

The fiavored CP asymmetries are given by the following expression [87] : 



2 



Xi 



3(x,M-l) ^"' [^''c.Mh'^h),,] I , (3.24) 



where ^{x) was defined in Eq. (2.17), and we recall that Xi = Mf /Ml. 

In general, the final asymmetry will depend on all 18 see-saw pa- 
rameters. As explained in the introduction (Section 1.2.2), until now we 
have only measured two mass-squared difi^erences and two mixing angles 
in neutrino oscillation experiments. This is essentially the only informa- 
tion on the 18 see-saw parameters we have. Thus, we can write that 
'Hb — VBiiTii, U, Mi, a;2i, wsi, a;32). It is interesting that including fiavor effects 
there is a potential dependence of the final asymmetry also on the unknown 
parameters contained in the PMNS mixing matrix U [119, 131], namely the 
mixing angle 9is and the CP- violating phases 6 and $1^2 [cf. Eq. (A. 5)]. 

Let us assume that the mass spectrum of the heavy neutrinos is hierar- 
chical. Ml <^ M2 <^ M3. In this case the general expression Eq. (3.24) for 
the CP asymmetries reduces to 



3.3 In practice, what changes? 



47 



-hm[hl,K^{h^h)u]Y (3.26) 
^3a - - 8^(^\/,)33 E^"^ [K3hajih^h)js\ . (3.27) 

Expressing the flavored CP asymmetries eia in terms of the orthogonal 
parametrization, one obtains [129] 

3Mi sr^ mi ^Jrn^n^, 
£ia = - E ^ ^^[^"'^ ^'i] ' (3.28) 

where the explicit dependence on the PMNS matrix can be seen. Similar 
expressions can be found for £2a and £30,, but as they are slightly more com- 
plicated, we do not show them here. 

In the following, we shall assume no rotation in the plane 23, i.e. i?23 = 1 
[cf. Eq. (2.22)]. Under these conditions, both total CP asymmetries £2 and 
£3 are suppressed like ~ Mi/M2,3. On the other hand, it is interesting to 
notice that the £20 's and the £30 's are not necessarily suppressed. This can 
potentially lead to a scenario where the final asymmetry is produced by the 
decays of the two heavier RH neutrinos, provided M2,3 < 10^^ GeV in order 
for the flavored regime to apply. Here we do not pursue this possibility (see, 
however. Sections 5.2 and 6.2) and focus on a typical A'"i-dominated scenario 
where the dominant contribution to the final asymmetry comes from the 
decays of the hghtest RH neutrino, so that 

N'b-l - N'b_l\n, = E '^i" • (3.29) 

a 

It will prove important for our discussion that both the total CP asym- 
metry £1 and the flavored ones Sia cannot be arbitrarily large. The total 
CP asymmetry is indeed upper bounded by Eq. (2.45), and each flavored 
CP asymmetry Sia is bounded by [129] 

|£i«| < £(Mi) y^^maxj [|t/„,|] , (3.30) 

where we recall that e{Mi) = 3 Mi rriatm/ (167r v^). Therefore, while the total 
CP asymmetry is suppressed when mi increases, the single-flavor CP asym- 
metries can be enhanced. The existence of an upper bound on the quantity 

ria = eiJe{Mi) (3.31) 
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independent of Mi, implies, as in the unflavored analysis, Eq. (2.52), the 
existence of a lower bound on Mi given by 

Ml > Mf^'iKi) = -r-, — — - , (3.32) 

where we will always use the 3(J lower value for Mi [cf. Eq. (2.53)], and 
we indicated with k\{Ki) the efficiency factor in the unflavored case, corre- 
sponding to with P^^ = 1. We also defined 

Ci^E^i- ^ith ^,^^!2^l4r^- (3-33) 

a J. \ / 

This quantity represents the deviation introduced by fiavor effects compared 
to the unflavored treatment in the hierarchical fight neutrino case. We could 
then use ria < ^/Pi^rn^/msx.ra to maximize ^i. Notice however that, first, the 
ric's cannot be simultaneously equal to \/^^ because of the bound on the 
total asymmetry, and, second, there are sign cancellations in Therefore, 
the bound (3.32) is more restrictive than this possible estimation, and we 
prefer to keep it in this form, maximizing in each particular situation. 

As usual, the lower bound on Mi implies an associated lower bound on the 
initial temperature Tin of leptogenesis and hence on the reheat temperature 
T 

reh- 



3.4 Dependence on the initial conditions and 
lower bounds 

From Eq. (3.20), extending an analytic procedure derived within the unfla- 
vored treatment [80], one can obtain simple expressions for the flavored effi- 
ciency factors K^Q,. In the case of a thermal initial iVi-abundance (A^jvi = 1); 
one has 

~ k{Ki^) , (3.34) 

where the function n{x) was given in Eq. (2.31). Notice that, in the particu- 
larly relevant range 5 < Ki^ < 100, this expression is well approximated by 
Eq. (2.35) replacing Ki — > Ki^. 

In the case of vanishing initial A^i-abundance {N'^^ = 0), one has to take 
into account two different contributions, a negative and a positive one, 

«L = ^-{Ku PL) + 4(i^i, PL) . (3.35) 
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The negative contribution arises from a first stage when N^^ < N^^, for 
z < Zeq, and is given approximately by 



exp 



3 TT Kia 



exp(^7V^,(^eq)) -1 



, (3.36) 



where A^Ar^(-2cq) was defined in Eq. (2.38). 

The positive contribution arises from a second stage when N^-^ > N^^, 
for z > Zeq, and is approximately given by 



exp 



Zeq, 



. (3.37) 



It is interesting to notice that N]\f-^^{zeq) is still regulated by Ki [cf. Eq. (2.38)], 
since the RH neutrino production is not affected by fiavor effects, contrarily 
to the washout, which is reduced because regulated by Kia. 

These analytic expressions make transparent the two conditions to have 
independence from the initial conditions. The first is the thermalization of 
the A^i-abundance, such that, for an arbitrary initial A^i-abundance, one 
has Nf^^{zeq) — 1. The second is that the asymmetry produced during the 
non-thermal stage, for z < Zeq, has to be efficiently washed out, leading to 
<^ They are both realized for large values Ki ^ 1. More quantita- 
tively, we introduced in Section 2.2 the quantity K^, such that, for Ki > K^, 
the final asymmetry calculated for a thermal initial A^i-abundance differs 
from the one calculated for a vanishing initial abundance by less than some 
quantity S. This can be used as a precise definition of the strong washout 
regime. Let us consider some particular cases, showing how flavor effects 
tend to enlarge the domain of the weak washout at the expense of the strong 
washout regime. 



3.4.1 Alignment 

The simplest situation is the alignment case, realized when the A^i-decays 
are just into one flavor a, so that Pia = Pia = 1 and Pip^a = Pip^a = 0, 
implying eia = Si. Notice that we do not have to worry about the fact 
that the lightest RH neutrino inverse decays might not be able to wash 
out the asymmetry generated from the decays of the two heavier neutrinos, 
since we are assuming negligible 52/3 and £3/3 anyway. In this case the general 
set of kinetic equations (3.17) and (3.18) reduces to the usual unflavorcd 
equations and all results coincide with those in the unflavored analysis [119]. 
In particular one has N^b-l — ^1 'Aa- 



50 



Chap. 3: Adding flavor to vanilla leptogenesis 




Figure 3.1: Efficiency factors (upper panel) and lower bounds on Mi (lower 
panel) in the alignment (solid lines), semi-democratic (dashed lines) and 
democratic (short-dashed lines) cases. 



In the case of alignment we are considering, we obtain ~ 3.5 for 
S — 0.1, as shown in Fig. 3.1. The value of plays a relevant role since 
only for Ki > one has predictions from leptogenesis on the final baryon 
asymmetry resulting from a self-contained set of assumptions. On the other 
hand, for Ki < K^, leptogenesis has to be complemented with a model for 
the initial conditions. Additionally, the calculation of the final asymmetry 
in the weak washout regime requires a precise description of the RH neu- 
trino production, potentially sensitive to many poorly known effects. It is 
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then interesting that current neutrino mixing data favor Ki to be in the 
range Kgoi — 8.2 ^ -f^i ^ 48 ~ -ft'atm [80,88,93], where one can have a 
mild washout assuring full independence from the initial conditions, as one 
can see in Fig. 3.1, but still successful leptogenesis. In this case one can 
place constraints on the see-saw parameters which do not depend on specific 
assumptions for the initial conditions and with reduced theoretical uncer- 
tainties. 

Since one has = 1 here [cf. Eq. (3.33)], the general lower bound on Mi 
in Eq. (3.32), like all other quantities, becomes the usual lower bound derived 
from an unflavored treatment. In the lower panel of Fig. 3.1, we have plotted 
it both for N^^_^ — 1 and N^^^ — (solid lines). This corresponds precisely 
to the lower bound in the unflavored treatment shown in Fig. 2.5. One can 
see that the dependence on the initial conditions in k,\^ translates into a 
dependence on the initial conditions in M™™. The lowest model-independent 
values are then obtained for = i^^* ~ 3.5 and are given by 

Ml > 4 X 10^ GeV and T.^h > 1-5 x 10^ GeV . (3.38) 

We did not show the lower bound on T^ch in Fig- 3.1, not to overload the plot, 
but it is precisely given by the dashed line in Fig. 2.5. Another typically 
quoted lower bound on Mi is the one obtained for a thermal initial Ni- 
abundance in the limit Ki 0, given by Mi > 5 x 10^ GeV [81]. 

3.4.2 Democratic and semi-democratic cases 

Let us now discuss another possibility. For definiteness, we assume a three- 
flavor regime; the extension of the results to the two-flavor regime is straight- 
forward. Let us assume a democratic situation where Pia = Pia = 1/3 for 
any a and consequently APi^ = 0. This case was also considered in [120]. 
From Eq. (3.20) it follows that the three efficiency factors k\^, like the 
three CP asymmetries Sia, are all equal and thus Eq. (3.29) simplifies into 
N^B-L — ^1 'Aa-i ^ i^ usual unflavored treatment. However, the washout is 
now reduced by the presence of the projector, so that Ki Ki = Ki/3. 

The result is that, in the case of a thermal initial A'^i -abundance, the effi- 
ciency factor, as a function of Ki, is simply shifted. The same happens for 
vanishing initial A^i-abundance in the strong washout regime. However, in 
the weak washout regime, there is not only a simple shift, since the RH 
neutrino production is still depending on i^i. A plot of is shown in the 
upper panel of Fig. 3.1 (short-dashed lines). One can see how the reduced 
washout increases the value of to ~ 10, approximately I /Pia — 3 larger, 
thus compensating almost exactly the washout reduction by a factor ~ 3^'^ 
[cf. Eq. (2.35)]. 
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In this way the lowest bound on Mi in the strong washout regime, at Ki — 
K^, is almost unchanged. On the other hand, for a given value Ki ^ K^^, the 
lower bound gets approximately relaxed by a factor 3 [120]. The lower bound 
for the democratic case is shown in the lower panel of Fig. 3.1 (short-dashed 
lines). It is apparent that the lower bound for — > and thermal initial 
A^i-abundance does not change with respect to the alignment case (or the 
unflavored treatment). One can then say that flavor effects simply induce a 
shift of the dependence of the lower bound on Ki. 

The semi-democratic case is intermediate between the democratic and 
the alignment cases. It is obtained when one projector vanishes, for example 
Pip — 0, and the other two are 1/2. In this case, K^, ~ 7. The correspond- 
ing plots of the efficiency factor and of the lower bound on Mi are also 
shown in Fig. 3.1 (dashed lines). The semi-democratic case can actually be 
identified with the two-flavor regime, where the two projectors are equal. 



3.4.3 One-flavor dominance 

There is another potentially interesting situation that motivates an extension 
of the previous results to arbitrarily small values of Pi^. This occurs when 
the final asymmetry is dominated by one fiavor a, and Eq. (3.29) can be 
further simplified into 



analogously to the alignment case but with <ti 1. Notice that this cannot 
happen due to a dominance of one of the CP asymmetries, for example with 
Eia being close to its maximum value, Eq. (3.30), much larger than the other 
two that are strongly suppressed, simply because one has APi^ = 0. One 
has then to imagine a situation where the CP asymmetry Sia is comparable 
to the sum of the other two, but Kip^a ^ ^la ^ 1, so that ^i^^ kI^. The 
dominance is then a result of the much weaker washout. 

The analysis of the dependence on the initial conditions can then be 
performed as in the previous cases calculating the value of for any value 
of Pi°Q,. The result is shown in Fig. 3.2. The alignment case corresponds to 
P^^ = 1, the semi-democratic case to P^^ — 1/2 and the democratic case 
to Pi°„ = 1/3. Notice that the resuh is very close to the simple estimation 
K^,{Pl^) = K^{l)/P^^ which would follow if k.{^ were just depending on Kia- 
In Fig. 3.3 we have plotted the values of the lower bounds on Mi and Tj-eh 
for hierarchical light neutrinos, implying ms = matm in Eq. (3.30). These 



can be obtained plugging ^i = ^/P^ K{^{Kia) / k{{Ki) < 1 in Eq. (3.32). 



They correspond to the lowest values in the strong washout regime, when 




1/2. 



B-L — ^la l^la ) 



(3.39) 
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Figure 3.2: Values of defining the strong washout regime, for 5 = 10% 
(solid line) and 6 = 50% (dashed line). 



There are two possible ways to look at the results. On the one hand, 
flavor effects can relax the lower bounds for fixed values of i^i ^ 1. Indeed, 
for each value i^i ^ 1, one can choose P^^ = K^{P^^ = l)/Ki such that 
Ki = K^,{P^^), thus obtaining the highest possible relaxation in the strong 
washout regime. This is shown in the right panel of Fig. 3.3. It is important 
to say that this relaxation is potential. A direct inspection is indeed necessary 
to determine whether it is really possible to achieve at the same time not 
only small values of P^^ but also a single-flavor CP asymmetry eia that is 
not suppressed compared to ei^^a- 

On the other hand, as a function of P^^ the bounds get more stringent 
when P^^ decreases, so that the minimum is obtained in the alignment case, 
corresponding to the unfiavored case. This is clearly visible in the left panel 
of Fig. 3.3. Therefore, it is important to emphasize here that flavor effects 
cannot help to alleviate the conflict of the lower bound on the reheat tem- 
perature Treh from successful leptogenesis with the upper bound on Tj-eh in 
order not to overproduce gravitinos (see Section 2.3 for a discussion). In par- 
ticular, the bounds on Mi and Tj-eh that are usually quoted in the literature 
[cf. Eq. (3.38)] are not changed by flavor effects. 

Notice that together with the one-flavor dominance case, one can also 
envisage, in the three-flavor regime, a two-flavor dominance case, where two 
projectors are equally small and the third is necessarily close to one, while 
all the three flavored CP asymmetries are comparable. 

In the next section we consider a specific example that illustrates what 
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Fi gure 3.3: Lower bounds on Mi and Treh calculated choosing P^^ — 
Ki,{Pi^ = l)/Ki such that Ki = i^'^(P°„) (thick lines) and compared with 
the usual bounds for P^^ = 1 (thin lines). In the left panel they are plotted 
as a function of Pi^, while in the right panel as a function of Ki. 

we have discussed on general grounds. At the same time, it will help to 
understand which are realistic values for the projectors and their differences, 
given a specific set of see-saw parameters and using the information on the 
PMNS mixing matrix we have from neutrino oscillation experiments. 

3.5 Study of a specific example 

The previous results have been obtained assuming no restrictions on the 
projectors. Moreover, in the one-flavor dominance case, where there can be a 
relevant relaxation of the usual lower bounds derived following an unflavored 
treatment, we have assumed that the upper bound on eia, Eq. (3.30), is 
saturated independently of the value of the projector. 

This assumption does not take into account that the values of the projec- 
tors depend on the different see-saw parameters, in particular on the neutrino 
mixing parameters, and that severe restrictions could apply. Let us show 
a definite example considering a particular form of the orthogonal matrix, 
Q = i?i3 [cf. Eq. (2.21)]. This case is particularly meaningful, since it realizes 
one of the conditions (^21 = 0) to saturate the bound (2.45) for ei. Moreover, 
the decay parameter is given by 

Ki = K^,,, \l- colli + i^atm ^31 1 > (3-40) 
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where Xmin = iTii/m^, and we recall that i^^atm = tTiatm/iTii,- The expres- 
sion (3.22) for the projector gets then specialized as 
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while, specializing Eq. (3.28) for eia, one obtains 
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where m^/m^tm = a/1 + mf/m^ 

If we first consider the case of fully hierarchical light neutrinos, mi 
then 
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and 
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0. 



(3.43) 



For the PMNS matrix U we adopt the parametrization Eq. (A. 5). One then 
finds 

P{',<0.03, P{;^P°,^ 1/2. (3.44) 

The fact that the projector on the electron fiavor is very small and AP°^ = 
for all fiavors implies that the asymmetry generated in the electron fiavor 
will be small, while the muon and tauon contributions will be equal, since 
the projectors on these flavors are equal. Summing Eq. (3.18) over a, one 
then simply recovers the unflavored case with the washout reduced by a 
factor of 2. This is a realization of the semi-democratic case that we were 
envisaging at the end of Section 3.4.2 where ~ 7. In this situation fiavor 
effects do not produce large modifications to the usual results, essentially 
a factor of 2 reduction of the washout in the strong washout regime with a 
consequent equal relaxation of the lower bounds (see dashed lines in Fig. 3.1). 
Moreover, there is practically no difference between a calculation in the two- 
or three-flavor regime. 

Let us now consider the effect of a non- vanishing but small lightest neu- 
trino mass mi, for example mi = O.lmatm- In this case the results can also 
depend on the Majorana and Dirac phases. We will show the results for uj^^ 
purely imaginary, the second condition that maximizes the total CP asym- 
metry if mi <^ ?7T.atm [93]. Noticc that this is not in general the condition 
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Figure 3.4: Dependence of different quantities on Ki for mi/matm = 0.1 
and real U. Left panel: projectors P^^ and normalized CP asymmetries ri^; 
central panel: and as defined in Eq. (3.33) for thermal (thin) and 
vanishing (thick) initial A'"i-abundances; right panel: lower bound on Mi for 
thermal (thin solid) and vanishing (thick solid) abundances compared with 
the unflavored result (dash-dotted line). 



that maximizes ria for mi > matm; however, we will also use it in this case 
for simplicity. 

We first consider the case of a real U. The results are only shghtly sensi- 
tive to a variation of ^13 within the experimentally allowed 3a range 0-0.2. 
Therefore, we shall set ^13 = 0, corresponding to Ue3 = 0, in all examples. 
In the left panel of Fig. 3.4 we show the values of the projectors P{J^ and of 
the normalized CP asymmetries function of Ki . The calculations are 

performed in the two-flavor regime since we obtain that successful leptogen- 
esis is possible only for Mi > 10^ GeV, where the two-flavor regime apphes. 
Now a difference between the tauon and the sum of the muon and electron 
projectors and asymmetries arises. On the other hand, for Ki ^ 100, this 
difference tends to vanish and the semi-democratic case is again recovered. 

In the central panel we show the quantities [cf. Eq. (3.33)] and their 
sum ^1. We recall that ^1 gives the deviation of the total asymmetry from 
an unflavored calculation for hierarchical light neutrinos. One can sec how 
the contribution to the total asymmetry from the tauon flavor is now twice 
as large as from the electron plus muon flavors. For Ki ^ 100, the semi- 
democratic case is restored, the two contributions tend to be equal to the 
unflavored case, and the total flnal asymmetry is about twice larger. Finally, 
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Figure 3.5: Same quantities as in the previous figure but with one non- 
vanishing Majorana phase: $i = — tt. 



in the right panel we show the lower bound on Mi, and compare it with 
the results from the unfiavored analysis (dash-dotted line). At ~ 14, the 
relaxation in the strong washout is maximum, a factor ~ 3. For Ki ^ K^, 
the relaxation is reduced to a factor 2, as in the semi-democratic case. 

Let us now study the effect of switching on phases in the U matrix, again 
for mi = O.lmatm- The most important effect arises from one of the two 
Majorana phases, namely $i. The other Majorana phase $2 is irrelevant in 
the particular model we consider, Q = Ri^. As for the Dirac phase 5, it plays 
a role similar to $1, but its effect is suppressed by the small angle ^13. In 
Fig. 3.5 we show, again in three panels, the same quantities as in Fig. 3.4 for 
$1 = — TT. One can see how this further increases the difference between the 
e + ii and the r contributions and further relaxes the lower bound on Mi . The 
effect is small for the considered value mi/matm = 0.1. However, considering 
a much larger neutrino mass mi while keeping $1 = — tt, the effect becomes 
dramatically bigger. In Fig. 3.6 we show the same quantities as in Fig. 3.4 
and Fig. 3.5 for mi/matm = 10. In the left panel one can see that now, for 
Ki -ftTmin, kie+/i| — I^itI — 2, much larger than in the previous case. This 
means that the dominant contribution to the fiavored CP asymmetries comes 
now from the APi^ term [cf. Eq. (3.6)]. At the same time, very importantly, 
Pi^ <^ Pie+^ and in this way, as one can see in the central panel, the dominant 
contribution to .^i is given by ^i^-. This case thus finally realizes a one-flavor 
dominance. The final effect is that the lower bound on Mi is about three 
orders of magnitude relaxed compared to the unfiavored case. 
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Figure 3.6: Same quantities as in the previous two figures but with 
TTii/mgjtm — 10 and one non- vanishing Majorana phase: $i = — tt. 

In Fig. 3.7 the dependence on mi/matm of the lower bound on Mi is 
summarized, showing both the case with zero Majorana phase and the case 
with $1 = — TT. One can notice how the effect of the phase in relaxing the 
lower bound increases with mi/m^t^. 

Finally, we want to study the interesting case of a real orthogonal matrix 
Jl, implying Si = 0. In the particular model we are considering, namely 
= it!i3, there is no asymmetry produced if mi = 0, since £ia oc £i = [cf. 
Eq. (3.42)]. For a non- vanishing mi and a non-real U, we have that APi^ 7^ 
and consequently eia 7^ 0. In Fig. 3.8 we show the results for mi/matm = 0.1 
and $1 = 7t/2. Wc present again the same quantities as in Figs. 3.4, 3.5 
and 3.6. One can see that an asymmetry can still be produced, as envisaged 
in [119]. However, successful leptogenesis is possible almost only in the weak 
washout regime. In the strong washout regime, for values Mi < 10^^ GeV, 
there is a small allowed region only for K^^ — 14 < i^^i < 30. 

Let us summarize the main results of this section, where we assumed the 
Q matrix to have the simple form Q = Ri^. 

We have seen that the relaxation of the lower bounds compared to the 
unfiavored case is only of order one for mi <S matm, but grows to several 
orders of magnitude for mi > 10 matm — 0.5 cV, i.e. in the quasi-degenerate 
limit. The reason is that for quasi-dcgcncratc neutrinos (mi ~ m2 ~ m^) a 
partial cancellation in one of the projectors, Eq. (3.22), can be more easily 
obtained, simply because all terms in the numerator arc of the same order. 
When a partial cancellation occurs, one projector will be suppressed, leading 
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Figure 3.7: Lower bound on Mi. The solid lines arc for vanishing phase, while 
the dashed lines are for $1 = — tt. The results for both vanishing (thick lines) 
and thermal (thin lines) initial abundances are presented. 




Figure 3.8: Same quantities as in Figs. 3.4, 3.5 and 3.6 in the case of real 
Q for mi/matm = 0.1 and $1 = 7r/2. The dot-dashed fine still refers to the 
unfiavored case for = it!i3 and purely imaginary 00^^. 
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to a one-flavor dominance and thus to a large relaxation of the lower bound 
at large values of Ki. This is exactly what happened in the situation depicted 
in Fig. 3.6. 

The role of the phases in the U matrix to achieve a cancellation in one 
of the projectors is crucial. For Q = R13, only one Majorana phase, $1, is 
relevant, and the specific value it takes leads to very different results. In 
Fig. 3.6 we used the value $1 = — tt which approximately maximizes the 
asymmetry. The Dirac phase S can also play a similar role to $1 for Q = R13, 
but its effect is always suppressed by the small angle 6*13. It should be noted 
that for more general cases like Q = i?i2-Ri3, the cancellation in the projector 
leading to a one-flavor dominance is also possible when mi = 0, because 
there is more freedom in the fl matrix to achieve a cancellation. Moreover, 
both Majorana phases in U can play a role in principle, not only $1. 

It seems therefore that when flavor effects are taken into account, one 
obtains an opposite result compared to the unfiavored analysis we performed 
in the previous chapter, where one had a suppression of the final asym- 
metry for growing absolute neutrino mass scale, leading to the stringent 
upper bound Eq. (2.57). This upper bound seems now to hold only for 
Ml > 10^^ GeV [120]. However, as we shall argue in the next chapter, this 
issue requires a full quantum kinetic calculation, in order to keep track of the 
correlations in flavor space and of partial losses of coherence [133]. 

We wish also to emphasize that, when flavor effects arc included, leptoge- 
nesis is an interesting example of phenomenology, beyond neutrinoless double 
beta decay, where Majorana phases play an important role. The Dirac phase 
is also relevant, although it appears in combination with the small angle ^13. 
The importance of the CP-violating phases in U shows up at two levels. First, 
as we have discussed above, they arc crucial to achieve a cancellation in one 
of the projectors, leading to a one-flavor dominance. Second, they can pro- 
vide the unique source of CP violation required for successful leptogenesis, 
as illustrated in Fig. 3.8. We shall come back to this interesting possibihty 
in Chapter 6. 



Chapter 4 



Prom classical to quantum 
kinetic equations 

We saw in the last two chapters that two quite different pictures of Icpto- 
genesis are vahd in two different temperature regimes. Roughly speaking, 
when the temperature relevant for leptogenesis is above 10^^ GeV, the un- 
flavored picture should hold, whereas at lower temperatures, T < 10^^ GeV, 
the fully flavored regime should apply, cither with two independent flavors 
(10^ GeV < T < 10^2 GeV) or with three (T < 10^ GcV). 

The final baryon asymmetry was estimated in both regimes, unfiavored 
and fully flavored, by means of classical Boltzmann equations. However, one 
could expect that quantum effects, such as correlations in flavor space and 
partial losses of coherence, might play some role in the intermediate regime 
around 10^^ GeV. In such a case, the classical Boltzmann equation for the 
Aq asymmetry given in Eq. (3.18) is not expected to describe correctly the 
generation of asymmetry. One should indeed turn to a density matrix equa- 
tion, which provides the right formalism to keep track of correlations in flavor 
space and partial losses of coherence. 

In this chapter we re-visit the conditions of validity of both the unfla- 
vored and the fully flavored pictures. We shall give physical arguments why 
one actually expects to have a substantial part of the parameter space where 
quantum effects are important. Interestingly, this part of the parameter space 
is exactly the one relevant for a discussion about the upper bound on the 
absolute neutrino mass scale, which is evaded if one believes a fully flavored 
analysis. We shall argue that only solving the relevant density matrix equa- 
tion will give a deflnite answer to that question. In the flnal section we discuss 
the density matrix formalism and give tentative equations for leptogenesis in 
the transition region. 



62 



Chap. 4: From classical to quantum kinetic equations 



4.1 Validity of the different pictures 

Let us first discuss in some detail tlie quantity that will be the angular 
stone of the subsequent analysis, namely the washout rate by inverse decays 
[of. Eq. (2.12)]. As usual, we restrict our analysis to the lightest RH neutrino 



It was shown in [80] that the inverse-decay rate reaches a maximum 
Wl^{zjaax) — 0.3 Ki at Zmax — 2.4. In the weak washout regime, when 
-^1 ^ 3.3, one has Wl^{z) < 1 for any value of z. In this case the washout is 
negligible, and the final asymmetry depends on the initial conditions. In the 
strong washout regime, when Ki > 3.3, there is an interval [zon, Zos] where 
Wl^ > 1. The asymmetry produced at z ^ Zos is very efficiently washed out, 
and thus the final asymmetry is essentially what is produced around zb — Zos 
by the out-of-equilibrium decays of the residual RH neutrinos, whose number 
corresponds approximately to the final value of the efficiency factor. 

In Fig. 4.1 we show the washout term from inverse decays for three dif- 
ferent values of Ki. For Ki = 100, we show the interval [zon, Zofi]. The max- 
imum value, Wl^{zma.x) — 33, is reached at ^max — 2.4. For Ki ~ 3.3, one 
has Zon — ^max — Zqs- This can be taken as the threshold value distinguish- 
ing between the strong and the weak washout regimes. For Ki — 10~^, one 
has Wl^ <^ 1 for any value of z. Notice however that even in this case the 
weak washout can be important for successful leptogenesis if the initial iVi- 
abundance is zero, since it prevents a full cancellation between two different 
sign contributions to the final asymmetry [80]. 

4.1.1 When are flavor efl'ects important? 

As discussed in Section 3.1, flavor effects are caused by charged-lepton 
Yukawa interactions [127] that occur with the rate given in Eq. (3.1). The 
largest rate is for a = r where 



Therefore, if T > 10^^ GeV, charged-lepton Yukawa interactions are not ef- 
fective, and all processes in the early Universe are flavor blind, justifying the 
unflavored treatment. For T < 10^^ GcV. the r- Yukawa coupling is strong 
enough that the scatterings £r Cr ^ are in equilibrium. However, this 
condition is not necessarily sufficient for important flavor effects to occur, 
because we need to compare the speed of the Yukawa interactions with that 
of the RH neutrino decays and inverse decays. To this end, we study the 




(4.1) 
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Figure 4.1: Comparison between the washout term Wl^{z) (thick sohd hues), 
defined in Eq. (2.12) and plotted for the three indicated values of Ki, and the 
charged-lepton Yukawa interaction term F^, defined in Eq. (4.5) and plotted 
for the indicated values of Mi. 



weak and strong washout regimes separately and consider only a two-flavor 
case, because the r-lepton Yukawa coupling causes the main modification. 

In the weak washout regime, assuming a vanishing initial A^i-abundance, 
the production of RH neutrinos by inverse decays occurs around T ~ Mi. 
At this epoch, inverse decays are by definition slower than the expansion 
rate. Therefore, the condition T < 10^^ GeV is sufficient to conclude that the 
charged-lepton Yukawa interactions are faster than the inverse decay rate. 
This translates into the condition Mi < 10^^ GeV because the RH neutrino 
production occurs at T ~ Mi, in agreement with the previous literature [119, 
120]. 

However, this condition does not guarantee that fiavor effects indeed have 
an impact on the final asymmetry, because this impact depends on washout 
playing some role. For a vanishing initial A^i-abundance, this is the case in 
that washout effects prevent a full sign cancellation between the asymmetry 
produced when Nj^^ < iV^ and the asymmetry produced later on. On the 
other hand, for a thermal initial iVi-abundance, no such effect arises from 
the weak washout and fiavor effects do not modify the final asymmetry. We 
will come back to this point later on. 
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In the strong washout regime, the situation is very different. The rate of 
RH neutrino inverse decays at T ~ Mi is larger than the expansion rate. 
Therefore, we need to compare the charged-lepton Yukawa rate with the 
RH neutrino inverse-decay rate F^^. For the unflavored treatment to be vahd 
for z < Zfi < zb then requires 

^ 10^2 GeV 

where Zfi is that value of z where the two rates are equal, i.e. F^°(^fl) = Tr{zf[). 
This condition guarantees that at temperatures T > = Mi/zfi flavor 
effects will not be able to break the coherent propagation of the Icpton states. 
The final asymmetry is dominantly produced around z ^ zb- Therefore, the 
condition for flavor effects to be negligible is 

Ml > 5 X 10^^ GeV , (4.3) 

similar to the weak washout regime. However, the corresponding condition 
on the temperature 

_ ^ 10^2 GeV 

T > (4.4) 

is now less restrictive. 

If one starts with a non-vanishing initial A'^i-abundance, then the flnal 
asymmetry is also determined by how efficiently the initial value is washed 
out; this is described by the integral in Eq. (2.29). In this case even a value of 
Zfi < zb would be important to determine the flnal asymmetry since washout 
in the unflavored regime would be effective for z < Zfl < Zb, while a reduced 
washout would apply in the ffavored regime at lower temperatures such that 

We conclude that the condition (4.2) obeys the intuitive expectation that 
there is always a threshold value for Ki above which the unflavored case is 
recovered. In this case the temperature below which flavor effects play a role 
indeed becomes smaller and smaller. The situation is illustrated in Fig. 4.1 
where we compare Wid with 

^ _ 1 F^ 5 X 10" GeV ,^ , 

(4.5) 



2 Hz Ml 



the analogous quantity for the charged-lepton Yukawa interactions. For any 
value of Ml and Ki, there is a value Zf[ such that Fr > Wid for z > z^. 
If Ml < 2 X 10^2 GeV/i^i and Ki > 3.3, corresponding to Fr > Wwiz^a^) 
in the strong washout regime, then Zd — 0, meaning that flavor effects are 
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important during tfie entire thermal history. On the other hand, for a fixed 
value of Ml, one has 2;^ — > cx) for Ki 00, implying that flavor effects 
tend to disappear for sufficiently large values of Ki. Notice however that if 
Ml > 5 X 10^^ GeV, then zq > Zos — z-q for any value of Ki. This confirms 
that only for Mi > 5 x 10^^ GeV flavor effects can be neglected and the 
unflavored regime is recovered. 



4.1.2 Maximum flavor effects 

We now turn to the opposite extreme case when flavor effects are maximal, 
the "fully flavored regime." In other words, the charged-lepton Yukawa inter- 
actions are now taken to be so fast that the lepton flavor content produced 
in iV — ^ £ + $ on average fully collapses before the inverse reaction can take 
place, i.e., the I density matrix in flavor space is to be taken diagonal in 
the charged-lepton Yukawa basis. In this case each single-flavor asymmetry 
has to be calculated separately because generally the washout by inverse 
decays is different for each flavor. Moreover, the single-flavor CP asymme- 
tries now have an additional contribution compared to the total, as shown 
in Eq. (3.6) [119,120]. Finally, the inverse decay involving a lepton in the 
flavor a does not wash out as much asymmetry as the one produced by one 
RH neutrino decay. The reduction is quantified by the probability Pf^, aver- 
aged over leptons and antileptons, that the lepton li produced in the decay 
of Ni collapses into the flavor eigenstate I a- Focusing on the lightest RH 
neutrino iVi, the fully flavored Boltzmann equations are given by Eqs. (3.17) 
and (3.18), with i — 1. Since we are dealing with the two-flavor case, here 
a — T or e + II where the latter stands for a suitable superposition of the e 
and flavors, as explained after Eq. (3.18). 

As in the unflavored case, we next identify the condition for the fully 
flavored approximation to hold. The flnal asymmetry in the flavor a is domi- 
nantly produced at 2; ~ ZBa = ZB{Kia), where Kia = Pi^ Ki. Therefore, one 
must require that Fq, > Fi° holds already at z ~ ZBa, or else the washout 
reduction takes place too late. We stress that flavor effects modify the fl- 
nal asymmetry only if the flavor projection takes place before the washout 
by inverse decays freezes out. Otherwise the washout epoch is over and the 
unflavored behavior is recovered. It is easy to verify that if the projectors 
are set to unity and the equations are summed over flavors, the kinetic equa- 
tion (2.28) for Nb-l holding in the unflavored regime is recovered. Therefore, 
we require 

as an approximate condition for the fully flavored behavior. 
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Figure 4.2: Relevance of flavor effects in schematic regions of parameters 
Ki and Mi. The region above the horizontal dashed line corresponds to the 
condition (4.3). The vertical dot-dashed hne is the border between the weak 
and the strong washout regime. The region below the inclined dotted line 
corresponds to the condition (4.6) for ZBa = ^max- 



In Fig. 4.2 we summarize the different possible cases in the plane of pa- 
rameters Ki and Mi. For Mi > 5 x 10^^ GeV, above the dashed line, flavor 
effects arc not important independently of Ki. The condition (4.6), in the 
most restrictive case when ZBa = ^max and Wl^ ~ 0.3^1, is satisfied below 
the inclined dotted line. This case typically occurs in a one-flavor dominated 
scenario, as we explain below. The vertical dot-dashed hne is the border that 
separates the weak from the strong washout regime in the unflavored case. In 
the flavored case the condition Ki < 3.3 still implies a weak washout regime 
because flavor effects can only reduce the washout. However, the condition 
for the strong washout regime can be more restrictive than Ki > 3.3, as 
discussed in [123]. For Ki < 3.3, flavor effects modify the final asymmetry 
only marginally and more specifically only if the initial A^i-abundance van- 
ishes, as indicated in Fig. 4.2. On the other hand, for Ki > 3.3 and below 
the diagonal line, flavor modifications of the final asymmetry can be large, 
especially in the one-flavor dominated scenario. 

There is a region in parameter space where neither condition (4.2) 
nor (4.6) holds. This intermediate regime can become very large in the case 
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of a one-flavor dominated scenario, where several orders of magnitude en- 
hancement of the flnal asymmetry compared to the unflavored calculation 
are possible (see Fig. 3.6). In this case one of the two projectors is very small 
compared to the other, and so the washout is very asymmetric in the two 
flavors. On the other hand, if the two flavored CP asymmetries are compa- 
rable, then the flnal asymmetry is dominantly produced into one flavor and 
deviations from the unflavored regime can become very large. This scenario 
is realized, in particular, when the absolute neutrino mass scale increases, 
relaxing the traditional neutrino mass bound. 

However, the condition (4.6) strongly restricts the applicability of the 
one-flavor dominated scenario. Even though the reduction of the washout is 
driven by Ki^ <^ Ki, implying z^^ ^ ^b, the possibility for flavor effects to 
be relevant relies on the dominance of the charged- Icpt on Yukawa interaction 
rate compared to the RH neutrino inverse-decay rate, which however is still 
driven by Ki. Therefore, increasing Ki, one can enhance the asymmetry in 
the one-flavor dominated scenario compared to the unflavored case, if smaller 
and smaller values of the projector P^^ are possible. On the other hand, the 
inverse-decay rate increases so that the fully flavored behavior may no longer 
apply. In particular notice that the maximum enhancement of the asymmetry 
is obtained when Ki » — 1, when z^a — -^max and the condition (4.6) 
is maximally restrictive. 



4.2 Neutrino mass bound 

One possible consequence of flavor effects is to relax the traditional upper 
bound on the neutrino mass that is implied by successful Icptogenesis. In 
order to explore the impact of our modified criteria, we first recall the ori- 
gin of this bound in the unflavored case. Maximizing the flnal value of the 
asymmetry over all see-saw parameters except Mi and mi yields [80] 
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(4.7) 



where we have approximated k{Ki) ~ 0.5 i^f ' [cf. Eq. (2.35)], and wc 
have neglected the dependence of zb on Ki in the derivative. This constraint 
translates into mi < m™'^'^(Mi) shown by the curved solid line in the upper 
part of Fig. 4.3 where the unflavored behavior obtains. This curve sports an 
absolute maximum, mi < 0.12 eV, for Mi ~ 10^^ GeV. 
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Figure 4.3: Relevance of flavor effects similar to Fig. 4.2, now mapped to 
schematic regions of parameters mi and Mi . The region above the horizontal 
dashed line corresponds to the condition (4.3) for the applicabihty of the un- 
flavored regime. The region below the inclined thick dashed line corresponds 
to the condition (4.6) calculated for that value of ZBa that maximizes the fi- 
nal asymmetry in the one- flavor dominated scenario and for Ki = mi/m*. In 
this same case, the lower inchned dot-dashed hues includes also the effect of 
scatterings in the condition (4.6), while the upper inclined and the horizontal 
dot-dashed lines include the effect of oscillations. The area between the two 
inclined dot-dashed lines gives an estimation of the uncertainty on the con- 
dition for the fully flavored regime to hold. The area between the horizontal 
thin dot-dashed line and the horizontal thick dashed line gives an estima- 
tion of the uncertainty on the condition for the unflavored regime to hold. 
The two thick solid lines borders the region where successful leptogenesis is 
possible: on the left in the unflavored regime and above in the fully flavored 
regime. The thin solid line is a more restrictive border obtained for a specific 
choice of the see-saw orthogonal matrix [Q = Ri^) and the thin dashed line 
is the corresponding condition (4.6). In this case one has Ki > mi/m*. 
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The possibility that flavor effects could relax this bound is based on the 
observation that the flavored CP asymmetries, for rui > matm, are propor- 
tional to mi [cf. Eq. (3.30]. On the other hand, the total CP asymmetry 
is suppressed like m^^, contributing to the upper bound in the unflavored 
regime. However, if mi increases, then Ki has to increase as well, so that it 
is not guaranteed that the fully flavored treatment remains justified. 

Quantitatively, the value of Ki is bounded from below by [92] 

K,>^. (4.8) 
m* 

For Ki > mi/m* ^ 1, the final asymmetry is maximized when a one-flavor 
dominated scenario is realized. In this case the final asymmetry is approxi- 
mately N^B-L — ^lai^ia- Thc bouud on the single-flavor CP asymmetry was 
given in Eq. (3.30). It is then possible to find the value of P^^ that maxi- 
mizes the asymmetry as a function of Ki and the corresponding value of ZBa- 
Imposing rj^^ > Vb^^ imphes a lower bound on Mi as a function of Ki. 

This limit can be translated into a lower bound on Mi as a function of 
mi by replacing Ki with its minimum value mi/m^. In this way the washout 
is always minimized and the final efficiency factor and the final asymmetry 
are maximized. Notice that the single-flavor CP asymmetries, like the total, 
vanish for Ki — mi/my^. Therefore, this lower bound cannot be saturated. 
Notice moreover that for mi < m^ the one-flavor dominated scenario does 
not necessarily hold because it is possible that Ki < 1. Actually, for Ki — * 0, 
flavor effects disappear and one recovers the usual asymptotic value of the 
lower bound obtained in the unflavored case for thermal initial A^i-abundance, 
Ml > 5 X lO^GeV [cf. Eq. (2.53)]. For intermediate values of Ki, one can 
use a simple interpolation. The flnal result is shown in the bottom part of 
Fig. 4.3 as a thick solid line. 

In Fig. 4.3 we also show the condition (4.6) calculated for the same value 
of ZBa that maximizes the flnal asymmetry , but replacing Ki with its min- 
imum value mi/m^t (thick dashed line). Since Wf^ increases with Ki, this 
produces a necessary, but not sufficient, condition in the mi-Mi plane for the 
fully flavored behavior. This condition matches the validity of the unflavored 
regime at mi ~ 3 x 10~^ eV and the lower bound on Mi at mi ~ 2 eV. This 
means that for mi > 2 eV the fully flavored behavior does not obtain. No- 
tice also that this upper limit is quite conservative because the lower bound 
on Ml has been obtained neglecting that, for Ki = mi/m*, the flavored 
CP asymmetry vanishes and thus the bound cannot be saturated. Moreover, 
we have assumed that P°q, can always assume the value that maximizes the 
asymmetry. 

In Fig. 4.3 we also show (thin solid line) the lower bound Mi (mi) in 
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the specific scenario considered in Section 3.5, namely Q = R13 and uj^i is 
taken purely imaginary [cf. Eq. (2.21)]. In this case the value of ZBa is not 
necessarily the same as the one that maximizes the asymmetry in the one- 
flavor dominated scenario, and Ki > mi/m*. Therefore, a specific calculation 
is necessary in order to work out correctly the condition (4.6). The result is 
shown in Fig. 4.3 with a thin dashed line. 

In this case the upper limit on nii for the applicability of the fully flavored 
regime is much smaller, mi ~ 0.1 cV. Allowing for a non- vanishing real part 
of Wgi, slightly larger values are possible. It should be however kept in mind 
that these values are indicative since they rely on a condition for the fully 
flavored regime that comes from a simple rate comparison. 

4.3 Limitations of a simple rate comparison 

We have exploited a somewhat qualitative rate comparison for the determi- 
nation of the region where the fully flavored regime obtains. While we believe 
that our approach nicely illustrates the modifications that derive from our 
more restrictive criterion for the significance of flavor effects, there are also 
important shortcomings. First, we have simply compared the inverse-decay 
rate with the charged-lepton Yukawa interaction rate, ignoring flavor os- 
cillations caused by the flavor-dependent lepton dispersion relation in the 
medium. If the oscillations are much faster than the inverse-decay rate, they 
also contribute effectively, together with inelastic scatterings, to project the 
lepton state on the flavor basis. Therefore, including oscillations will tend 
to enlarge the region where the fully flavored behavior obtains (the inclined 
upper dot-dashed line in Fig. 4.3) and to reduce the one where the unfla- 
vored behavior obtains (the horizontal upper dot-dashed line in Fig. 4.3). 
In our case, the oscillation frequency is comparable to F^, and so the two 
estimations are not too far off. 

Moreover, we have also neglected AL = 1 scatterings. They also con- 
tribute, like inverse decays, both to generate the asymmetry and to the 
washout and hence, together with inverse decays, contribute to preserving 
the flavor direction of the leptons. At the relevant z ~ ZBa ^ 2, the AL = 1 
scattering rate is actually larger than the inverse-decay rate and thus tends 
to reduce the region where the fully flavored behavior obtains (the lower in- 
clined dot-dashed line in Fig. 4.3). Therefore, the effects of oscillations and 
of AL = 1 scatterings may partially cancel each other. In Fig. 4.3 the region 
between the two inclined dot-dashed lines gives then an indication of the 
theoretical uncertainty on the determination of the region where the fully 
flavored regime holds. It can be seen that current calculations cannot estab- 
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lish whether the upper bound holding in the unflavored regime is nuUified, 
just simply relaxed or still holding, when flavor effects are included. 

Only a full quantum kinetic treatment can give a final verdict on the 
effectiveness of flavor effects in leptogenesis and its impact on the neutrino 
mass limit. While we have verified that our rate criteria are borne out by the 
quantum kinetic equations stated in Ref. [120] (see [134]), these equations 
are not necessarily complete in that the term describing the generation of the 
asymmetry has been added by hand. Moreover, the "damping rate" caused 
by the flavor-sensitive Yukawa interactions ultimately derives from a collision 
term in the kinetic equation [135]. Extending the pioneering treatment of 
Ref. [120] to allow for a complete understanding of flavor effects remains a 
challenging task. In the next section we show how such a quantum kinetic 
equation could be derived. 

4.4 Density matrix equation 

Here we aim to derive a density matrix equation for leptogenesis follow- 
ing [135], where the general formalism was introduced. We consider exclu- 
sively the decay of the lightest RH neutrino, A^i, with a mass between 10^ GeV 
and 10^^ GeV where only the r- Yukawa interactions may be faster than the 
Hubble rate. We are therefore either in the unflavored or in the two-flavor 
regime, the two flavors being r and a combination of e and /j, which we call 
P in the following. 

We follow the definition of the matrices of densities for leptons and an- 
tileptons given in [135]: 



lation (creation) operator for a lepton and b (b^) the annihilation (creation) 
operator for an antilepton. All matrix components are implicitly given in 
a portion of comoving volume that would contain one heavy neutrino in 
ultra-relativistic thermal equilibrium. Note that we introduce here matrices 
of number densities, not occupation numbers as in [135]. We shall explain 
later on how the integration over all modes can be performed here. As usual 
in leptogenesis, one assumes kinetic equilibrium, an assumption which should 
be good at the 10% level in the strong washout regime [136]. Following this 
assumption, all interaction rates are thermally averaged. 




(4.10) 



(4.9) 




denoting the annihi- 
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Let us first study the two simplest - nevertheless very important - pro- 
cesses, namely decays and inverse decays: A^i <-> £i + and A^i <-> + 
The state was defined in Eq. (2.1). 

The idea is to write the equation for the generation of lepton number 
in the basis {|^i), |^±)} = {^i}, and then rotate it into the fiavor basis 
{\^t), I-^/?)} = {^o}- The unitary matrix for the change of basis is given by 



u{{h} - {4}) 



11 



(4.11) 



The same can be done for antileptons, and one gets the following matrix for 
the change of basis: 



u'm - {4}) 



(4.12) 



where it should be remembered that I\ is not the CP-conjugate of i\ at the 
one-loop level. 

Decays and inverse decays are just source and sink terms, so one can 
proceed by analogy with [135] with the conventions of [80] on leptogenesis 
quantities. In the basis {^i}, the equation tracking the density matrix for 
leptons is given by 



1 




-2(l-£i)W^iD 



1 




Pl^ Pi 



where z — Mi/T, and for antileptons: 



1 

(4.13) 



PI 



1 \ ^ [10 
)P~i + Pi[o 

(4.14) 

where both equations are written to first order in the CP asymmetry param- 
eter El, and we neglected Pauli blocking effects. 

Let us now rotate these equations to the flavor basis applying the 

transformations introduced above, Eqs. (4.11) and (4.12). It will prove useful 
to define the following matrices which will appear in the new equations: 



and 



U 



U' 



1 




1 




Wr)Wl) Pip 



{A\lr){l',\h) PW 



r, 



= v, 



(4.15) 



(4.16) 
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where the projectors are defined as Pia — |(-^aKi)P and Pia — 

[cf. Eqs. (3.3) and (3.4)]. At tree level, the two matrices just defined are 

equal and given by 

In the flavor basis, Eqs. (4.13) and (4.14) become 
d 1 X 

—pe =2(1 + ^i)DNnJ> - -(1 - £1) WiD [Ppi + piV] , (4.18) 

and 

^p,- = ^(1 - e,)DNj,^V - ^(1 + e,)WiT, [Vpi + piV\ . (4.19) 

Subtracting the flrst equation with the second, and deflning p(^—pi = Pi_i^ 
one obtains 

^^p,_l = s,D {N^, + N^) V' + {N^, - N'^J - ^W,n{V', p,^i}, 

^ _ (4.20) 

to flrst order both in £1 and V — V. Note that we used the fact that 
P£ + Pg^2Npi, where we neglected contributions of 0{ei), because the 
whole term is already at first order in Ei. As usual in leptogenesis, one no- 
tices that there is an asymmetry production even in thermal equilibrium [82] . 
This is due to the missing contribution from the on-shell AL = 2 processes 
Properly accounting for them yields [120, 127] 

= e,D {Nr,, - K^^) (v' + - VID{P^ (4.21) 

which looks intuitively correct, with both the reduction of the washout by 
the projector (a matrix now!) and the additional contribution to the source 
term, of the type AP (also a matrix). If the off-diagonal elements are quickly 
damped (see below), one recovers exactly the fully fiavored equations on the 
diagonal entries [cf. Eq. (3.18) with i = 1]. 

Actually, the matrix of CP asymmetries, which can be defined as 

was proposed in [120] for the three-fiavor regime to be given by 
^"''^ ^ 327r(/itl)^i XI {Kiih^h)ijh0j - hf3i{h^h)jihlj) , 



(4.22) 
(4.23) 
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Figure 4.4: Flavor diagonal contribution to the refractive index for the lep- 
tons. 

in the HL, Mi <^ <^ M3, where = e,fi,T. It is straightforward to 
check that the diagonal elements correspond to the single-flavor CP asym- 
metries shown in Eq. (3.25). 

Let us now examine the effect of the term fa^ia^Ra^, which we said 
in Section 3.1 is responsible for flavor effects to manifest themselves. This 
interaction is flavor diagonal, and we expect it to affect Eq. (4.21) in two 
ways. First, it induces a contribution to the refractive index for the lepton 
doublet ia, and, second, it makes the lepton doublet interact with the lepton 
singlet Ca- 

Starting with the index of refraction effect, in perfect analogy with [135], 
one has a contribution in the form of a commutator both for leptons and 
antileptons: 

^p, = [RHSofEq (4.18)]-z[A:,p,] (4.24) 
Ap^- = [RHSofEq (4.19)]+z[A:,p,-], (4.25) 

where 

Hz[ J - HzQA [00) ^ ' 

The effective mass Ma aquired by the lepton in the medium due to the 
Yukawa coupling /„ was calculated in [137] from the real part of the diagram 
shown in Fig. 4.4 and used in the second equality. 

It is important to be aware that the possibility to use a density matrix of 
number densities, i.e. to have been able to integrate over the 3-momentum 
of the leptons, is non-trivial and comes from [138], where it was noticed 
that when gauge interactions are very fast, all modes oscillate at the same 
frequency, given by the thermal average (M^/2p) in this case. 
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Similar to the effect of a vacuum mass, leptons and antileptons oscillate 
in opposite directions, so that subtracting Eq. (4.25) to Eq. (4.24) leads to 

^Pe-i= [RHS of Eq. (4.21)]-z[A:,p,+p,-] = [RHS of Eq. {4.21)]-t[Al, p_], 

^ (4.27) 
where, in the second equality, we defined p_ = pi + pi—2Npi, which has all 
elements of 0{ei) since the diagonal part of p£ + pg has a part proportional 
to 2Np which drops out of the commutator. In order to have a closed set of 
equations, we must give the equation of evolution of p_, 

^P_ = z[A:,P,_,-]. (4.28) 

As for the new interactions between lepton doublets and singlets, they 
yield again two terms, one source and one sink term, which describe how 
the lepton asymmetry is shared between the doublets and the singlets. As 
we already mentioned, only the r- Yukawa coupling is relevant for the range 
of RH neutrino masses we consider, and thus the relevant interaction rate is 
given by F^, which we defined in Eq. (4.5). Assuming a negligible asymmetry 
in the Higgs field, the resulting term can be written in the form 

+ - ( J °)-iF.{(J °).ft_,-}.(4.29) 

The asymmetry in the lepton singlets is found by solving the following 
Boltzmann equation: 

^(iVe^ - iVe-J = [{Ne^ - iV,-J - 2(iVe, - iVgJ] , (4.30) 

where Ng^ — N^^ is the upper left component of the density matrix 
Clearly, this equation will tend to equalize Ni^ — Ni^ and 2{Ne^ — Ng^). It is 
interesting to notice that when the latter two quantities are equal, Eq. (4.29) 
can be written as 



1 



1 




1 




:PL 



(4.31) 



where the total lepton number Nl is defined as {Ni — Nj) + (iVg — iVg). 
In particular, this means that the upper left component of pl is given by 
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2{N£ — Ni) /3. The coefficient 2/3 translates the fact that some of the lepton 
asymmetry, being stored in the right-handed fields, escapes the washout by 
inverse decays. This is nothing else than an effect caused by one spectator 
process, as discussed in Section 3.2, with the difference that we consider here 
a simplified situation where the asymmetry stored in the Higgs is assumed 
to be zero. 

Eq. (4.31) is very illustrative in that it includes all the well-known con- 
tributions discussed in [135] in a completely different context, namely the 
production (source term), the washout (sink term), the double commutator, 
which damps the off-diagonal elements, and the commutator, which drives 
the oscillations in flavor space. 

In the transition region T ~ 10^^ GeV, one cannot assume N^^ — N^^ — 
2{Ne^ — Ng^), and one has to solve simultaneously four coupled equations, 
Eqs. (3.7), (4.28), (4.29) and (4.30). They constitute the system of equations 
to solve in order to have a correct description of leptogenesis at the transition 
between the fully flavored and the unflavored regime, with an approximate 
treatment of spectator processes. 

The system of equations we propose differ from what can be found in the 
present literature [120, 134] in that 

1. we explicitly include the effect of oscillations in flavor space, which 
adds one equation, since leptons and antileptons oscillate in different 
directions in flavor space; 

2. we keep track of the asymmetry in e^, which adds another equation, as 
should be done in the transition regime around T ~ 10^^ GeV. 

The second point is actually the addition of one spectator process to the 
picture. In principle, all spectator processes should be taken into account, and 
a large number of Boltzmann equations should be solved simultaneously. In 
particular, sphalerons should be included and the final asymmetry obtained 
would then be in Aq, = S/3— L^, as required. We did not include these effects 
here for the sake of simplicity and clarity, as well as because the difference is 
expected to be of order 20-30%. 

The conditions for the validity of the different pictures of leptogenesis 
which we discussed in the previous sections can be put in perspective by 
looking at Eq. (4.29). Roughly speaking, if Wid » -^t, A^, one expects the 
unflavored regime to be recovered, whereas in the opposite situation, the 
fully flavored regime should hold. To follow up on the comments made in the 
last section about the inclusion of AL = 1 scatterings, it should be noted 
that, even though we did not include them for simplicity in our analysis, 
their addition is straightforward in that one has to replace Di ^ Di + Si 
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and Wl^ Wl^ + W^^^^ — j{z)WiB, with the corresponding analytic 
expressions given in Eqs. (2.6) and (2.11), respectively. 

Finally, let us note that there is another quantum limit which can be 
considered for the Boltzmann equations relevant for leptogenesis. It relies 
on the closed time path formalism for non-equilibrium quantum field theory. 
Within such an approach, the so-called "memory effects" are for instance 
accounted for. In the case of leptogenesis, these effects translate into a time- 
dependent CP asymmetry. This quantum limit was studied in [139, 140], and 
the main conclusion reached there is that memory effects have some impact 
in the case of resonant leptogenesis in the weak washout regime. In all other 
cases, they can be safely neglected. 
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Chapter 5 



Going beyond vanilla 
leptogenesis 



We have always assumed up to now that the final B — L asymmetry was 
predominantly produced by the hghtest RH neutrino, A^i. A A^i-dominated 
scenario typically follows from the assumption of hierarchical heavy neu- 
trino masses, Mi <tC -C M3. Actually, it can even be enforced by setting 
R23 = 1, since a large contribution from N2 is then unprobable, as we ex- 
plained in Section 2.1. This is "vanilla" leptogenesis, where successful lepto- 
genesis requires high-scale values of the lightest RH neutrino mass. Mi, and 
of the reheat temperature, T^eh- Moreover, the upper bound on the absolute 
neutrino mass scale implied by successful leptogenesis was discussed in Sec- 
tion 2.2 within an unflavored treatment and in Section 4.2 including flavor 
effects. 

In this chapter we would like to go beyond this picture. Actually, flavor 
effects themselves allow that to some extent. For example, they might modify 
the upper bound on the absolute neutrino mass scale, and they open new 
ways to have N2- and A^3-dominated scenarios. First, we want to relax the 
assumption of a hierarchical heavy neutrino mass spectrum, and concentrate 
on quasi-degenerate spectra. Mi ~ M2 ~ M3. Then, we discuss the impli- 
cations of having R23 7^ 1, opening the way to the iV2-dominated scenario. 
Actually, we also analyze this possibility within a flavored perspective and 
notice that the window for A^2- and A'a-leptogenesis opens up compared to 
the unflavored case. Finally, we discuss the possible effects coming from one 
element of the Q, matrix, namely 0,22- 
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5.1 Degenerate limit for the heavy neutrinos 

On theoretical grounds it is not difficult to motivate models with a quasi- 
degenerate heavy neutrino mass spectrum. For example, a slightly broken 
U{1)l lepton flavor symmetry is enough to achieve this goal [141-143]. Note 
that this possibility has other implications as well, as we shall see below. 
Such a spectrum can also be motivated in the context of "radiative leptoge- 
nesis" [144,145]. 

In leptogenesis a quasi-degenerate mass spectrum for the heavy neutrinos 
leads to important qualitative and quantitative differences with respect to 
the vanilla picture described in Chapters 2 and 3, because all heavy neu- 
trino contributions must be taken into account, both at the production and 
washout levels. 

It is straightforward to generalize, including flavor effects, a result ob- 
tained in [94] for the efficiency factors within an unffavored treatment in the 
degenerate hmit, (M3 — Mi) /Mi < 0.1, and for Ki S> 1, for all i. Indeed, one 
can now approximate dNi/dz' ~ dN^'^/dz' in Eq. (3.20), obtaining that 

Kl^C^K{Kia, + K2a + K3a). (5.1) 

The function k,{x) was deflned in Eq. (2.31), and it approximates k.\^ in the 
hierarchical limit when x = Ki^. In the degenerate limit (DL) one has then 
only to replace Kia with the sum Kia + i^2a + K^a- The number of efficiency 
factors to be calculated reduces from 6 to 2 in the two-flavor case and from 
9 to 3 in the three-flavor case, i.e. one for each flavor, like in the hierarchical 
limit. If, instead of a full degeneracy, one has only a partial degeneracy. 
Ml ~ M2 -C M3, then 

4a < «L - 4a - '^(^la + K^a)- (5.2) 

It is interesting to notice that, as a consequence of the orthogonality of 
O, one has 

Ki^ + i^2a + i^3a = V — \U^k? ■ (5.3) 

This means that the sum over the decay parameters will typically be in the 
strong washout range, Ki^ + -f^2a + -f^sa ^ 3. This is a nice property of 

leptogenesis with degenerate heavy neutrinos which still holds when flavor 
effects are included. Consequently, we do not need to introduce analytic ex- 
pressions for a vanishing initial number of heavy neutrinos, since it will give 
the same result as K,{Kia + i^2a + -^sa), which was derived for a thermal 
initial A^i-abundance. 
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Concerning the flavored CP asymmetries, one can rewrite the general 
expression (3.24) in the DL, 

'''' ^ ImW ^ ^^"^ KhaAh^h),,] + Im [hl,ho,,{h^h),,] } SJ^ , (5.4) 
where we used the fact that ^{xj/xi) ~ l/(35jj) and we defined 

As it is clear from Eq. (5.4), there can be an enhancement of the CP asymme- 
try proportional to 5~l^ [87] for quasi-degenerate heavy neutrino masses. This 
enhancement originates from the one-loop self energy contribution. Note how- 
ever that an enhancement of the CP asymmetry might require that all three 
heavy neutrinos are quasi-degenerate. As a matter of fact, when Q = i?i3, 
there is no enhancement when M2 — > Mi, but only when M3 — > Mi [94]. 
We can now write in general for the final B — L asymmetry 

where the sums over i and k run from 1 to 3 in the case of full DL and from 
1 to 2 in the case of partial DL. When computing the final asymmetry and 
the lower bound on Mi in the DL, there will be a competition between two 
effects, namely the increase of the washout because of the sum of the decay 
parameters in the efficiency factor, and the enhancement of the CP asym- 
metry. The second effect will become dominant once the first one saturates, 
which occurs for 5ji ~ 0.01 [94]. From this point going to higher degeneracies 
will lower the bounds proportionally to Therefore, as we pointed out in 
Section 2.3, it is possible to relax the lower bound on the reheat temperature 
in such a way that the gravitino problem is avoided. 

One can go even further and look for the most extreme relaxation one can 
achieve. This occurs when the resonance in the CP asymmetry is reached, 
for Mj — Mi ~ Vj/2 [117, 118]. We shall use a slightly modified condition. 



ci£(M,)/3, (5.7) 



where we introduced the uncertainty parameter d = 1 10 because of the 
claim in [85] that it is not possible to be exactly on resonance without break- 
ing perturbation theory, contrary to what was stated in [76, 117]. We do not 
aim here at a resolution of this discrepancy and therefore introduced the 
parameter d. 
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It is straightforward to check that the resonance condition (5.7) imphes 
that Si — 1/d in the unflavored case with maximal phase. In other words, 
the resonance condition simply cancels out the dependence on the heavy 
neutrino mass scale in the CP asymmetry and therefore also in the final 
baryon asymmetry. Hence, in resonant leptogenesis there is essentially no 
lower bound on Mi and Treh and one can go down to the TeV scale [76, 117]. 
The only requirement is actually to let the heavy neutrinos decay before the 
freeze-out of the sphalerons (at around 100 GeV), in order to produce not 
only lepton number but also baryon number. In [76] the authors followed an 
unflavored treatment where the third heavy neutrino is decoupled and the 
other two are quasi-degenerate. In such a scenario one cannot avoid the fact 
that the Yukawa coupling must be quite small in order to have a successful 
leptogenesis. 

Accounting for flavor effects in the case of resonant leptogenesis, it was 
realized in [132, 146] that an even more dramatic situation could be envisaged. 
If three quasi-degenerate RH neutrinos are considered, it is possible to have 
some of the Yukawa couplings of order one, leading to phenomenological 
implications, such as a rate for the lepton flavor violating decay — > 67 
within reach of future experiments, as well as successful leptogenesis at the 
TeV scale. The crucial point there was to have one flavor, the r, that is 
very weakly washed out thanks to a very small projector Pi^-, and one heavy 
neutrino, A^3, that is weakly coupled. However, the latter condition does not 
necessarily imply that the CP asymmetry £3 is suppressed because it receives 
other contributions when M3 ~ Mi^2, as we will explain in Section 5.3. It 
must be said, however, that the extreme situation needed seems to occur in 
a region of the parameter space where the condition of validity (4.6) which 
was derived in the previous chapter is not satisfled. 

Actually, strictly speaking, the possibility of explaining the smallness of 
neutrino masses with large Yukawa couplings and relatively light RH neutrino 
masses (TeV scale) does not involve a conventional "see-saw mechanism", 
even though the same matrix for the light neutrinos applies, Eq. (B.5). The 
point is that the neutrino masses are not in this case small because of the 
suppression due to the heavy neutrino mass, but because of a very large 
cancellation due to the O, matrix. The elements of the Q matrix can indeed 
be arbitrarily large, as long as they satisfy Q^Q = QQ^ = 1. This can 
be theoretically motivated by a slightly broken U{1)l lepton flavor symme- 
try [141-143]. It is interesting that this type of model may lead to possible 
signatures of heavy neutrinos at future coUiders [143, 147, 148] . 



5.2 A^2-dominated scenario 



83 



5.2 Ar2-dominated scenario 

In this section we consider again a hierarchical heavy neutrino mass spectrum, 
Ml <^ M2 <^ M3. In the unflavored regime, this assumption typically implies 
a A^i -dominated scenario, where the final asymmetry is dominated by the 
contribution from the lightest RH neutrino decays, Eq. (3.29). Indeed, in 
general, in the HL one has two effects. The first effect is that the asymmetry 
production from the two heavier RH neutrinos, N2 and A^3, is typically later 
on washed out by the A^i inverse processes, implying k^, k.^ -C k.{. The second 
effect is a consequence of the fact that the total CP asymmetries vanish in 
the limit when all particles running in the loops become massless, and this 
yields typically \es\ ^ 1^2 1 ^ 

However, for a particular choice of the see-saw parameters, Q = R23 
[cf. Eq. (2.22)] and mi -C [cf. Eq. (1.27)], the contribution to the fi- 
nal asymmetry from the next-to-lightest RH neutrino N2 is not only non- 
negligible but even dominant, giving rise to a iV2-dominated scenario [93]. 
Indeed, for Q = R23, different things happen simultaneously. First, N2, even 
though decoupled from A^^i, is still coupled to N3 and in the HL the total CP 
asymmetry £2 not only does not vanish, since it receives an unsuppressed con- 
tribution from graphs where A^3 runs in the loops, but can even be maximal 
[cf. Eq. (2.25)]. On the other hand, one has now £1 = 0, since A^i is essen- 
tially decoupled from the other two RH neutrinos. At the same time, one has 
Ki = mi/m^, <^ 1, so that the washout from A^i inverse processes is negligi- 
ble. The final result is that 1^2 '^2! 3> 1^17^2 '^^^2 1' ^^d the final asymmetry is 
dominantly produced from A^2-decays. 

A nice feature of this model is that the lower bound on Mi does not 
hold anymore, being replaced by a lower bound on M2 that, however, still 
implies a lower bound on Tj-ch [93]. If one switches on some small -R12 and R13 
complex rotations, then the lower bounds on M2 and Treh become necessarily 
more stringent. Therefore, there is a border beyond which this scenario is 
not viable and one is forced to go back to the usual A^i -dominated scenario 
for successful leptogenesis. 

Within an unflavored treatment and in the HL, the condition 'u;32 ~ in 
the r2- matrix parametrization [cf. Eq. (2.19)], implying i?23 — 1, is sufficient 
to have a negligible asymmetry production from the two heavier RH neutrinos 
and to guarantee that the A^i-dominated scenario holds. This condition is 
even not necessary for mi m^, since in this case, due to the fact that 
m-i > mi, one has necessarily Ki ^ 1 and the washout from A^i inverse 
processes is strong enough to suppress a possible contribution to the final 
asymmetry from A^2-decays. 

When flavor effects are taken into account, the domain of applicability 
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of the A^i-dominated scenario reduces because the importance of N2 and 
increases. There are two aspects to be considered. 

The first aspect is that the washout from A^i inverse processes becomes 
less efficient. Indeed, the projectors Pi„ can considerably reduce the washout 
of the asymmetry produced in the flavor a from iV2-decays [149]. This turns 
the condition rrii ^ into a looser condition toi ^ rrii^/Pia- Another 
effect is that A^i inverse processes can be fast enough to quickly destroy the 
coherence of £2- Then a statistical mixture of £\ and of the state orthogonal 
to ii builds up, and hence part of the asymmetry produced in iV2-decays 
is protected from the washout from Ni inverse processes [6, 127, 150]. These 
effects may occur in the unflavored regime (Mi > 10^^ GeV) or in the two- 
flavor regime (10^ GeV < Mi < 10^^ GeV), but disappear in the three-flavor 
regime, where the full flavor basis is resolved, since no direction in flavor space 
is protected from the A^i-washout. Recently, it has been also pointed out that 
the off-diagonal elements in the matrix introduced in Eq. (3.15) and which 
encodes the effects of all spectator processes can lead to a reduction of the 
washout from A^i inverse processes as well [151]. One can therefore conclude 
that the assumption K2a <S i^ia is not valid in general, even when Mi <S M2. 

The second aspect concerns the flavored CP asymmetries. Even though 
the results obtained at the end of Section 2.1 in the particular cases Vl — 
i?23 and VL = Ris are still valid in the flavored case, i.e. eia = for fl — 
i?23, and £2a = for Q = it!i3, there are other effects to keep in mind. As 
we shortly mentioned in Section 3.3, the flavored CP asymmetries £2a are 
not necessarily suppressed by factors M1/M2 compared to sia, as it is the 
case in the unflavored regime for £2 compared to Si. This observation [123] 
can also potentially contribute to enlarge the domain of applicability of the 
A^2-dominated scenario when flavor effects are taken into account. Another 
related observation is that the Ssa's, contrarily to £3, are not suppressed in the 
HL. This could open the door even to a Ai'a-dominated scenario, although this 
is possible only for M3 < 10^^ GeV, when flavor effects affect the generation 
of asymmetry by N^. 

Therefore, when flavor effects are taken into account, the conditions of 
applicability of the A^i-dominated scenario become potentially more restric- 
tive than in the unflavored case. There is a clear choice of the parameters, 
for Q = 7?i3 and M3 > 10^^ GeV, where the iVi-dominated scenario holds. 
Indeed, in this case, one has that e2a = and £3 is suppressed as M1/M3. 
This can be considered somehow opposite to the case fl = R23, where the 
A^2-dominated scenario holds [93]. 

In general, one can say that the asymmetry produced from the two heavier 
RH neutrinos is non-negligible if two conditions are satisfied: 
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(i) The asymmetry generated from A^2,3-decays at T ~ M2,3 is non- 
negligible compared to the asymmetry generated at T ~ Mi from A^i- 
decays. This depends on an evaluation of the CP asymmetries £23 

of the washout due to the same N2^3 inverse processes. 

(ii) The asymmetry produced from A^2,3-decays is not afterwards washed 
out by A^i inverse processes. Notice that this second condition is sub- 
ordinate to the first condition. 

5.3 Effects of |ri22| 

The general formula for the total CP asymmetry Si was given in Eq. (2.16) 
with i = 1. It can be rewritten in the illustrative form Eq. (2.40), where 
the second term was said to be usually negligible for hierarchies larger than 
3Mi < M2. In fact, there is a situation where the term [^{x^) — ^(0:2)] A^i in 
Eq. (2.40) becomes dominant over ^{x2)e^^ and one has to require a stronger 
hierarchy than 3Mi < M2 to recover the results presented in Section 2.2 [98]. 
From another perspective one can say that the second term in Eq. (2.40) 
offers a way to evade the bound on the CP asymmetry Eq. (2.45) for mild 
hierarchies. 

The term proportional to Asi is maximized when ^ 1, implying 
^(xa) ~ 1. Moreover, if for definiteness one imposes Q21 — and Re[Q|i] = 0, 
so that sin^L = 1 [cf. Eq.(2.47)], then 

= = C(X2) + [1 - ^{x2)]{Re[nl2] + mjm^tm M^l^]). (5.8) 

e[Mi) 

One can immediately see that when Re[022] = 1 ^^[^22] — 0) corre- 
sponding to r2 = i?i3, one recovers ^^-^ = 1, independently of the value of X2- 
Notice also that when M2 = M3, the Aei term vanishes. However, one can 
now perceive another possibility: if 1^222] S> 1, then the CP asymmetry can 
be enhanced, i.e. > 1, even when X2 ^ 1- 

This possibility is interesting because having large |f222| values does not 
imply that the washout is large. As a matter of fact, the decay parame- 
ter Ki does not depend on fl22- So there can be an enhancement of the 
CP asymmetry without any enhancement of the washout. This observation 
was made in [98] and supplemented with theoretical motivations for models 
with large \0,22\ in [152]. These models make possible to evade the bound 
on the CP asymmetry Eq. (2.45) without reducing the efficiency factor, thus 
relaxing the lower bounds Eqs. (2.54) and (2.55). It was found that values of 
Ml ~ 10^ GeV for mild hierarchies M2/M1 ~ 10 were possible, hence avoid- 
ing the gravitino problem. Note moreover that models with large IQ22I can 
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yield observable signals at experiments searching for lepton flavor violating 
decays n ej, t ^ fx'j and t —>■ ej [5, 153, 154]. 

There is another important consequence of the term A^i, which concerns 
the upper bound on the absolute neutrino mass scale, Eq. (2.57). First, when 
the absolute neutrino mass scale increases, it is well known that the upper 
bound on the CP asymmetry decreases [cf. Eq. (2.45)]. Thus, a dominance 
of the Asi term over the first term in Eq. (2.40) can be more easily ob- 
tained when mi is larger, even for moderate values 1^22! ~ 1- The term 
A£i is indeed not suppressed at large mi. Second, it can be argued that 
when nil ~ 0.1 eV, the light neutrinos are quasi-degenerate, so that quasi- 
degenerate heavy neutrinos with a mild degeneracy are more natural. Under 
the assumption of 'natural' mild degeneracy, the effect of the A^i term has to 
be included and leads to a relaxed upper bound rrii < 0.6 eV [98]. Note that 
the washout from the heavier RH neutrino A^2 must be taken into account 
for degeneracies S21 < 0.1, as emphasized in Section 5.1. 

The bottom line is that there are situations where the Aei term in 
Eq. (2.40) is crucial and hence should not be neglected. The first situation 
is when a large 11^221 allows to evade the upper bound on the CP asymme- 
try even when mi — > and for moderate hierarchies M2/M1 ~ 10, without 
reducing the efficiency factor. The lower bounds on Mi and Tj-eh can then be 
relaxed by a few orders of magnitude. The second situation is when a quasi- 
degenerate mass spectrum for the heavy neutrinos is considered, 5ji < 1, and 
\^22\ ~ 1- Then the contribution from the Asi term can be large, especially 
for quasi-degenerate light neutrinos. 

We think it is judicious to recall that Eq. (2.40) is only a convenient 
rewriting of the general expression (2.16). All the effects discussed in this sec- 
tion are of course present in the general expression, but they can be overseen 
when the basic assumptions behind some simplified expressions like in 
Eq. (2.40) are forgotten, such as a sufficiently large hierachy Mi <C M2 <C M3. 



Chapter 6 

Leptogenesis from low-energy 
CP-violating phases 

The possibility of relating the CP violation required for successful leptoge- 
nesis with the one that could be seen in future neutrino experiment is very 
attractive. However, in the context of unflavored leptogenesis (see Chapter 2), 
we saw that the PMNS matrix, which includes the observable CP-violating 
phases [cf. Eq. (A.5)], cancels out in the general case. Yet, even in this unfa- 
vorable situation, there have been attempts to relate low-energy CP- violating 
phases with the baryon asymmetry of the Universe produced through lep- 
togenesis [79, 155-160]. It must be however said that a link could only be 
made for specific textures of the Dirac mass matrix. 

As discussed in Chapter 3, flavor effects have modified the conventional 
picture of leptogenesis in a number of ways. In particular, the fact that the 
final asymmetry in the fully flavored regime depends explicitly on the PMNS 
mixing matrix is interesting. Thanks to the new source of CP violation im- 
plied by flavor effects, namely the AP contribution, Eq. (3.6), one can imag- 
ine a situation where the only source of CP violation comes from the PMNS 
matrix or, equivalently, where the contributions from the unobservable phases 
in the Vt matrix are negligible. This leads to the exciting possibility of ex- 
plaining the baryon asymmetry of the Universe thanks to CP- violating phases 
that are accessible in neutrino experiments. This scenario has attracted some 
attention recently [123,161-166]. 

In this chapter we first introduce the concepts of low-energy CP violation 
due to the Dirac phase and to the Majorana phases. Then, we study the 
possibility that the Dirac phase, which offers the best prospects for a possible 
measurement in the next years, provides the unique source of CP violation 
required for leptogenesis. This can be considered as the most conservative 
case, since we know that the Dirac phase comes always together with the 
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small angle 9is [cf. Eq. (1.8)]. Afterwards, we shortly comment on the role 
played by the Majorana phases as sole sources of CP violation in leptogenesis. 
Finally, we discuss the theoretical relevance of models where fl is real. 



6.1 CP violation in neutrino physics 

6.1.1 Neutrino oscillations and the Dirac phase 

Searching for CP-violating effects in neutrino oscillations is the only prac- 
tical way to get information about Dirac CP violation in the lepton sector, 
associated with the phase 6 in the PMNS mixing matrix U [cf. Eq. (A. 5)]. A 
measure of CP and T violation is provided by the asymmetries [167-169] 

where a ^ a' = ii, r. For three-neutrino oscillations in vacuum, which 
respect the CPT symmetry, one has [170] 

_ A{ti,r) _ 4(e,r) _ j pvac /i _ /i 
■^T ~ ~ ~ 'J CP ^osc ■> -^GP ~ ' 

J CP = Im {UeiU^2U*2U*i) = ^ sin 26112 sin 26123 cos^ 6^13 sin 6113 sin S, 

where x is the distance travelled by the neutrinos, and E is their common 
energy [cf. Eq. (1.11)]. Thus, the magnitude of the CP-violating effects in 
neutrino oscillations is controlled by the rephasing invariant associated with 
the Dirac phase S, the so-called Jarlskog invariant Jqp [171]. The existence 
of Dirac CP violation in the lepton sector would be established if, e.g., 
some of the vacuum oscillation asymmetries A^^p^ are proven experimen- 
tally to be non-zero. This would imply that Jcp 7^ and, consequently, that 
sin6'i3sin5 7^ 0. Without doubt, the search for CP-violating effects due to 
the Dirac phase in U represents one of the major goals of the future exper- 
imental studies of neutrino oscillations [172], in experiments like T2K [173] 
or NOz/A [174] as well as in the (far) future neutrino factories and/or beta 
beam experiments [175]. 
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6.1.2 Neutrinoless double-beta decay and the Majo- 
rana phases 

Theories with neutrino mass generation of the see-saw type predict the mas- 
sive neutrinos Uj to be Major ana particles. Determining the nature of mas- 
sive neutrinos is one of the most formidable and pressing problems in today's 
neutrino physics (see, e.g., [172, 176, 177]). Even if neutrinos are proven to 
be Majorana fermions, getting information about the Majorana CP-violating 
phases in U, $i and $2 [cf. Eq. (A. 5)], will be very difficult. Moreover, the 
oscillations of flavor neutrinos, i^a — J^a' a-nd — > Pa', cc, a' = e,ii,T, are 
insensitive to the Majorana phases [168, 178]. On the other hand, they can 
affect significantly the predictions for the rates of lepton-fiavor-violating de- 
cays 11 ^ e + J, T ^ 11 + J, etc. in a large class of supersymmetric theories 
with type-I see-saw mechanism (see, e.g., [179-181]). 

In the case of three-neutrino mixing under discussion, there are, in prin- 
ciple, three independent CP violation rephasing invariants. The first is Jcp, 
the Dirac one, associated with the Dirac phase S, which we discussed in the 
previous subsection. The existence of two additional invariants, 5*1 and S2, 
is related to the two Majorana CP-violating phases in U. The invariants Si 
and S'2 can be chosen as [182] 



The rephasing invariants associated with the Majorana phases are not 
uniquely determined. Instead of the flavor r involved in both and 5*2, 
one could have chosen e or Note also that CP violation due to the Majo- 
rana phases imply that both = Im {U*-JJ-,-2) 7^ and Re {U*-JJt2) ^ 0, and 
similarly for S'2. 

The only feasible experiments which, at present, have the potential of es- 
tablishing the Majorana nature of light neutrinos and of providing informa- 
tion on the Majorana CP- violating phases in U are the experiments searching 
for neutrinoless double-beta (Oz//3/3) decay, (A, Z) {A, Z + 2) + 2e~. The 
Oi//3/3-decay effective Majorana mass, as defined in Eq. (1.18), which 

contains all the dependence of the 0z//9/9-decay amplitude on the neutrino 
mixing parameters, is given by the following expressions for a normal hier- 
archical (NH, mi -C 1712 -C ma), inverted hierarchical (IH, mi -C m2 -C ma, 
but the U matrix changes, as explained in Appendix A) and quasi-degenerate 
(QD, mi^2,3 ^ m>0.1 eV) neutrino mass spectra: 



Si = Im {U*,Ur2) , S2 = Im {U:^Ur3) . 



(6.1) 



|(m)| ~ -^A^sin^ ^126^*2 + yA;;4;;sin2 ^i3e-'*\ (NH), (6.2) 
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|(m)| ~ m\cos^9i2^^' +e'^^sm^9i2\, (QD). (6.4) 

Obviously, 1(171)1 depends strongly on the Majorana phases: the CP- 
conserving values of ^2 — ^1 = 0,±n [183], for instance, determine the 
range of possible values of | (m) | in the case of IH and QD spectra, while 
the CP-conserving values of $2 = 0, ±7r, can be important in the case of NH 
spectrum. 

The planned Oz//3/3-decay experiments of the next generation such as 
GERDA [49], MAJORANA [50] or CUORE [51] are aiming to probe the QD 
and IH ranges of |(m)| (see, e.g., [176, 177]). If the Oi//3/3-decay is observed in 
these experiments, the measurement of the Oi//3/3-decay half-life might then 
allow to obtain constraints on the combination of Majorana phases $2 — ^1- 

6.2 Dirac phase leptogenesis 

In this section we would like to investigate in detail a scenario of leptogenesis 
where the CP violation necessary for leptogenesis is exclusively provided by 
the Dirac phase S, hence the name (5-leptogenesis. Following closely [164] we 
start analysing the case where the heavy neutrino mass spectrum is hierar- 
chical, i.e. Ml ^ M2 ^ Afs, which, in the fully flavored regime, docs not 
necessarily lead to a A'^i-dominated scenario. Therefore, we shall not only 
consider the decay of the lightest RH neutrino, but that of all three RH neu- 
trinos. Then, we will assume a quasi-degenerate mass spectrum for the heavy 
neutrinos. Mi ~ M2 ~ M3, and perform the same analysis. 

6.2.1 The hierarchical hmit 

Throughout this section, we will always assume that the fully flavored 
regime holds. However, the condition of validity of the fully flavored regime, 
Eq. (4.6), is qualitative. One should not expect it to be exact. For example, 
this simple condition neglects the effect of AL = 1 scatterings and of refrac- 
tive effects, the first contributing with inverse decays to preserve the quantum 
state coherence, the second, conversely, in projecting it on the flavor basis. 
Both of them can be as large as the effect from inverse decays. Moreover, in 
a rigorous quantum kinetic description, it is likely that other subtle effects 
will contribute to the determination of the exact value of Mi below which 
the fully flavored regime can be assumed. Only making a precise study of the 
transition region between the unflavored and the fully flavored regimes solv- 
ing the full system of equations, which includes the density matrix equation 
(4.29), will allow to derive a better condition of validity. 
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In the plots showing the lower bound on Mi, we will then distinguish 
four regions. All plots will be cut at Mi = 10^^ GeV, since above this value, 
according to the condition (4.3), the unflavored regime is recovered and the 
asymmetry production has to switch off, since only the phases in the fl matrix 
can make the total CP asymmetries Si non-zero. On the other hand, when 
the condition (4.6) is satisfied, one can expect the fully fiavored regime to 
hold. There is an intermediate regime where a transition between the fully 
flavored regime and the unflavored regime takes place. This regime will be 
indicated in all plots with a squared region. This signals that, even though 
we still show the results obtained in the fully flavored regime, important 
corrections are expected, especially when Mi gets close to ~ 10^^ GeV. Since 
this region describes a transition towards the unflavored regime, where the 
asymmetry production has to switch off, these corrections are expected to 
reduce the final asymmetry, making more stringent the lower bounds shown 
in the plots. Furthermore, since large corrections to the condition (4.6) cannot 
be excluded, we will also indicate, with a hatched region, the area where the 
condition (4.6) holds but a very conservative condition. 



10^^ GeV 



does not. In this region some corrections to the results presented cannot be 
excluded, but the fully fiavored regime should represent a good approxima- 
tion. 

We anticipate that, in the A^i-dominated scenario, successful leptogenesis 
always requires Mi > 10^ GeV, where the two-flavor regime applies. There- 
fore, considering that we are assuming Si = Sir + £i,e+/i = 0, Eq. (3.29) can 
be specialized into 

showing that, in order to have a non- vanishing flnal asymmetry it is necessary 
to have that ^ -Pi°e+^- Indeed, the efficiency factors crucially depend on 
Kia- Useful analytical expressions for the flavored efficiency factors, both for 
a vanishing and a thermal initial A'^i-abundance can be found in Section 3.4. 
As for the fiavored decay parameters in the usual orthogonal parametrization, 
Eq. (2.18), they were given in Eq. (3.23). 

The parameter tiq, defined in Eq. (3.31) is convenient to describe the be- 
havior of the CP asymmetry eia- Using Eq. (3.28) for real fl, one obtains [129] 

= - > ^hi ^11 Im t/aft U^il . (6.7) 
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Taking a — t and specifying the matrix elements Uaj from Eq. (A. 5), one 
has 

ri, = -!^[A2 + A3 + ^23], (6.8) 

mi 

where, in the case of normal hierarchy for the light neutrinos^, 

V^ima (ma - mi) 

Ai2 = 2 ^^11 ^^21 Ini[(Si2 S23 " C12 C23 Sl3 6 ) 

X (Ci2 S23 + Si2 C23 S,, e"^^) C"^ (*-*^)] , 

Vmima (m3 - mi) i5x ic&n 

^13 = 2 Sill i^3i C23 ci3 Im[(si2 S23 - C12 C23 si3 e )e2 ij , 

^atm 

^23 = _ V"^2?7^3 ^02lf^3lC23Cl3lni[(Ci2S23 + Si2C23Si3e'^)e^*2] 

^atm 

In the case of 5-leptogenesis ($i = $2 = 0), these expressions further spe- 
cialize into 

Vmim2 (ma -mi) 

>ll2 = 2 ^ ^11 ' ^21 S23 C23 A , 

mf 



''atm 

_ Vmim3 (m3-mi) o ^2 ^ . 

/I13 — 2 i ^11 i ^31 C23 C12 Ci3 l\ , 



^atm 

_ V"^2 m3 (m3 - m2) 2 a 

^23 — 2 21 i '31 C23 -S12 Ci3 l\ , 

^atm 

where we defined A = sin ^13 sin 5 . 

It is now instructive to make some general considerations. Looking at the 
expression Eq. (6.6), one can see that in order for the final B~L asymmetry to 
be non-zero, two conditions have to be simultaneously satisfied : Sir 7^ and 

^ '^i,e+/<- These two conditions are a specialization of two of Sakharov's 
necessary conditions for baryogenesis to the case of 5-leptogenesis. 

The first condition is the requirement to have CP violation and, as one 
could expect, from the expressions found for the terms A^j, one can have 
Sir + only if A 7^ 0. 

The second condition is a specialization of the condition of departure 
from thermal equilibrium in quite a non-trivial way. Indeed, in the case of 
5-leptogenesis, in a fully out-of-equilibrium situation where only decays are 
active, no final asymmetry is generated since e\ = 0, implying that there is an 
equal number of decays into leptons and antileptons. However, the presence 
of inverse processes can remove this balance, yielding a different washout rate 



^For an inverted hierarchy, the elements Aij must be computed using a modified U 
matrix, as explained at the end of Appendix A. 
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for the T asymmetry and for the e + // asymmetry, so that, if Kit 7^ -f^i,e+M) 
one has a dynamical net lepton number generation. From the flavored decay 
parameters Eq. (3.23), it can be seen that this is possible independently of 
the value of the Dirac phase, which is therefore directly responsible only for 
CP violation and not for lepton number violation, exactly as in neutrino 
mixing, where lepton number is conserved. 

It should also be noticed that the ei^s are expressed through quanti- 
ties lm[Uah U*i] that are invariant under a change of the PMNS matrix 
parametrization [161,182]. Actually, the CP asymmetry itself is an invari- 
ant quantity [184]. Therefore, the final B — L asymmetry depends only on 
physical quantities, as it should be. 

Maximizing the asymmetry over all involved parameters at fixed Mi and 
Ki and imposing r^^^^ > 77^^^ [cf. Eqs. (2.15) and (1.6)], a lower bound 
on Ml is obtained, which can be conveniently expressed as in Eq. (3.32). 
Notice that ri^- cx A, implying N^^-l oc A as well. Therefore, the maximum 
asymmetry is obtained for \S\ — 7r/2 and S13 = 0.20. 

As we explained in Section 5.2, accounting for flavor effects it is not clear 
whether the contribution from the lightest RH neutrino, iVi, is the only one 
that matters. Since wc want our analysis to be as general as possible, we will 
always estimate the asymmetry produced by the two heavier neutrinos N2 
and A^3. 

The calculation of the contribution to the asymmetry from A^2-decays pro- 
ceeds in an analogous way. It can be again calculated in the two-flavor regime, 
since, in the HL, successful leptogenesis always implies M2 > 10^ GeV. There- 
fore, one can write an expression similar to Eq. (6.6) for the contribution to 
the final asymmetry from A^2-decays, 

^i-LU,-(4.-4e+M)^2.. (6.9) 

The difference is now in the calculation of the efficiency factors as they are 
suppressed by the washout from the A^i inverse processes. In the HL this 
additional washout factorizes and [94, 123, 149] 

i4a - i^{K2a) exp (y-^ Ki^ , (6.10) 

where = -P2°a-^2- The tree- level projectors P^a can be readily evaluated 
using the general expression (3.22). 

The calculation of the contribution to the final asymmetry from A^s-decays 
proceeds in a similar way and analogous expressions hold. The only non- 
trivial difference is that now, in the calculation of the efficiency factors, one 
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has also to include the washout from N2 inverse processes, so that 



4a - «(^3a) exp 



3 TT 

-— (Xl« + X2a) 



(6.11) 



Notice that in the calculation of K^2a 



'^3^) we are not including a possible 



effect where part of the asymmetry in the flavor a = e + /x produced in N2- 
or iVa-decays is orthogonal to Ni inverse decays [127, 150] and is not washed 
out. This washout avoidance does not apply to the asymmetry in the r-flavor. 
Since in all cases we shall consider, a r-dominated scenario will be realized, 
we shall simply neglect this effect. For a short discussion about these effects, 
see Section 5.2. 

Let us now calculate the final asymmetry in some interesting cases. 



Q — i?i3 

The first case we consider is = i?i3 [cf. Eq. (2.21)], implying A12 = A23 = 
in Eq. (6.8). It is easy to check from Eq. (3.26) that £2t = 0, and therefore 
there is no asymmetry production from A^2-decays even if M2 < 10^^ GeV. 
On the other hand, one obtains 

2 A/mi rris (rris — mi) / ~ 2 a / \ 

'^ar = - o ~ ^31 V 1 - UJu C12 C23 Ci3 A , (6.12) 

essentially the same expression as for ri^ but with mi replaced by fhs. There- 
fore, for M3 < 10^^ GeV, one has to worry about a potential non-negligible 
contribution from A^a-decays. However, when the washout from A^i and N2 
inverse processes is taken into account [cf. Eq. (6.11)], we always find that 
the contribution from A^a-decays is negligible and the A^i-dominated scenario 
holds. 

The results are shown in Fig. 6.1 for S13 = 0.20, S — — 7r/2 and mi/matm = 
0.1, a choice of values that approximately maximizes the final asymmetry 
and yields the lower bound Mf™(i^i). In the left panel we show the tree- 
level projectors P^^ and the normalized CP asymmetries ria [cf. Eq. (3.31)]. 
It can be seen how for Ki ^ 10 one has P^^ ~ ^le+n — ^I'^i while for 
Ki ~ 10 one has <^ ^te+n' ■'■^ central panel ^1 and the ^i^'s are 
plotted [cf. Eq. (3.33)], and one can see how for ~ 10 a r-dominance is 
realized. Finally, in the right panel, we show M™™(ii'i) and we compare it 
with the lower bound from an unfiavored calculation where uj^^ is taken purely 
imaginary [93]. One can see how, at Ki ^ 10, the asymmetry production 
rapidly dies, so that ^1 — > and M"""(i^ri) — > 00. Notice that we have plotted 
the lower bound both for thermal and vanishing initial A^i-abundances. 
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Figure 6.1: Dependence of different quantities on Ki for mi/matm = 0.1, 
Si3 = 0.2, 6 = — 7r/2 and real Q = i?i3 with uj^i < 0. Left paneh projectors 
and normalized CP asymmetries ri^; central panel: ^iq, and as defined 
in Eq. (3.33) for thermal (thin) and vanishing (thick) initial A^i-abundance; 
right panel: lower bound on Mi for thermal (thin solid) and vanishing (thick 
solid) initial iVi-abundance compared with the unflavored result (dash-dotted 
line) obtained for complex Q = R^. In the squared region the condition (4.6) 
is not satisfied, and in the hatched region even the more conservative con- 
dition (6.5) is not satisfied. The dotted lines (thick for vanishing and thin 
for thermal initial A^i-abundance) correspond still to a real Q = Ris but this 
time 5 — while the only non-vanishing low-energy phase is the Majorana 
phase $1 = — 7r/2. 
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We also indicated X*, defined as that value of Ki such that for Ki > 
the dependence on the initial conditions can be neglected and the strong 
washout regime holds. One can notice that the intermediate regime between 
a fully flavored regime and the unfiavored regime, the squared quite 
extended. In this regime corrections to the results we are showing, obtained 
in the fully flavored regime, are expected in a way that the unfiavored regime 
should be recovered for Mi — * 10"^^ GeV. In this limit the asymmetry pro- 
duction has to switch off and thus one expects that the lower bound on Mi 
becomes more restrictive and eventually, for Mi — > 10^^ GeV, the allowed re- 
gion has to close up. One can see that there is essentially no allowed region in 
the strong washout regime outside the squared area. The hatched area, where 
corrections cannot be excluded within current theoretical uncertainties, cuts 
away almost completely any allowed region even in the weak washout regime. 

In conclusion, the allowed region where one can safely rely on the fully fla- 
vored regime according to current calculations is very restricted and confined 
only to a small region in the weak washout regime. 

M3 > 10^^ GeV 

The second case we consider is the hmit M3 » 10^^ GeV. This corresponds to 
an effective 2-RH neutrino model, where the third RH neutrino is decoupled. 
In this limit one has necessarily mi <^ mgoi, implying ms ~ matm, and the fl 
matrix takes the special form [160, 185, 186] 

/ M 

n = ^r^T4 -1^31 . (6.13) 
V ^31 ^/^^ 0/ 

Notice that this form of Q corresponds to set 0032 — 1 and 0021 = 1 in 
Eq. (2.19). In the expression (6.8) for ri^, one has now A12 — Ais — 0, 
and hence 

nrc^'^.f^ (1-^) n3iA/r^4ci3.i2A. (6.14) 

mi V ^atm V m^tmj ^ 

If M2 > 10^^ GeV, there is no contribution from the next-to-lightest RH 
neutrino decays, since these occur in the unfiavored regime where £2 — 0. On 
the other hand, if M2 < 10^^ QgV, then one has to worry about a (flavored) 
asymmetry generation from A^2-decays. A calculation of £20 shows that the 
first term in Eq. (3.26) vanishes while the second term gives 

r2r--l^.f^(l-^) 1^31 0^4^12 A. (6.15) 
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Figure 6.2: Same quantities as in Fig. 6.1 but for the case M3 ^ 10^'' GeV, 
corresponding to the special form of Q in Eq. (6.13). Here we are moreover 
assuming M2 > 10^^ GeV and a normal hierarchy for the hght neutrinos. The 
lower bound M^^'^{Ki) is obtained for lo^i > and 5 — 7r/2. 



This is an example of how the second term in Eq. (3.26) is not suppressed in 
the HL like the first term. However, like for the contribution from A^3-decays 
in the case Q = R13, when the washout from A^i inverse processes is taken 
into account, one finds N^_l\j^^ <C Nq_]^\j^^ and a A^i-dominated scenario 
is realized. 

Notice that there is a strong dependence whether one assumes a nor- 
mal or an inverted hierarchy. For normal hierarchy the results are shown in 
Fig. 6.2 for u-^i > and S = 7r/2. For inverted hierarchy the asymmetry is so 
suppressed that there is no allowed region. This means that for any choice 
of the parameters one always obtains M™'" > 10^^ GeV. 

Notice that some results for 5-leptogenesis in the particular case where 
M3 » 10^^ GeV have been recently presented in [163] for vanishing initial 
iVi-abundance. For example the authors obtain a lower bound sin ^13 > 0.09 
imposing the existence of an allowed region for Mi < 5 x 10^^ GeV, while 
we would obtain sin ^13 > 0.05. Different reasons might explain this discrep- 
ancy. First, we are using a (~ 30%) more conservative lower bound on Mi 
[cf. Eq. (2.53)]. Second, we do not include the effects of AL = 1 scatterings 
on the efficiency factors, which might produce some difference when Ki^ ~ 1, 
even though in the strong washout regime this effect should be negligible. Fi- 
nally, another likely minor source of difference is that we are not accounting 
for the effect of spectator processes encoded in the matrix [cf. Eq. (3.15)] 
that relates the B/3 — La asymmetries to the La asymmetries [127]. However, 



98 Chap. 6: Leptogenesis from low-energy CP-violating phases 



notice that we do not want here to emphasize too much a precise value of this 
lower bound on sin ^13, since we believe that it is anyway affected by much 
larger theoretical uncertainties on the validity of the fully flavored regime. 



= R12 

The third case we consider is = R12 [cf. Eq. (2.20)]. This time one has Ai^ = 
A23 = in Eq. (6.8). In the case of normal hierarchy, the CP asymmetry, 
compared to the case Q — R13, is suppressed by a factor (msoi/^riatm)^^^, 
while it is essentially the same for inverted hierarchy. The projectors present 
very similar features to the case Q = R13. One can also again calculate, for 
M2 < lO"*^^ GeV, the contribution from A^2-decays to the final asymmetry, and 
one finds again that the first term in Eq. (3.26) vanishes, while the second 
produces a term oc Mi, so that 

2 y/mi m2 (ma - mi) r - 

= o ~ a;2i A / 1 - S23 C23 A . (6.16) 

3 m2 matm ^ 

When the efficiency factors are taken into account, one finds that only in 
the case of normal hierarchy the contribution to the final asymmetry from 
Ai'2-decays can be comparable to the one from A^i-decays. However, in this 
case both productions are suppressed and there is no allowed region in the 
end. In the case of inverted hierarchy, the contribution from A^2-decays is 
always negligible compared to the one from A'^i-decays. Notice, moreover, 
that e^a = for r2 = R12, and hence there is no contribution from iVs-decays. 

In conclusion, for Q, — R12, the lower bound on Mi for normal hierarchy 
is much more restrictive than in the case ft — R12, while it is very similar for 
inverted hierarchy. A production from the two heavier RH neutrinos can be 
neglected and the A^i-dominated scenario always applies when the asymmetry 
is maximized. 



n^R. 



23 



The last interesting case is ^2 = R23 [cf. Eq. (2.22)]. From Eq. (3.28) one can 
easily check that Sia — 0. One can also check that, contrarily to the case 
Q = Ri2: the second term in Eq. (3.26) vanishes while the first term does not 
and yields 

£2r _ V"^2 m3 (m3 - m2) 



= i(M2) = V 1 - ^12 C23 CIS A . (6. 17) 

Notice that this time e2T oc M2 and actually, more generally, one can see 
that this expression is obtained from Eq. (6.8) for rir in the case Q = it!i3. 
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Figure 6.3: Dependence of different quantities on K2 for mi = 0, S13 = 0.2, 
6 = n/2 and real fl = R23 with UJ22 > 0. Left panel: projectors quan- 
tities rgQ,; central panel: ^2a ^-nd ^2 for thermal (thin) and vanishing (thick) 
initial A^i-abundance; right panel: lower bound on M2 for thermal (thin sohd) 
and vanishing (thick solid) abundance compared with the unflavored result 
(dash-dotted line) as obtained in [93]. 



just with the replacement (Mi, mi) (M2,m2)- At the same time, one has 
Ki — mi/m*, so that the washout from A^i inverse processes vanishes for 
mi — > 0. For M3 < 10^^ GeV, one has to worry about a possible contribu- 
tion to the asymmetry from A^s-decays. A straightforward calculation shows 
that £30; = (2/3)e2Q, and therefore an asymmetry is produced at T ~ M3. 
However, we verified once more that the washout from N2 inverse processes 
is always strong enough for the contribution to the final asymmetry from 
A^3-decays to be negligible. 

In complete analogy with the unflavored case [93] , one has that the lower 
bound Mi™(Xi) is replaced by a lower bound M^''^{K2) obtained for a;32 > 
and shown in the right panel of Fig. 6.3. One can see that also in this case, 
within the validity of the condition (4.6), the allowed region is constrained 
to a small portion falling in the weak washout regime. Assuming the very 
conservative condition of validity for the fully flavored regime outside the 
squared and hatched regions, there is essentially no allowed region even in 
the weak washout regime. 

One can wonder whether there is some choice of Q beyond the special cases 
we analyzed where the final asymmetry is much higher and the lower bound 
on Ml much more relaxed, especially in the strong washout regime. We have 
checked different intermediate cases, and we can exclude such a possibility. 
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Therefore, the lower bound shown in Fig. 6.1 has to be considered, with good 
approximation, the lowest bound for any choice of real ft. 

Another legitimate doubt is whether going beyond the approximations we 
made the lower bound in Fig. 6.1 can be considerably relaxed. The inclusion 
of non-resonant AL = 2 or AL = 1 scatterings is not expected to produce 
large corrections. Recently the effect of the off-diagonal terms in the matrix 

[cf. Eq. (3.15)] has been considered, but it has been shown that it does 
not produce any relevant change in the final asymmetry [130] . 

Relevant corrections, as already pointed out, can only come from a full 
quantTim kinetic treatment, which should describe accurately the transition 
between the unflavored regime and the fully flavored regime. 

The same kind of considerations holds for the iV2-dominated scenario, 
realized for Q = R23. As soon as Q deviates from R23, the washout from A^i 
inverse processes comes into play suppressing the final asymmetry and, at 
the same time, S2t gets also suppressed. Therefore, the lower bound on M2 
is necessarily obtained for Q = R23 in complete analogy with the unfiavored 
treatment [93]. 

In conclusion, 5-leptogenesis in the HL is severely constrained, confirming 
the conclusions of [123] and [122]. In particular, imposing independence from 
the initial conditions, then not even a marginally allowed region seems to sur- 
vive. Notice moreover that all plots have been obtained for S13 = 0.2, the cur- 
rent 3(7 upper limit. Assuming that for values of Mi above the condition (4.6) 
the unflavored regime is quickly recovered and therefore that the asymmetry 
production quickly switches off, then a one-order-of-magnitude improvement 
of the upper limit on sin^^i3 would essentially rule out 5-leptogenesis in the 
HL, even the marginally allowed regions falling in the weak washout regime. 

In the next section we consider the effect of close heavy neutrino masses 
in enhancing the CP asymmetries and relaxing the lower bounds on Mi , M2 
as well as the related one on Tret- But before concluding this section, we 
want to mention that in the more general case of real Q with non-vanishing 
Majorana phases, an upper bound mi < 0.1 eV has been obtained in the fully 
ffavored regime [162] and for a hierarchical heavy neutrino mass spectrum. 
This bound clearly applies also to 5-leptogenesis, but in this case, considering 
the results we have obtained and the expected quantum kinetic corrections 
to the fully flavored regime, the issue is actually whether an allowed region 
exists at all in the HL, even for mi = 0. Therefore, we do not even try to 
place an upper bound on mi in the HL. In the next section we show that an 
upper bound on mi from successful 5-leptogenesis actually holds even in the 
resonant limit, where the CP asymmetries are maximally enhanced. 
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6.2.2 The degenerate limit 

We would like to show that going beyond the HL, the lower bound on Mi 
(or on M2) can be considerably relaxed. Nevertheless, we shall see that some 
interesting constraints on the involved parameters still apply. For simplicity, 
we assume from the beginning that Mi (or M2) <^ 10^ GeV, so that the 
three-flavor regime applies, where the muon- Yukawa interactions are also 
faster than inverse decays. Therefore, when we show the flavor index a, we 
mean a — e, ii,t. This assumption simplifies the calculation, since we do not 
have to describe a transition between the two and the three-flavor regime, 
and because we can completely neglect the effect envisaged in [127, 150] where 
part of the asymmetry produced from iV2-decays is not touched by A^i inverse 
decays. 

In order to go beyond the HL, it is convenient to use the quantity 5ji 
deflned in Eq. (5.5). We are interested in the degenerate limit (DL), where 
at least one 6ji is small enough that both the asymmetry production from 
iVj^j-decays and the washout from the corresponding inverse processes can 
be approximately treated as if they occurred at the same temperature, so 
that they can be simply added up. The DL is a good approximation for 
\Sji\<OM [94]. 

If i,j 7^ 3 and Mi ^ M2 <^ M3, then one has a partial DL, where 
the efficiency factors can be approximated as in Eq. (5.2) for thermal ini- 
tial A'"i-abundance. In all the cases we shall consider, we will always have 
+ Kja ^ 1, so that the strong washout regime applies, and there is no 
need to consider the case of vanishing initial iVi-abundance. 

Another possibility is to have a partial DL with i,j 7^ 1, so that 
Ml -C M2 — M3. In this case one has to take into account the washout from 
the lightest RH neutrino, implying 

<a ^ 4a ^ + Kja) exp ("y^ia) • (6-18) 

Finally, in the full DL, Mi ~ M2 — M3, and the efficiency factor is given by 
Eq. (5.1). 

Let us now calculate the flavored CP asymmetries. In the case of real fl, 
which implies real {h^h)ij = {h)h)ji, the expression (5.4) becomes 

^ ^ ^^"^ ^"'^ ^^'^^^ 

We can again express the neutrino Yukawa coupling matrix through the 
orthogonal representation. This time the presence of the factor does not 
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allow to remove the sum on j', as it has been possible in the HL in order to 
derive Eq. (3.28). However, considering the same special cases as in the HL, 
only one term j ^ i survives, and we can write 

2e{Mi) ^ mn^/mPmi n \n n doit \tt* tt i (ann\ 

Ha ^ 2^ ^ni ^nj Phi ^Ij-^li ^hj] l^[Uah ^al] ■ (6.20) 

■' n,h<l 

The same expression holds for eja simply exchanging the i and j indexes. 
We can always choose j > i such that Mj > Mi. In all the particular cases 
we shall consider, we will have that Eka — for k ^ i,j, and moreover the 
following simplifications apply: 

rrin flni ^nj = {mg - nip) Qji fljj (6.21) 

n 

Vm^i [Qhi ^ij - ^li ^hj] = y/mgrup , (6.22) 

h<l 

with q > p. Except for the case M3 » 10^^ GeV, we will always have q — j 
and p — i. The final B — L asymmetry can then be expressed as 

^h-L - + ^J") '^ai^i^ + Kja, K^a) = ^(mi, Viji, 613, S) A , 

O Oji 

(6.23) 

where 

g{mi, ^13, S) = — — % Jl - 

'"'atm ^ 

X E ^aiK^c. + i^,a, K,^) 5^^!^^^ (6.24) 
a 

and where K.l^{Kia + Kja,Kka) = i^la = i^^ja given by one of the three 
expressions Eq. (5.1), Eq. (5.2) or Eq. (6.18) according to the particular 
case. 

It is interesting to notice that because of the unitarity of U, the sum over i 
of the flavored decay parameters, Eq. (5.3), tends to m/m^ in the degenerate 
limit for the light neutrinos (m is the common mass scale) independently of 
the flavor. Therefore, the sum over the flavor in Eq. (6.24) tends to vanish. 
This will contribute, as we shall see, to place a stringent upper bound on the 
absolute neutrino mass scale in the full DL. 

It is also worthwhile to notice that the sign of A cannot be pre- 
dicted from the sign of the observed baryon asymmetry, since the sign of 
g{mi,i}ji,9is,5) depends on the sign of which is undetermined. Notice 
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also that Im[C/*^ C/q,;]/A does not depend on A, but nevertheless there is a 
dependence of g{'mi,Q,ji,9i3,S) on S and on ^13 coming from the tree-level 
projectors in the sum Ki^ + Kj^- However, in any case, for A — one 
has that g{'mi,flji,9i3,6) A — >• 0, since the final asymmetry has to vanish 
when sin ^13 or sin 5 vanishes. 

The function \g{mi,D,ji,9is,5)\ can be maximized over fi^j. Indeed for 
mi = 0, since k < 1 and Ki + Kj < K-^tm, one has g^rrii = 0, Ki, 613, S) < 4. 
Increasing mi there is a suppression due to the fact that Ki > mi/m^, and 
5'max(^i, 6'i3, 5) decreases monotonically. Therefore, for any mi, there is a 
lower bound on Mi given by 

Ml > Mr (mi, ^13, 5) ^ 1^^' ^ . (6.25) 

gm3x{mi,9i3,5) |A| 

The CP asymmetries, and consequently the final baryon asymmetry, arc max- 
imally enhanced in the extreme case of resonant leptogenesis [76,117,118], 
when the heavy neutrino mass degeneracy is comparable to the decay widths. 
An approximate resonance condition was given in Eq. (5.7), where the un- 
certainty parameter d — 1 -j- 10 was introduced because of a discrepancy in 
the literature [85,117]. When the resonance condition is satisfied, one has 
El = 1/d in the unfiavored case with maximal phase. This can be taken as a 
conservative limit that implies, maximizing over 6, a lower bound 

sin ^13 > sin ^r = ^-^^ n ■ • .M ■ 6.26 

Osph maxs[gniax{mi, 9^"^, d) sm d)\ 

Let us now specialize the expressions for the four special cases we have 
already analyzed in the HL. 



M3 > 10^^ GeV 

Remember that in this case one has {h^h)3j = 0, implying e^a = 0, a conse- 
quence of the fact that the heaviest RH neutrino decouples. Moreover, one has 
mi <^ msoi, so that terms oc mi can be neglected, ~ matm and 1712 — mgoi 
for normal hierarchy or m2 — matm\/^ — '"^soi/''^atm fo'^ inverted hierarchy. 
Therefore, there is actually no dependence on mi in g'(mi, Jljj, ^13, S), which 
is given by Eq. (6.24) with = (1, 2) and (p, q) = (2, 3), 

^^^^^''-'^ " — lui^. — ['-^J V^^^^v^-^^^ 

^.(Xi. + i^2.)^^^^i%^. (6.27) 



X 

a 
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Figure 6.4: Case M3 ^ 10^"^ GeV for normal hierarchy in the DL. Plot of the 
function \g{Ki,9i2,5)\ in the limit A — > 0. The maximum gives the lower 
bound on Mi [cf. Eq. (6.28)] and on sin^ig [cf. Eq. (6.30)]. 

In the case of normal hierarchy, |5'(f^2i; ^13, ^)| slightly decreases when A 
increases, and so the maximum is found for A = and the dependence on 
^^13 and 6 disappears. Replacing the dependence on ^221 with a dependence 
on Ki, we have plotted \g{Ki,A = 0)| in Fig. 6.4 for central values of mgoi 
and matm- Including the errors, one finds g^max — 0.160 ± 0.005. 

The (3(7) lower bound on Mi for normal hierarchy, from the general ex- 
pression (6.25), is then given by 

Ml > 0.9 X 10^° GeV ^ . (6.28) 

In the case of inverted hierarchy, the situation is somehow opposite, 
since for ^13 = the electron flavor contribution vanishes in Eq. (6.27) 
and there is an exact cancellation between the r and /i contributions. Con- 
sequently, the asymmetry increases for increasing values of ^13 and the 
maximum is found for sin ^13 = 0.2 and S ~ n/A. In this case one has 
that maxgjg ,5[(7jnax("^i = 0,^13, 5) A] ~ (9 ± 2) x 10~^, which, plugged into 
Eq. (6.25), gives at 3a 



Ml > 6 X 10^^ GeV 5; 



(6.29) 



6.2 Dirac phase leptogenesis 



105 



It should be remembered that these conditions have been obtained in the 
three-flavor regime and in the DL, i.e. they are valid for Mi 2 ^ 10^ GeV. 
This implies that ^21 ^ 10~^ |A| for normal hierarchy and ^21 ^ 10~^ for 
inverted hierarchy. 

Analogously, the general expression (6.26) gives the following (Scr) lower 
bounds on sin ^13 for normal and inverted hierarchy, respectively: 

sin ^13 > 3.3 X 10-^ d and sin ^13 > 0.06 d. (6.30) 

Q — i?i3 

In this particular case, the next-to- lightest RH neutrino is decoupled from 
the other two heavy neutrinos, which implies that £20 — for any a and 
that Eia does not depend on M2; in particular, it does not get enhanced if 

^21 — ^ 0. Therefore, one has necessarily to consider 531 < 0.01, implying a 
full DL with all three RH neutrino masses quasi-degenerate. The function 
g{mi, flji, 6*13, 6) is now obtained from Eq. (6.24) with j = q = 3 and i = p — 
1 and = K,{Kia + K^a + K^a), or explicitly 

g{mi,n3i,0i3,5) = — — ^ flsiJl-nii 

X J2 <Ki^ + + i^3a) ^^^^^^^^ . (6.31) 

a 

It is interesting to notice that in this case an e-dominance is realized. More- 
over, one has that the dependence of \g{mi, Vt^i, 613, 6) \ on 9is and 5 is slight, 
and the maximum is for A = and mi = 0. We find 5'max(0) = 0.24±0.01 for 
normal hierarchy and g'max(O) = (3.1 ±0.2) x 10~^ for inverted hierarchy. The 
lower bound on Mi in Eq. (6.25) yields then at 3a for normal and inverted 
hierarchy, respectively. 

Ml > 5.5 X 10^ GeV ^ and Mi > 5 x lO" GeV ^ , (6.32) 

while in the case of resonant leptogenesis Eq. (6.26) yields 

sin ^13 > 2.3 X 10"^ d and sin ^13 > 1.5 x 10"^ d . (6.33) 

Increasing mi, the value of gmax{mi) decreases, and the lower bound on sin ^13 
in resonant leptogenesis becomes more and more restrictive. This dependence 
is shown in Fig. 6.5 both for normal (left panel) and inverted (right panel) 
hierarchies and for d = 1 (solid line) and d = 10 (short-dashed line). In- 
terestingly, imposing the experimental {3a) upper limit sin ^13 < 0.20, one 
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Figure 6.5: Case Q — Ris in the full DL. Lower bound on sin versus uii 
obtained in resonant leptogenesis for 0? = 1 (solid line) and d — 10 (short- 
dashed line). Values sin ^13 > 0.20 are excluded at 3a by current experimental 
data. 



obtain the upper bound mi < 0.2-0.4 eV, depending on the value of d. This 
upper bound will become more stringent if no signal for a non-vanishing 
mixing angle ^13 is seen in future neutrino oscillation experiments, since the 
experimental limit on sin ^13 will then go down. Assuming no discovery, the 
most stringent experimental upper limit is expected to be reached in neu- 
trino factories, where one could obtain sin ^13 < 10^"^ [187]. This asymptot- 
ical upper limit is also shown in Fig. 6.5 and would imply an upper bound 
mi < 0.05-0.1 eV for normal hierarchy and mi < 0.03-0.08 eV for inverted 
hierarchy. Therefore, an interesting interplay between two measurable quan- 
tities is realized, making 5-leptogenesis falsifiable independently of the RH 
neutrino mass spectrum. 

In the more conservative case of normal hierarchy (see left panel of 
Fig. 6.5), a good approximation for the upper bound on mi {d — 1) is given 
by the fit 

mi < 0.6 (sin ^13 - 2.3 x 10"^) eV . (6.34) 

It is interesting that this upper bound holds in the extreme case of resonant 
leptogenesis and therefore holds for any RH neutrino spectrum. However, we 
still have to verify if it holds also for a different choice of fl. 
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Figure 6.6: Case Q = Ru in the partial DL. Lower bound on sin versus 
mi obtained in resonant leptogenesis. Same conventions as in the previous 
figure. 



^ — R12 

The situation for Q = R12 is quite different compared to the previous two 
cases. One has now i = p = 1 and j = q = 2, and it is possible to have both 
a partial DL with Mi ~ M2 < M3 < lO^"^ GeV and a Ml DL. In the first 
case, the general expression (6.24) becomes 

/ n a x\ - 2 /sTatrn (^1 + ^2) (^2 - mi) f — 

X Y: + K.o) ^^^^^^^ . (6.35) 

a 

This time the contribution from the electron fiavor vanishes. Furthermore, 
for normal hierarchy, there is an almost perfect cancellation between the 
and the r contributions. In the left panel of Fig. 6.6, we show the lower bound 
on sin ^13 versus mi, and one can see that it is more restrictive than in the 
previous case, Q, = Ri^. In particular, imposing sin ^13 < 0.2, one obtains 
now a much more stringent upper bound mi < 0.06 eV. On the other hand, 
for inverted hierarchy, the cancellation between the /x and the r fiavors does 
not occur, and one has a lower bound on sin6'i3 for mi <^ 0.01 eV that is 
similar to what has been obtained in the case Q = it!i3 (see the right panel of 
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Figure 6.7: Case Q — R12 in the full DL. Lower bound on sin ^13 versus 
nil obtained in resonant leptogenesis. Same conventions as in the previous 
figures. 



Fig. 6.6). However, here there is no flavor cancellation for increasing values of 
mi because Kia + K2a does not tend to a common value like Kia + K2a + K^^- 
Therefore, one can see in Fig. 6.6 that this time the upper bound on mi is 
much looser, both compared to normal hierarchy and compared to = R13. 

In the full DL the flavor cancellation at large mi occurs. The results 
are shown in Fig. 6.7. One can notice that the upper bound on mi is very 
restrictive for normal hierarchy, and one has a situation similar to the case 
fl = R12 for inverted hierarchy. 



^ — R23 

When fl = R23, the lightest RH neutrino decouples and eia = indepen- 
dently of Ml . Therefore, there is no contribution to the flnal asymmetry from 
A^i-decays. On the other hand, £20 and £30, do not vanish and hence there 
is a contribution from the decays of the two heavier RH neutrinos. Still Ni 
inverse processes have to be taken into account since they contribute to the 
washout. There are two different possibilities. 

In the full DL the washout from A^i inverse decays just cumulates with 
the washout from the two heavier. Therefore, using expression (6.24) with 
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Figure 6.8: Case Q — R23 in the full DL. Lower bound on sin ^13 versus 
mi obtained in resonant leptogenesis. Same conventions as in the previous 
figures. 



i = p = 2 and j = ? = 3 and fi;^ = n{Kia + i^2Q + K^a), one obtains 

^(mi, 1^32, ^13, 5) = — — ^ 032^1-1) 

(6.36) 



1? 
'32 



atm 



In Fig. 6.8 we show the lower bound on sin ^13 versus mi for successful res- 
onant leptogenesis. This time there is a bigger suppression than in the case 
Vt = i?i3, especially for inverted hierarchy. 
In the case Mi <^ M2 ^ M3, one has 

^ o /) jcN 2 i^atm (^2 + K3) {tus - ms) .Jm^m^ r 
g[mi,ilz2,0i3,6) = — — 2 "32Vl- 

Ko Ki m,r^.-_ v 



^32 



X ^ K{K2a + K^a) exp f - ^ Ki^ j 



37r^ \\m[Ul^U^,] 



.(6.37) 



The lower bound on sin ^^13 versus mi is shown in Fig. 6.9 for normal hierarchy. 
One notices that the upper bound on mi is slightly less stringent than in the 
full DL. For inverted hierarchy the asymmetry production is so suppressed 
that there is no allowed region. 



110 Chap. 6: Leptogenesis from low-energy CP- violating phases 



10-^ 
10-2 
10"^ 

as 

■i 10-" 

10-^ 
10-^ 
10-^ 

0.001 0.01 0.1 

m, [eV] 

Figure 6.9: Case Q — R23 in the partial DL. Lower bound on sin ^13 versus 
mi obtained in resonant leptogenesis. Same conventions as in the previous 
figures. 

We can conclude this section noticing that our results show that 6- 
Icptogenesis can be falsified. In the case of normal hierarchy, the current 
upper limit sin ^13 < 0.2 implies mi < 0.1-0.3 eV, while, in future, a poten- 
tial upper limit sin ^13 < 10~^ would imply mi < 0.01-0.1 eV, with a more 
precise determination depending on the possibility of improving the current 
estimation of the parameter d in resonant leptogenesis. 

6.3 Leptogenesis from the Majorana phases 

It is apparent from the elements Aij in Eq. (6.8) that the Majorana phases 
in U also contribute to the CP violation necessary for leptogenesis. Actually, 
they can also play the role of unique source of CP violation, as illustrated in 
Fig. 3.8 for Q = i?i3 real and non-zero Majorana phase $1 ($1 = tt/2). One 
can even say quite generally that the Majorana phases can be more easily 
responsible for enough CP violation than the Dirac phase, simply because 
they are not associated with the small mixing angle ^^13. To illustrate this, 
we have plotted in the right panel of Fig. 6.1 with dotted lines the case of 
$1 = — 7r/2 and 5 = 0, compared with 5 = — 7r/2 and $1 = for sin ^13 = 0.2, 
i.e. the maximal allowed value. One notices that the lower bounds on Mi for 
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the case of non-vanishing Majorana phase is about a factor of 2 lower than 
in the case of non- vanishing Dirac phase, even for the maximal value of ^13. 

We do not aim here at making a thorough study of the role of the Ma- 
jorana phases as the sole source of CP violation for leptogenesis. We think 
that the case of the Dirac phase, which we analysed in detail in the previ- 
ous section, represents a more conservative situation, and the prospects for 
a measurement seem to be more encouraging. However, for completeness, we 
would like to report some results obtained in [163] and [166] concerning ex- 
clusively the Majorana phase and where only the limit of hierarchical heavy 
neutrinos and a vanishing initial A^i-abundance were considered. 

When a fully hierarchical light neutrino spectrum (mi <S mgoi) is con- 
sidered, the authors in [163] obtain two main results: for a normal hierarchy 
(real Vt matrix), they find the lower bound Mi > 3.6 x 10^° GeV for successful 
leptogenesis, whereas for inverted hierarchy (purely imaginary they 
find Ml > 5.3 x 10^° GeV. It should be noted that for inverted hierarchy and 
real Q, there is no allowed region below 10^^ GeV. 

In [166] the study was extended to arbitrary mi, and some new effects 
were found. In particular, even in the case of inverted hierarchy and real Q, 
an allowed range was found, with the lower bound Mi > 3 x 10^*^ GcV. In the 
case of normal hierarchy, no such relaxation occurs and the bounds quoted 
above still hold. 

6.4 Discussion 

We have discussed situations where the "observable" CP- violating phases, S, 
$1 and $2, act as the only source of CP violation responsible for the matter- 
antimatter asymmetry of the Universe. Such possibilities, especially for the 
Dirac phase, which can be realistically discovered in the future, represent 
by themselves a strong motivation. We may indeed soon be in the situation 
to probe the second of Sakharov's necessary conditions for baryogenesis (see 
Section 1.1). 

As we have seen, successful leptogenesis from low-energy phases is only 
marginally possible in the HL, Mi <C M2 <C M3, and with dependence on the 
initial conditions. This is especially true when the only source of CP violation 
is the Dirac phase S. We have also argued that a definite conclusion on the 
existence of such a marginally allowed region requires a quantum kinetic 
treatment, which is expected to shrink the already quite restricted allowed 
region. 

Therefore, 5-leptogenesis and more generally leptogenesis from low-energy 
phases motivate models with quasi-degenerate RH neutrino masses, the ex- 
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treme limit being resonant leptogenesis. Even in this extreme limit, imposing 
successful 5-leptogenesis we could derive interesting conditions on quantities 
accessible in low-energy neutrino experiment: sin ^13, the absolute neutrino 
mass scale, normal or inverted hierarchy, the Dirac phase itself. An interest- 
ing aspect of 5-leptogenesis is then that it is falsifiable independently of the 
heavy neutrino mass spectrum. 

There are however some objections to the scenario of leptogenesis exclu- 
sively from low-energy phases. At the moment, it still lacks a strong theo- 
retical motivation. In [162] a model where such a situation naturally arises 
was shortly discussed. There, it was said that the simplest way of restricting 
the number of CP- violating phases is through the assumption that CP is a 
good symmetry of the Lagrangian, only broken by the vacuum. For example, 
one can add to the standard type-1 see-saw framework three Higgs doublets, 
together with a Z3 symmetry under which the left-handed fermion doublets 
ipLj transform as ■0Lj — e'^'^^^'^i/'Li and the Higgs doublets as — > e"^^'^-'/^^^, 
while all other fields transform trivially. It can be readily shown that there 
is a region of parameters where the vacuum violates CP through complex 
vacuum expectation values (0°) = fjc'^'. Due to the Z3 restrictions on the 
Yukawa coupfings, the combination h'^h is real, thus implying a real Vt matrix, 
but keeping U complex. Of course, even though such a model might work, it 
is not the simplest and most economical one. 

There has been a recent claim that sequential dominance models [188] 
(see [189] for a more complete discussion) could represent a theoretical frame- 
work for leptogenesis from low-energy phases. In [184] it was shown that these 
models correspond to have an Q matrix that slightly deviates from the unit 
matrix or from all the other five that can be obtained from the unit matrix 
exchanging rows or columns. However, it has been noticed in [88, 93] that in 
the limit Im[r2] 0, the total CP asymmetries ei do not necessarily vanish. 
Writing Vt^- = |r2^^|e''^'^ , the correct condition to enforce £j — > is to take 
the limit — > 0. This is a more demanding limit than Im[fi] — > 0, and it is 
not currently motivated by generic sequential dominance models. This limit 
is not motivated either by radiative leptogenesis [144, 145] within the context 
of minimal flavor violation [190], as recently considered in [165, 191, 192]. It 
must however be said that the limit Im[fl] — > 0, when assuming a vanishing 
initial abundance of RH neutrinos, as done in the works cited above, can 
effectively mimic the condition ^ because the efficiency factor Ki — > 0, 
when Im[n] and Rc[0] 0, implying Ki ^ 0. 

Another possible objection to leptogenesis from low-energy phases is that 
it cannot be distinguished from the general scenario where both high- and 
low-energy phases are present. In particular, the Dirac phase will likely give 
in this case only a subdominant contribution, since its effect is always sup- 
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pressed by the small ^13 mixing angle. Following this approach, one can even 
say that the baryon asymmetry produced through leptogenesis is not sensi- 
tive to the phases in U [193], in the sense that the baryon asymmetry can 
be accounted for with the phases of U having any value. Conversely, if the 
phases in U are measured, the baryon asymmetry is still not constrained. 
The hope is then that the theoretical framework supporting leptogenesis 
from low-energy phases has some other testable predictions. An experimen- 
tal support for this model would then represent a support for leptogenesis 
from low-energy phases. 

However, even if a model supports leptogenesis from low-energy phases, 
how can one know if the source of CP violation comes from the Dirac phase or 
from the Majorana phases? Actually, it was noticed that the contribution to 
the final asymmetry from Majorana phases is in general dominant compared 
to the one coming from the Dirac phase (see right panel of Fig. 6.1). Thus, it 
will be only possible to tell if 5-leptogenesis really occurs once the Majorana 
phases are constrained from neutrinoless double-beta decay experiments to 
give small contributions. 

One can even imagine a situation where there is an exact cancellation 
between the Dirac and Majorana contributions. It would be however strange 
to think that nature disposes a sufficient source of CP violation, sets up a 
second source that exactly cancels the first one, and the observed asymmetry 
is explained by yet a third one, e.g. the phases in Jl. 
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Chapter 7 
Conclusion 



The amount of baryonic matter in the Universe which we infer, for instance, 
from the CMB temperature anisotropics represents one of the most impor- 
tant puzzles of modern cosmology. In order to explain this number, one needs 
a hary agenesis mechanism which generates dynamically a small baryon asym- 
metry in the early Universe. We discussed in the introduction that a solution 
to this problem necessarily leads to physics beyond the SM. It is very excit- 
ing that this puzzle of cosmology may actually be related to the existence of 
tiny but non-zero neutrino masses, which are now established. As a matter of 
fact, a simple extension of the SM naturally leads to small neutrino masses 
via the see-saw mechanism, and its cosmological consequence is leptogenesis, 
which elegantly yields the required baryon asymmetry. 

In the present thesis, we have thoroughly discussed the mechanism of 
leptogenesis, where a lepton asymmetry is produced by the decays of heavy 
right-handed (RH) neutrinos and then transferred to a baryon asymmetry by 
the non-perturbative sphaleron processes. Let us now summarize the main 
points discussed throughout the thesis. 

First, we have presented the unflavored treatment of leptogenesis, where 
the leptons produced in the decays of the heavy neutrinos have no flavor 
structure. It was shown that in the hierarchical limit for the heavy neutrino 
mass spectrum, Mi ^ M2 ^ M3, it is typically enough to consider only 
the decay of the lightest RH neutrino (iVi-dominated scenario), with a re- 
duction of the number of parameters relevant for the computation. In this 
minimal scenario, which we refer to as the "vanilla" scenario, the contribution 
from the lightest RH neutrino washes out all previous asymmetry and, for 
the theoretically favored values of the effective neutrino mass [cf. Eq. (1.26)] 
"^soi ^ "^1 ^ ^atm, where mgoi and matm stand for the solar and atmospheric 
neutrino mass scales, respectively, the strong washout regime is obtained, 
with no dependence of the flnal asymmetry on the initial conditions. Ad- 
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ditionally, the strong washout regime imphes that the simple picture with 
only decays and inverse decays, i.e. neglecting all scattering processes, and 
without including thermal corrections yields a very good estimation of the 
final asymmetry. 

Interestingly, in the A^i-dominated scenario, general constraints on a few 
parameters of the model can be derived. For successful leptogenesis, the mass 
of the lightest RH neutrino, Mi, cannot be smaller than about 4x 10^ GeV for 
the strong washout to be obtained [81, 96]. This lower bound leads to a related 
lower bound on the initial temperature of leptogenesis. Tin > 1.5 x 10^ GeV, 
which can be identified with a lower bound on the reheat temperature Tj-^h 
within infiation. Such a high value of the reheat temperature may be in 
confiict with locally supersymmetric theories due to an overproduction of 
gravitinos. 

The baryon asymmetry produced through leptogenesis is also sensitive to 
the absolute neutrino mass scale. In the context of vanilla leptogenesis, the 
stringent upper bound mi < 0.12 eV was obtained [75, 88, 99]. 

Then, wc introduced the "flavored" picture of leptogenesis, which has 
been understood only recently to be the correct one for a large fraction of 
the parameter space (roughly when Mi < 10^^ GeV) [119, 120]. Flavor effects 
introduce a dependence of the final asymmetry on essentially all parameters 
of the model. On the one hand, this makes the computation more involved, 
but, on the other hand, some interesting parameters, such as the CP-violating 
phases in the PMNS matrix, become accessible. About fiavor effects, one can 
say in general: 

• In most of the parameter space they lead to modifications of the pre- 
dictions by a factor 2-3 compared to the unfiavored analysis, due to a 
reduction of the washout by this factor. Incidentally, the region of inde- 
pendence from the initial conditions shrinks by the same amount [123]. 

• Large modifications are possible in two cases: i) when the high-energy 
phases in the ft matrix are zero or close to zero; ii) when a one-fiavor 
dominance is obtained. 

• The usually quoted values of the lower bounds on Mi and Treh (see 
above) do not change when fiavor effects are included [123] . 

Next, we re-analysed the conditions on the temperature T and on Mi 
for flavor effects to be important. We found that there should be a region 
in the parameter space where classical Boltzmann equations ( "fully flavored 
regime") are not enough to describe the generation of asymmetry [133]. Cor- 
relations in flavor space and partial losses of coherence might indeed be rele- 
vant there, so that a quantum kinetic equation in the form of a density matrix 
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equation should be used. Interestingly, the region concerned is exactly the 
one where the upper bound on the absolute neutrino mass scale seems to be 
evaded. Therefore, at the moment it is not clear whether the upper bound 
on mi from successful vanilla leptogenesis quoted above disappears, is simply 
relaxed, or still holds. 

Then, we went beyond the minimal picture where only the lightest RH 
neutrino is considered. In the quasi-degenerate limit, Mi ^ M2 — M3, the 
contributions from all RH neutrinos have to be included. Due to the enhance- 
ment of the CP asymmetry in this limit [87], the lower bounds on Mi and 
Treh can be lowered down to the TeV scale in the extreme case of resonant 
leptogenesis [76,117]. Accounting for flavor effects, this might even be pos- 
sible without having all Yukawa couplings unnaturally small [132, 146], even 
though it remains to be proven that the condition of validity of the fully 
flavored equations is satisfied in this case. 

Actually, the production of asymmetry from the heavier RH neutrinos 
N2 and A^3 is not only important when considering quasi-degenerate heavy 
neutrinos. Even for hierarchical heavy neutrinos, it was noticed in [93] that 
a particular choice of the Vl matrix leads to production of asymmetry by 
N2 instead of Ni. This implies that the lower bound on Mi does not apply 
anymore and is replaced by a lower bound on M2, which still yields a lower 
bound on Treh [93] . When flavor effects are included, the contributions from 
N2 and A^3 are potentially more important, due to the reduced washout from 
A^i [149,150]. Specifically, the domain of applicability of the A^2-dominated 
scenario is expected to be enlarged, and N2, might be important as well, even 
though it seems comphcated to avoid the washout both from A^i and N2. 

The last possibihty of going beyond the typical scenario that we dis- 
cussed is when one element of the Vt matrix, namely ^22, has a non-trivial 
value different from or 1. It turns out that making IO22I large opens the 
possibility of relaxing the lower bound on Mi by 3 orders of magnitude even 
for Ml <C M2 [152]. Moreover, with moderate values of 1^222 1 and moder- 
ate degeneracies Mi ~ M2, the upper bound on the neutrino mass can be 
evaded [98]. 

Finally, we studied the special case where the source of CP violation 
required for leptogenesis stems exclusively from the phases in the PMNS ma- 
trix. We focused on the Dirac phase, for which the prospects of measurement 
are the most promising. We found that for hierarchical heavy neutrino masses 
leptogenesis from the Dirac phase (5-lcptogcncsis) is only marginally allowed 
and in the weak washout regime. For quasi-degenerate heavy neutrinos, the 
strong washout is recovered, and we could even derive a upper bound on 
mi dependent on ^13 for successful resonant leptogenesis, which is the most 
favorable case one can imagine. Roughly speaking, for the 3(7 upper limit 
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sin ^13 = 0.2, resonant 5-leptogenesis only works for mi < 0.2 eV [164]. If 
the experimental upper limit on sin ^13 decreases in the future, so does the 
upper bound on mi for successful 5-leptogenesis. 

The see-saw mechanism and leptogenesis have very appealing features. 
However, in their vanilla form, they seem very difficult to prove or disprove. 
The scale of the heavy neutrinos necessary for leptogenesis, as well as to 
have a "natural" see-saw mechanism, is too large to be accessible at future 
colliders. The only possibility to produce heavy neutrinos at colliders would 
be to have at the same time TeV masses and large Yukawa couplings. Even 
though such a possibility relies on a cancellation mechanism rather than a 
see-saw mechanism [143, 148] , the question whether it is imaginable to have 
signals from heavy RH neutrinos at colliders is very exciting by itself [143, 
147, 148] . The inclusion of successful leptogenesis in the picture certainly 
deserves investigation. 

Another way of probing directly leptogenesis would be to measure a pri- 
mordial lepton asymmetry of the order of the baryon asymmetry. The relic 
neutrino background would carry such an information. Although the sole 
detection of relic neutrinos is in principle possible [194-196], measuring an 
asymmetry of order 10~^° in it seems hopeless. 

If a direct test of leptogenesis seems to be out of reach, then one has to 
wait for an accumulation of indirect hints in favor of this scenario. 

First, the observation of neutrinoless double-beta decay would establish 
the Majorana nature of light neutrinos, hence supporting the see-saw mecha- 
nism and leptogenesis. The next-generation experiments such as GERDA [49] 
or CUORE [51] are likely to observe a signal if the light neutrino mass hi- 
erarchy is quasi-degenerate or inverted. Furthermore, the discovery of the 
Majorana nature of neutrinos would immediately tell that lepton number is 
violated, hence verifying the first Sakharov's condition, thanks to the pres- 
ence of sphalerons. 

Second, the discovery of CP violation in the neutrino sector in future 
long-basehne neutrino experiments such as T2K [173] or NOi/A [174] will 
certainly strengthen the case for leptogenesis. It will indeed tell that the 
second Sakharov's condition is verified. Moreover, the source of CP violation 
in neutrino mixing might be sufficient to explain the origin of the matter- 
antimatter asymmetry of the Universe without resorting to the unobservable 
high-energy phases in the Q matrix. 

Conversely, it is interesting to ask oneself if there are ways to disprove 
leptogenesis. A non-observation of neutrinoless double-beta decay in the next- 
generation experiments would not disprove by itself the see-saw mechanism 
and leptogenesis. The mass hierarchy might simply be normal, hence out 
of reach of the next-generation experiments. However, if neutrino oscillation 
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experiments determine the mass hierarchy to be inverted [197,198], but no 
signal is found in neutrinoless double beta decay experiments, then one would 
conclude that the see-saw mechanism is not at the origin of the neutrino 
masses, and that leptogenesis is not responsible for the baryon asymmetry 
of the Universe. 

As concerns CP violation, the non-observation of leptonic CP violation 
in future neutrino experiments will not weaken the case for leptogenesis in a 
significant way. Instead, it would mean that the Dirac phase and/or the angle 
^13 are too small to give an noticeable contribution to the final asymmetry. 
The source of CP violation responsible for leptogenesis can be still given by 
the remaining 5 CP- violating phases present in the model. 

The absolute neutrino mass scale may also provide a way to contrain 
significantly leptogenesis. As a matter of fact, we have seen that the latter 
scenario leads to an upper bound on the absolute neutrino mass scale of 
about 0.1 eV, which is however subject to modification due to fiavor effects. 
Assuming that the upper bound remains below roughly 0.3 eV, implying 
< 0.3 eV, a signal in the future KATRIN experiment [41], which claims 
a discovery potential down to about 0.35 cV, would severely constrain lep- 
togenesis in its minimal version. On the other hand, other versions, such as 
the one with quasi-degenerate heavy neutrinos, would still remain viable. 

Even though the production of heavy neutrinos at future colliders seems 
to be unlikely, the LHC experiment at CERN will still have an indirect impact 
on leptogenesis. Indeed, if supersymmetry is discovered and the parameters 
measured are consistent with successful electroweak baryogenesis (e.g. light 
stop and light Higgs [199,200]), then the case for leptogenesis will become 
weaker. On the other hand, if electroweak baryogenesis is ruled out at the 
LHC, then leptogenesis, as one of the remaining possibilities to explain the 
baryon asymmetry of the Universe, will become stronger. 

In conclusion, since leptogenesis is unavoidable when considering the see- 
saw mechanism for the generation of neutrino masses, it provides a very 
elegant explanation to one of the outstanding problems of modern cosmology, 
the origin of the matter-antimatter asymmetry of the Universe. Even though 
a direct test is challenging, there is no doubt that in the next few years more 
experimental evidence will become available to weaken or strengthen the case 
for leptogenesis. 
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Appendix A 

Neutrino mixing parameters 



The currently existing data from neutrino oscillation experiments can be well 
described by the Lagrangian, 

9 



C = --^hY^LWi, - -{uLYm^UL - iRmeh + h.c. , (A.l) 

which includes the weak charged current interaction in the lepton sector, a 
Majorana mass term for neutrinos, and a Dirac mass term for charged lep- 
tons. When diagonalizing the neutrino and charged lepton mass matrices via 
rriu = U*m'^^^^Ul and rrii = Vem'^^^^Uj, one obtains the lepton mixing ma- 
trix, known also as the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing 
matrix [17-19], in the weak charged lepton current 

U = UjU,. (A. 2) 

The charged current interaction can then be written in the mass basis as 

-^h^rU^ii^LiW^. (A.3) 

Note that it is conventional to choose the basis where the mass matrix for 
the charged leptons is diagonal, in which case U = U^, i.e. the PMNS matrix 
is the matrix that diagonalizes the neutrino mass matrix. Neutrinos with 
a given flavor are then related to neutrinos with a given mass through the 
PMNS matrix, 1 



= XI ' « = e, T . (A.4) 



^If one wants the relation between the flavor eigenstates and the mass eigenstates, it 



is given by = X]j=i a. = e,fi,T. At first sight, this may seem surprising, 

but it is simply due to the fact that, by convention, field operators create antiparticles, or 
annihilate particles! [23] 
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Note also that the number of neutrino flavor states with mass below is 
known to be three from the precise measurement of the Z- width at LEP I [5]. 

The standard parametrization of the PMNS matrix was given in Eq. (1.8). 
Carrying out the matrix product, one obtains 

C/ = y xdiag(e'^,e^^,l), (A.5) 

(Cl2 Ci3 "S12C13 SisC " 

"-512 C23 ~ C12 >S23 >Si3 e"' C12 C23 — S12 S23 <Si3 e"^ S23 Ci3 
S12 S23 - C12 C23 si3 e"^ -C12 S23 - S12 C23 si3 e"^ C23C13 

We shall use in all calculations Ou — n/6 and 623 — 7r/4, compatible with 
the results from neutrino oscillation experiments discussed in the introduc- 
tion, Eqs. (1.13) and (1.15). For the third angle, we shall use the 3a range 
Si3 = 0-0.2 obtained from Eq. (1.16). 

Concerning the mass spectrum of light neutrinos, we shall use the con- 
vention that nil < ^i2 < iTT-s, whatever the hierarchy is. This means that 
when a normal hierarchy is considered, one has mg — m| = Am^^m ^^<i 
m| — mf = AmgQj [cf. Eqs. (1.12) and (1.14)], whereas for an inverted hier- 
archy, one has ml — m| = ■^^soi ^2 — nT'i — ■^^atm- Defining the two 
convenient quantities 

matm = ^Aml^ + Aml, = (0.052 ± 0.002) eV, (A.6) 

and 

msoi = ^Am^^i = (0.0089 ± 0.0002) eV, (A.7) 

one has that the light neutrino spectrum is quasi-degenerate when 
nil ^ "^atrnj while for nii <^ mgoi it is fully hierarchical. 

Within the convention for the light neutrino masses we are using, it must 
be pointed out that the case of inverted hierarchy is obtained by performing 
a cyclic permutation of the columns in the PMNS matrix Eq. (A.5), such 
that the i-th column becomes the (i -|- l)-th [164]. 

Finally, let us note that we shall always neglect the effect of the running 
of neutrino parameters from high energy to low energy [90,91]. 
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The see-saw mechanism with 
three RH neutrinos 



The (type-I) see-saw mechanism is based on the following extension of the 



where a = e, jj,, r, three new fields Nfn, i = 1,2, 3/ a Majorana mass matrix 
Mm and a Yukawa coupling matrix h have been introduced. We chose here 
the basis where the charged-lepton Yukawa matrix and the Majorana mass 
matrix are diagonal. The superscript c denotes charge conjugation, defined 
as (j)'^ = C(f)^ , where the charge conjugation matrix C satisfies 



The subscripts L and R denote the left-handed and right-handed chiral pro- 
jections, respectively: Pl^r = (1± 75)72. Finally, the SU{2)l x U{1)y charge 
assignments are as follows: 



It should be stressed that the new fields A^^, often called right-handed (RH) 
neutrinos, are singlets under SU{2)l x U{1)y- 

^In principle, it would be possible to consider the addition of only two RH neutrinos, 
since only two light neutrinos are known to be massive. But, in the following, we want to 
be slightly more general and therefore keep three RH neutrinos. 



SM: 




c-\c = -7J, = -c, = c-\ 



(B.2) 
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In Eq. (B.l) the Yukawa-type term is simply the analog of the other mass 
terms for fermions in the SM. Concerning the Majorana mass term, it is 
new and unique, but in full generality there is no reason why it should be 
absent. This particularity for neutrinos is due to the fact that they are the 
only neutral fermions in the SM. 

After spontaneous symmetry breaking, a Dirac mass term mo = hv, is 
generated by the vacuum expectation value of the Higgs field, v = 174 GeV. 
Using the identity V^Nji — {Nii)'^{ulY, one can then rewrite the mass term 
in a more compact form: 

The sec-saw mechanism then assumes that the entries of the Majorana 
mass matrix matrix are much larger than all Dirac matrix elements, i.e. 
Mm 3> m-Q. Under this assumption, the mass matrix in Eq. (B.3) can be 
block-diagonalized with a 6 x 6 unitary matrix V: 



where 



and 



nil, = rnji^—m^ , (B.5) 
Mm 



(B.7) 



The K matrix only makes sure that the mass eigenvalues are positive. 

The matrix in Eq. (B.5) corresponds to the light neutrino mass matrix, 
which has naturally suppressed entries due to the heavy scale in the denomi- 
nator. This is the reason why this mechanism is called the see-saw mechanism. 
As explained in Appendix A, the matrix that diagonalizes the mass matrix 
rriy for the light neutrinos is the PMNS mixing matrix [cf. Eq. (A. 5)], so that 

U^mi,U* — diag(mi, m2, rus), (B.8) 

which are the masses of the three light neutrinos. The corresponding mass 
eigenstates are given by 

l^i^Y. (^^)ia ii^'La - [l^LaT] " ^Da^Mr^ [(A«,)' ' ^^1} ■ (B.9) 
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It can be easily checked that Vi — —i^f, which means that the hght neutrinos 
are Major ana particles. 

On the other hand, the lower right block on the right-hand side of 
Eq. (B.4) was diagonal before the diagonalization and, to leading order, it 
remains diagonal, i.e. M ~ Mm- The entries are Mi < M2 < M3, corre- 
sponding to three heavy neutrinos. The corresponding mass eigenstates are 
given by 

Ni^J2 ^^\'<U [^La + + [{Nr^)' + Nr^ . (B.IO) 

a 

It can be checked again that iV? = Ni, which means that the heavy neutrinos 
are Major ana particles. 

Finally, let us note that the extension of the SM in Eq. (B.l) introduces 
18 new parameters: 6 masses, 6 mixing angles and 6 CP-violating phases. 
The number of parameters in principle accessible in low-energy neutrino ex- 
periments is 9: the masses of the 3 light neutrinos and the 6 parameters (3 
mixing angles and 3 CP- violating phases) in the PMNS matrix [cf. Eq. (A. 5)]. 
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